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Abstract 

Various aspects of the theory of quantum integrable systems are reviewed. 
Basic ideas behind the construction of integrable ultralocal and nonultralocal 
quantum models are explored by exploiting the underlying algebraic structures 
related to the Yang-Baxter equations. Physical meaning of abstract mathe- 
matical notions like universal 7?.-matrix, quantized algebras, Sklyanin algebra, 
braided algebra, Hopf algebra etc. and the role played by them in integrable 
systems are highlighted. Systematic construction of quantum integrable lat- 
tice as well as field models and their exact excitation spectra are presented 
through examples. The coordinate and algebraic formulations of the Bethe 
ansatz are illustrated with comparison, along with the description of nested 
and functional Bethe ansatzes. The techniques for deriving quantum Hamil- 
tonians from the Lax operators are demonstrated on concrete models. The 
exposition of this review is kept in a fairly elementary level with emphasis on 
the physical contents. 
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1 Introduction 



1.1 Background and major concepts 

Search for unifying theories continues in all scientific fields. In recent years we have 
witnessed again the emergence of deep interconnections between seemingly diverse 
subjects like quantum and statistical systems, knot and braid theories, Yang-Baxter 
equation, quantized and braided algebras etc. centered around the theory of inte- 
grable systems in (1 + 1) dimensions. The linkage between integrable systems and 
the conformal field theory, string theory as well as with topological models have also 
been revealed and a wide varieties of review articles are dedicated to these aspects 
d], ^, I, H, 1^, H, 0. The aim of the present review work however, is to focus on the 
subject of the integrable systems itself, concentrating mainly on the role played by the 
recently discovered algebraic structures underlying this theory and their applications 
in constructing and solving physical models. We attempt to present this beautiful 
interplay between the abstract mathematical objects in one hand and the physical 
problems in integrable quantum models in the other, in an elementary way. 

The theory and applications of nonlinear Integrable systems, an immensely im- 
portant subject in mathematical physics of today, has made profound influence in 
many branches of physics as well as in other disciplines. The seed of its development 
can possibly be traced back to 1934, when on an August day of a British engineer- 
historian named Scott Russell had a chance encounter with a strange stable wave in 
the Union canal of Edinburgh p|. Next evidence of such paradoxically stable solu- 
tion was found in the famous computer experiment of Fermi, Ulam and Pasta [^]. 
However only in the mid-sixties such phenomena were understood to be the exact 
localized solutions of integrable equations like Korteweg-de Vries (KdV) equation 

ut{x, t) + Un,xx{,x, t) + t)u{x, t)x = 0, (1.1) 

solvable through inverse scattering method (ISM) |jlO| and were named as Solitons 
[|Tl|]. Subsequently, this was revealed to be an universal feature of all nonlinear 



integrable systems in (1 + 1) dimensions |1T2| , p!3i [14 



Mathematical basis of classical integrable systems was laid down mainly through 
the works of Sofia Kawalewskaya, Fuchs, Painleve, Liuoville and others |T5[. There 
exist in fact various definitions of integrability |T^. We however, for subsequent 



generalization to the quantum models, shall adopt the notion of integrability in the 
Liuoville sense, where integrability means the existence of action-angle variables. 
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That is, if in a Hamiltonian system H[p{x, t), q{x, t)] given by the nonhnear equation 

5H . 5H 

P=-^, 1=^^ (1-2) 
dq op 

it is possible to find a canonical transformation {p{x, t), q{x, t)) (ci(A), b{X, t)) , 
such that the new Hamiltonian becomes dependent only on the action variables, i.e. 
H[a{X)] , then the system may be called completely integrahle. In this case the dynam- 
ical equations d = — ^ = 0, h = ^ = uj, can be trivially solved and as a bonus 
we also get a(A) as the generator of the conserved quantities ln(a(A)) = YfjLo cjX^ , or 
ln(a(A)) = J2'jLoC{-j)^~'' with C(-i-j) = . The number of such independent set of 
conserved quantities {cj} coincides with the degree of freedom of the system, which 
in case of field models with infinite degrees of freedom also becomes infinite. One of 
these conserved quantities may be considered as the Hamiltonian of the system. For 
example, for the KdV equation C(_3) is usually taken as the Hamiltonian. For dis- 
crete systems with finite number of particles, only number of conserved quantities 
(usually first A^) among an infinite set turns out to be independent. Therefore the 
existence of infinite set of conserved quantities itself is not enough, their involution 
(commutativity) must also be established for concluding the complete integrability. 
We shall see later that the validity of the famous Yang-Baxter equation usually 
guarantees such involution, even in the quantum case. 

The major development in the theory of integrable systems has been achieved 



in seventies and eighties with the invention and application of ISM [|T^, |14|, P). The 
important feature of this method is that, instead of attacking the nonlinear equation 
(|1.2|) directly, one constructs corresponding linear scattering problem 

%{x, A) = L{q{x,t),p{x,t), A) T(x, A), (1.3) 

where the Lax operator L{q,p, A) with the fields q,p contains all information about 
the original nonlinear system. Some concrete examples of Lax operators are given 
with the list of integrable models in sect. 2. The parameter A, known as the spectral 
parameter, acts like the momentum or the rapidity of the scattering waves, while 
the original field q becomes the scattering potential. Interestingly, the function a(A) 
related to the transition and 6(A, t) to the reflection coefficient of the scattering 
process are linked with the action and the angle variables, respectively. The aim 
of the ISM is presizely to find the canonical mapping between the original field 
and the action-angle variables and using this to construct the exact solutions for 
the original field. Soliton solution is a special solution, which corresponds to the 
reflectionless (6 = 0) potential and the form of a(A) as a(A) = jE^^ with. 
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zeroes (poles) at discrete points Aj (A*) in the complex spectral parameter plane. 
We shall see in sect. 7 that even in quantum problems some bound states (called 
string solutions) exhibit such solitonic structure. The Hamiltonian is identified from 
the infinite set of conserved quantities {cj}, generated by Ina(A), by knowing the 
Poisson bracket and the dynamical equation. The Lax operator L{q,p, A) therefore 
is the true representative of an integrable system and together with the Poisson 
bracket can generate the whole integrable hierarchy. This fact makes the systematic 
construction of the Lax operators, to which a major portion of the present review is 
devoted, so important. 

A new era started with the formulation of quantum inverse scattering method 



(QISM) by extending the notion of integrability to the quantum domain |16, 17, IS]. 
To highlight various aspects of quantum integrable systems is the main purpose of the 
current review. The action-angle variables in the quantum case become operators 
with generalized commutation relation 

A(A)5(/i) = /(A-^)5(/i)A(A) (1.4) 

For concluding complete integrability of a quantum system, the infinite set of con- 
served quantities, which now become operators must be in involution: Cm] = 0. 
This is guaranteed by [r(A),r(yu)] = , since the transfer matrix r(A) generates 
the conserved operators : Inr(A) = J^jCjX^. The existence of commuting transfer 
matrices in turn is assured by the matrix relation known as the Quantum Yang- 
Baxter equation (QYBE) R{X, n) Ti(A) T2(/i) = T2(/i) Ti(A) R{X,n), Ti = 
T ® /, T2 = / (S> T, involving the quantum i?(A, /i)-matrix and the monodromy 
matrix T(A) = Yl^iLi(X) , where Li{\) is the Lax operator of the correspond- 
ing lattice model at point i and r(A) = trT{X) . The QYBE associated with the 
quantum integrable systems and ensuring the commutativity of action variables for 
different A, exhibits amazingly rich algebraic properties. It represents in fact a set of 
generalized commutation relations like (|1.4| ) between action and angle like variables 
given in the matrix form. 

For the commonly discussed models known as ultralocal models, the Lax operators 
commute at different points: [Lu, L2j] = aand the above QYBE can be obtained 
as a consequence of its local version 

i?(A,/x) L,i(A) L,2(/i) = L,2(/i) L,i(A) i?(A,/x), j = l,...,iV. (1.5) 

On the other hand, for rest of the models constituting another important class, known 



as nonultralocal models |20, 19, 21] we have [Lij,L2j] 7^ , and the QYBE ([TT 
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has to be extended to its braided version (BQYBE) 

-Rl2(A - fl)Z2iLij{X)Z2lL2j{fi) = L2j{^)Zl2Lij{\)Ri2{\ - fi). (1.6) 

with additional matrices Z, Z. For describing the integrabihty one has to complement 
(|1.6|) also with commutation relations between Lu, L2j at different lattice points 
i 7^ j, known as braiding relations. As a result the local BQYBE can be raised to 
the global level giving braided QYBE for the monodromy matrix similar to (|1.6|) . 



1.2 Guide to the contents of the review 



We detail the basic idea behind the quantum and statistical integrable systems in sect. 
2, by introducing quantum i?-matrices, Lax operators, Yang-Baxter equations and 
their extensions [p^, [ 



16, 21 



26 



We also present here a list of well known quantum 
integrable models of both ultralocal and nonultralocal type along with their explicit 
Lax operators, belonging to the trigonometric RtHgW as well as to the rational 
RratW matrices. 

We focus on the algebraic aspects of the quantum integrable systems |2^ 
9i in sect 



27| , p8| , p9| in sect. 3. The significance of Hopf algebras in integrable systems is 
stressed by showing how the local and global QYBE relations are linked together 
through the related coproduct properties and how and why they differ for ultralocal 
and nonultralocal systems, depending on their underlying algebraic structures [pl|j . 
The origin of quantum algebras |3^, 25 1 in integrable theory is traced out and the 

H, its extensions [^, the 



fascinating structures like quadratic Sklyanin algebra 
Yangian algebra p5|, ^ and the braided and quantized braided algebra |^ 



are explored through quantum Yang-Baxter like equations. The notions of quantum 
group 1 139 1 , Faddeev-Reshitikhin-Takhtajan (FRT) algebra [Q, Hecke and Birman- 
Wenzl- Murakami algebra |Q are also introduced in this section. Abstraction of the 
quantum i?-matrix can also be raised to another level of algebraic sophistication 
by defining an Universal 7?--matrix intertwining between two different coproducts in 
an quasi-triangular Hopf algebra The Universal 7^-matrix for the general 

reductive Lie algebras [^, ^ is constructed here with the intention of its subsequent 
application in the next section. 

Such abstract algebraic notions have not only enriched the integrable systems 
alone, but also the algebraic aspects of the quantum integrable theory have made 
profound influence on pure mathematics, by introducing examples of noncocommut- 
ing and quasi-triangular Hopf algebras , giving FRT formulation of dual algebras 
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1^] and introducing various knot and link polynomials through Yang-Baxter algebra 
1- 

For systematic construction of the key objects involved in the quantum integrable 



ultralocal theory, a Yang-Baxterization scheme ^ is formulated in sect. 4 by 
exploiting the universal 7^- matrix. The richness of this abstract object is shown from 
the application point of view, by deriving from it the FRT relations, QYBE, YBE 
and most importantly the spectral parameter dependent quantum i?-matrices and 
the Lax operators. The models thus generated through realizations of the underlying 
quantum algebra and its undeformed limit include the well known sine-Gordon (SG) 
equation, nonlinear Schrodinger (NLS) equation and the XX Z, XXX spin-| models 



as well as the multicomponent Toda field equations |Q. Twisting transformation 



of the Hopf algebra also allows to construct from the TZ the colored and twisted 



generalizations of the models. 

For covering wider class of models an extension of the trigonometric Sklyanin al- 
gebra (ETSA) is considered in sect. 5, which through Yang-Baxterization leads 
to the construction of some integrable ancestor model. Different realizations of the 
underlying ETSA generate variety of quantum integrable models in a systematic 
way. Such models include Liuoville model P6[, massive Thirring model [ffH], rela- 



tivistic Toda chains Ablowitz-Ladik model quantum derivative NLS |^ 



, g-oscillator models etc. apart from the quantum algebra related models found in 
sect. 4. All these descendant models are associated with the same quantum Rtrig 
matrix of trigonometric type. At g ^ 1 limit linked to the Yangian algebra, one can 
construct similarly the related ancestor model for generating its integrable family 
of descendant models associated with the rational i^^at-niatrix solution. An added 
advantage of this approach is that, one can get directly the exact lattice regularized 
Lax operators, which are otherwise difficult to guess from the known field models. 
Such discrete Lax operators, apart from their own interest, yield also the correspond- 
ing field models at vanishing lattice spacing. In this approach one starts from a 
known i?-matrix solutions and constructs subsequently the Lax operator of an an- 
cestor model involving generators of quantized algebra, determined by the QYBE. 
Different reductions of this ancestor model through suitable realizations in physical 



variables like bosons, spin operators or g-oscillators ^J] generate in turn wide range 
of integrable descendant models. The result also brings out the hidden interrelations 
between diverse models and classifies them according to the associated i?-matrix. 
Various realizations of the extended quadratic algebra provides important criteria for 
recognizing integrable nonlinearities. To focus on such application aspects of abstract 
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algebraic objects in quantum integrable theory is the main objective of this section. 

The monodromy matrix is formed by multiplying local Lax operators and the 
conserved quantities are found by expanding the trace of the monodromy matrix in 
the spectral parameter. However, the explicit construction of the conserved operators 
including the Hamiltonian is often a difficult task in quantum models. We take up 
this issue in sect. 6. For the fundamental models [0], which include spin chains 
and Hubberd model, we derive the Hamiltonians in explicit form starting from the 
associated Lax operators. The construction of twisted, colored, impurity as well 
as inhomogeneous models are also discussed here. The generalized spin model and 
the projector methods for constructing the Hamiltonian of nonfundamental models 
[13; HOl; ^-g- bosonic models are described on the example of lattice NLS and the 



related difficulties for the lattice SG model are pointed out. Some special examples of 
nonfundamental models like Toda chains, a simpler lattice NLS allowing explicit 
construction of Hamiltonian are highlighted. 

A problem of significant physical interest in quantum models is to find the energy 
spectrum by solving the eigenvalue problem H \ m >= Em | m >, where H is 
a nonlinear Hamiltonian operator and | m > is the m-particle eigenstate. In case 
of integrable models this problem is usually solvable in an exact and nonperturba- 
tive way. Moreover, one can even find the eigenvalues A(A) for all other conserved 
quantities from the transfer matrix: Inr(A) | m >= InA(A) \ m > by expanding 
InA(A) = Y^'^qCjX^ , similar to the classical case. Interestingly, in this formulation 
the angle variable -B(A) acts like the creation operator and the m-particle state 
I m > can be obtained by applying m times the operators B{Xi),i G [l,m] on the 
pseudovacuum | >. Therefore using the generalized commutation relation (|1.4| ) 
and the properties of the pseudovacuum, one can calculate the required eigenvalues 
A(A). We discuss in sect. 7 the basic concepts of the Bethe ansatz method, a pow- 
erful method for solving the eigenvalue problem exactly. Its algebraic formulation 
p^ , p!6| , |53| , |23| , [1[] is explained on examples of spin models like XXZ and XXX 
chains and bosonic models like the SG and the NLS. We consider in parallel the 
coordinate formulation ||5^, |5^ introduced by Bethe way back in 1931 [0, on 



the same examples, mainly to bring out the linkage between these two approaches. In 



considering more general algebraic structures, the nested Bethe ansatz method [57 



is described for the models like the Hubberd and the vector NLS and the functional 



Bethe ansatz method |]58| , [59| is demonstrated on the examples of relativistic and 
nonrelativistic Toda chains. 

Our next concern in sect. 8 is another major and important class of integrable 
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models, namely quantum nonultralocal models, discussions of which are usually 
avoided in general reviews, due to the nonavailability of their well framed theories 
apart from the studies in Such systems, examples of which are well known 

models like KdV and modified KdV |13|, nonlinear a- model |^0|5 WZWN model 
complex SG model |^ etc. are difficult to tackle, mainly due to the presence 



of derivative 5-function type terms in their commutation relations. This makes 
the transition from local to the global QYBE difficult and the application of the 
standard Yang-Baxter relations related to the quantum algebra fails for such systems. 
However, based on the concept of quantized braided algebra [0, the integrability 
of nonultralocal models may be formulated |2^, much in parallel to the ultralocal 



case, through braided extensions of the QYBE and a trace factorization |]63|, |6J]. The 
applicability of this theory is demonstrated on nonultralocal models like modified 
KdV [Q, WZWN model nonabelian Toda chain |6^, quantum mapping [|67|] , 
Coulomb gas picture of CFT |^ etc. It is remarkable that the supersymmetric |^ 
and anyonic |^ integrable models as well as the systems related to the reflection 
equation [B^ also flt well in this formulation. 



Finally we discuss in sect. 9 some recent developments in this subject and few 
open problems, the starting concepts of which are mostly contained in the present 
review. Our discussions are focused on the multicomponent extension of the ETSA 



3^] with possible applications, Cologero-Sutherland spin model |£T|, its possible R- 
matrix formulation []72| and deformed Yangian symmetry |7^, flnite size correction 



and conformal properties of the integrable systems quantum mKdV and its 
application for an alternative formulation of CFT |75|, reaction-diffusion equation as 
integrable system We discuss further some newly discovered ladder symmetric 
models f?^, |149|| , recent development in elliptic quantum group |j78|, the possibility 
of generating coupled spin chains and extension of Hubberd model through twisting 
transformation |^ and flnally the problem in constructing quantum KdV, nonlinear 
(T-model and complex SG model as quantum nonultralocal systems. 

In sect. 10 we conclude with a few remarks. 

The arrangement of this paper is outlined in the contents. 



2 Quantum integrable systems: basic features 
and properties 



The Lax operator, a key object in classical integrable systems is generalized to the 
quantum case, where it continues to play a major role. Its matrix elements, classically 



12 



being functions of field variables, become now quantum operators. Examples of such 
quantum Lax operators as 2 x 2 matrices representing a variety of integrable models 
are listed in sect. 2.4. The explicit form of the Lax operator determines the associated 
quantum model. We show here how the integrability of such models are ensured by 
certain relations, known as the quantum Yang-Baxter equation (QYBE), satisfied 
by their Lax operators. As a result, using the Lax operator one can construct the 
commuting family of conserved quantities including the Hamiltonian. 

Quantum integrable systems (QIS) can be divided into two broad classes: ul- 
tralocal and nonultralocal depending on an important property of their representative 
Lax operators Lai{X), where a and i denote the auxiliary and the quantum spaces, 
respectively and A is the spectral parameter. For ultralocal QIS, which include well 
known models like Nonlinear Schrodinger equation, sine-Gordon model, Toda chain, 
etc. the Lax operators exhibit a common ultralocality property, i.e. they commute at 
different lattice points i ^ j: [Lij(A), L2j{fi)] = (see fig 2.1). We shall see that this 
fact is actively used in constructing the QYBE. On the other hand, the nonultralo- 
cal models, which include also famous models like quantum KdV , supersymmetric 
models, nonlinear a model, WZWN model etc. are characterized by more general 
property with [Lij(A), L2j{fi)] 7^ 0. This makes the formulation of their integrability 
theory difficult and the associated QYBE needs to be extended along with specifying 
additional braiding relations at different lattice points. We introduce basic equations 
for quantum integrability of both ultralocal and nonultralocal models, which depend 
crucially on the algebraic properties of the associated Lax operators. We should men- 
tion however, that compared to ultralocal models the integrability of nonultralocal 
systems is more involved and not yet well established. 

2.1 Integrability of quantum ultralocal models: Yang- 
Baxter equations 

For better understanding of algebraic structures of Lax operators for ultralocal mod- 
els, let us consider first usual matrices A, B, C, D satisfying the obvious relation 

{A^B^) = {B^A^) or [A,fi2] = 0, (2.7) 

and the standard multiplication rule 

{A®B){C®D) = {AC®BD) (2.8) 

where Ai = A®1, B2 = 1 ® B . If we choose now A = Li{X),B = Li{fi) as 
Lax operators at the same lattice points and check for the relation (|2.71 ), we see 
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immediately its nonvalidity due to the noncommuting nature of the matrix elements 
of the quantum Lax operators. The equation in effect turns into the famous QYBE 
(see fig. 2.2 for i = 3) 

Ri2{X,fi)Li,{X)L2j{fi) = L2j{fi) Lij{X)Ri2{X,fi), (2.9) 

where a matrix R12 with c-number elements appears to compensate for the non- 
commutativity of quantum L. This is the basic reason for the nontrivial algebras 
underlying such integrable systems. There are two main types of i?(A, /i)-matrices 
we will be interested in, which are given through trigonometric or rational functions 
of the spectral parameters X — /j,. The explicit forms of such matrices will be obtained 
below. 

Note however that for the choice A = Li^i{fi),B = Lj(A) at different lattice 
points, the trivial commutation ( p.7|) is indeed satisfied due the ultralocal nature 
of the Lax operators. Let us choose the matrices through Lax operators again at 
different lattice points a.s A = Lj+i(A),i? = Li^i{fi),C = Li{X),D = Li{fi) and 
turn to the multiplication rule ( |2.8| ) . We observe that this rule holds whenever B2 and 
Ci commute and therefore, it remains valid also for our choice due to the ultralocality 
condition on Lax operators resulting (see fig. 2.3) 

iLu+iiX)L2i+M)iLu{X)L2iiix)) = {Lu+l{X)Lu{X)){L2i+l{^l)L2^{^l)) (2.10) 



We show now that combining two commutation relations (|2.9| ) and (|2.10|) of the local 



Lax operators, one can derive the global QYBE essential for proving integrability, 
which in fact is a global property. Start from QYBE ( p. 91) at point i + 1 ( fig. 2.2), 
multiply both sides by (Lij(A)L2i(/u)) and apply in the rhs the same QYBE at i to 
flip the i?-matrix to the extreme right (see fig 2.4). 

We use next the relation ( p.lO|) (fig. 2.3) and the one obtained from it by inter- 
changing 1^2 and A /i to put together the Lax operators acting in the same 
space 1 or 2 in both the Ihs and rhs. As a result we arrive again at the same form 
of QYBE but for a nonlocal object Tf+^^^\X) = Lai+i{X)Lai{X) (see fig. 2.5). 

Repeating this globalization step times we finally obtain the global QYBE (fig 
2.6a) 

i?i2(A,^) Ti(A) T2(/i) = T2(/i) Ti(A) i?i2(A,/i), (2.11) 

for the monodromy matrix T(A) = Y{f Li{X) (fig 2.6b). We shall come back to this 
important transition from local to global QYBE related to the coproduct property 
of the underlying Hopf algebra in sect. 3. 
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Taking now the trace from both the sides of (|2.11| ) and using cychc rotation of 
matrices under trace, we get 



tri2 (Ti(A) TM) = tru {R-'TM Ti(A) r) = tr^ {T^ T,{X)) 

Further with the obvious property tri2{TiT2) = tri(Ti)tr2(T2) we obtain the trace 
factorization tr{T{X)) tr{T{fi)) = tr{T{fi)) tr{T{fi)) or for the transfer matrix r(A) = 
trT(A) the trivial commutation (see fig. 2.7) 

[r(A) ,r(/x)] = 0. (2.12) 

As in the classical case, r(A) generates the conserved quantities usually as In r(A) = 
EJ^qC^A"". Therefore the commuting transfer matrices at different spectral param- 
eters ( p.l2| ) lead to the involution [Cn, Cm] = 0, i.e. to the independence of the 
infinite number (in general) of conserved quantities, which thus proves the exact in- 
tegrability of the quantum system in the Liouville sense. Tracing back this logical 



chain, we see therefore that the integrability is guaranteed by the relation (|2.12|) 



which is originated from the global QYBE ( p.ll|) and that in turn is obtained from 



the local QYBE ( p.9| ) under the ultralocality constraint ( p.lO| ). Consequently, QYBE 
(p.9|) may be considered as the key equation for ultralocal integrable systems repre- 
sented by the Lax operator -Lj(A) and the quantum i?(A, /i)-matrix. For deriving an 
independent equation for the -R(A, /i)-matrix, one can start with the triple product 
Lij(A)L2i(/i)-L3j(z/) and revert the ordering by flipping successively two operators at a 
time using (|2.9|). Clearly this can be done in two different ways :(12) (13) (23) 
and (23) — > (13) (12), which due to associativity should be equivalent and yield 
the equation 

RuiX, yu) RniX, ly) i?23(/i, i^) = i?23(/^, i^) RniX, p) R^iX, //) ., (2.13) 

known as the Yang-Baxter equation (YBE). The equivalence of these two scattering 
paths can also be visualized as in fig. 2.8. 

It should be mentioned that the trigonometric solutions of YBE, we are interested 
in, in the limit A ±oo goes to -Ri2(A, /x) — >■ Rfi2 with q = e**^, the spectral parameter 
free solutions R^ of 

Ri2 Ri3 R23 = R23 Ri3 R12 (2-14) 
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given in the upper and lower triangular form as 



( 



\ 



1 q-q- 
1 



g ) 



R- = P (Rt)-' P 



-{q-q-') 1 



V 1 J 

(2.15) 

with permutation operator P. These solutions are related to the braid group repre- 
sentations as ^ = PR^. 

Note that the spectral independent equation ( p.l4| ) is obtained as the limiting 
case A — oo, /i ^ oo of YBE ( p.l3|) , when the QYBE (|2.9|) also becomes spectral 
parameter free. Taking therefore different sets of limits A— >oo, /i— s>oo, A— > 
— oo, fi — oo, A —>■ oo, fi — oo etc. with A — /i ^ oo we obtain from (|2.9| ) the 
relations 



R+ 



f?+ r+ T- 



^2j LijRgi2^ 



(2.16) 



where L(A)|a->±oo — ^ -^^ are spectral parameter free upper (lower)-triangular ma- 
trices. We shall return to such equations in sect. 3 for describing the underlying 
algebraic structures and in sect. 5 for the construction of Lax operator through 
Yang-Baxterization. The spectral parameter dependent Lax operator L{X) and the 
i?(A, /x)-matrix are sufficient, in principle, for constructing infinite number of con- 
served quantities and also for solving the corresponding eigenvalue problem exactly. 
We demonstrate these features in subsequent sections and show how exploiting the 
underlying algebraic structures the Lax operators can be generated systematically. 



2.1.1 Relation with statistical system, _R- matrix solutions 

The (1-1-1) dimensional quantum systems are intimately connected with 2 dimensional 
classical statistical systems and the notion of integrability is equivalent in both these 
cases. While in the quantum systems the QYBE (|2.9| ) is the central, for statistical 
models similar role is played by the YBE ( p.l3|) . Let us consider a 2-dimensional 
array of x M lattice points connected by the bonds assigned with -|-ve (-ve) signs 
or equivalently, with right, up (left, down) arrows in a random way (see fig. 2.9). 

For this classical statistical system known as vertex model, the partition function 
Z should give the probability of appearance of all possible arrangements on the whole 
structure. However, for analyzing this global property one should know first the local 
ones by finding the probability of occurrence of a particular arrangement at a fixed 
lattice point i. For two allowed signs on each bond, 2^ = 4 possibilities would appear 
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in the horizontal (H) and 4 in the vertical (V) direction. This accounts for the 
4 X 4 = 16 possible arrangements at each lattice point. Setting the corresponding 
Boltzmann weights as the matrix elements of a 4 x 4-matrix, we get the i?^2-matrix 
acting on the product space (product of and V spaces) Vi^V2. Parametrization of 
Boltzmann weights makes a crucial dependence of the R matrix on spectral parameter 
A. For finding out the configuration probability for a string of A^-lattice points, 
the i?(*)-matrices at different sites in a row should be matrix multiplied to give the 
transfer matrix r(a, P) = tr{Y[f -R^*-*)- It is important to note that in this expression 
the horizontal indices fii enter only as summation or repeated indices (see fig.2.9). 
The tr appears also for the same reason from the extreme ends due to the periodic 
boundary condition. As a result the transfer matrix r(a, /3) depends solely on the V- 
space (see fig. 2.9). \^-space is usually called the quantum space due to its significance 
for quantum models, while the if-space is named as the auxiliary space. For finding 
the partition function involving M such strings, one has to repeat the procedure M 
times to give Z = tr(n*^ T) = tr{T^^). 

The YBE ( |2.13| ) restricts the solution of the i?-matrix to integrable models. How- 
ever the R matrix with 16 different Boltzmann weights, representing in general a 
16-vertex model is difficult to solve. Therefore usually some extra symmetry and 
conditions are imposed on the i?-matrix to simplify its solution. For example, we 
shall require the charge conserving symmetry R]^i ^ 0, , only when /c + Z = i + j, 
i.e. the nontrivial elements to exist only when the sum of incoming charges = sum 
of outgoing charges, along with a charge or arrow reversing symmetry (see fig. 2.10). 

Using an overall normalization it leads to a 6-vertex model expressed as 

/ a(A) \ 
6(A) 1 
1 6(A) 

V «(^) J 

where a(A),6(A) are unknown functions of the difference A — /i, and the spectral 
parameters A and fi parametrizing the Boltzmann weights correspond to different 
pairs of spaces. Inserting this form of i?-matrix in the integrability condition YBE 
( p.l3| ), one gets simple equations like a(A — /u)6(A) — a(A)6(A — ji) — 6(/i) = , 
giving a solution through trigonometric functions in spectral parameters as 

a(A) = ^^"^^ + 6(A) = (2.18) 
sin a sin a 

Here a is the parameter fixed for a certain model, while the spectral parameter A is a 

variable one. Apart from ( 2. IS ) the YBE (|2.13|) allows also another simpler solution 



R{X) 



(2.17) 
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of (|2.17| ) given through rational functions in spectral parameters: 



a(A) = ^, 6(A) = -, (2.19) 
a a 

We do not take up more general elliptic i?- matrix solution |2^, ^ , except indicating 
only its relation with the Sklyanin algebra pi[ and the elliptic quantum group [ITS 



in sect. 3 and 9. Apart from the XYZ spin chain and the 8- vertex statistical model, 
there seems to be not many important models associated with the elliptic case and 
the underlying algebraic picture, except some recent developments [|78l, is also not 
much clear. 

2.1.2 Quantum Lax operators of ultralocal models 

To understand the structure of the Lax operators better and to illuminate the rela- 
tionship between the R and L matrices appearing in the QYBE ( |2.9| ), we consider in 
the above picture (fig 2.9) a situation, when the quantum (V) and the auxiliary (H) 
spaces become inequivalent. That is keeping the if-space to be of dimension 4, the 
\^-space is extended to arbitrary dimensions. The i?-matrix then turns into a 2 x 2- 
matrix: -Ri2(A) Li{X) G Vi (S> V, with its elements becoming operators acting in 
a Hilbert space, where to distinguish we denote the arbitrary dimensional quantum 
spaces by bold letters. This L-operator corresponds to the familiar quantum Lax 
operator, which similar to the i?-matrices of the statistical models, constructs mon- 
odromy and transfer matrices for the global description of the associated quantum 
system. 

If we choose for example, the space 3 to be quantum and nonisomorphic to 1 and 
2, then we have -Ri3(A) = i^i(A) and R2s{,lj) = L2{fi), with -Ri2(A — /x) remaining 
the same. As a result, the integrability equation YBE ( p.l3|) for the R matrix of the 



statistical systems turns into the QYBE ( p. 91 ) involving both R and L and describing 
the integrability of the quantum systems (compare fig. 2.2 and 2.8). Conversely, 
when the quantum and the auxiliary spaces become isomorphic, the L , R matrices 
become equivalent and the YBE and QYBE coincide. This happens in the case of 
fundamental models like spin models. 

The i?-matrix appearing in the statistical and the quantum systems however are 
equivalent and given by (|2l7|) , which can also be rewritten as 



with Ui{X) being the Boltzmann weights given by ti;o(A) = a(A) + 6(A), uJsi^X) = 
a(A) — 6(A), uji = UJ2 = 1 and cr* are the Pauli matrices. For inequivalent spaces we 
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may expect the associated Lax operator to be of the form 

L(A) = uJi{X)a' ® S\ a° = h (2.21) 

by replacing (T*'s with some abstract operators S** acting in the quantum space. We 
shall see in sect. 3, that such operators generate certain quadratic algebra intimately 
connected to the well known quantum algebra ||25|. Different realizations of (|2.21| ) 



yield different Lax operators representing a class of integrable quantum models, as will 
be demonstrated in sect. 4. Remarkably therefore, the same trigonometric i?-matrix 
(|2.18|) of the 6- vertex model is associated with a variety of integrable quantum models 



like the XXZ spin chain, sine-Gordon model, Liouville model, massive Thirring 
model, derivative nonlinear Schrodinger equation, Ablowitz-Ladik model etc. This 
intriguing fact will be elaborated further in sect. 4 and 5, while a list of such systems 
with some necessary information is given below. 

Before looking into relevant limits of the parameters, let us clarify certain points 
in view of some confusing remarks that are often encountered in the literature in this 
regard. First we note that, in quantum models the discrete Lax operator Ln(^,r7, A) 
like (|2.31|) contains the spectral parameter ^ = e^"^^, parameter t] and the lattice 



spacing A, but not directly h. On the other hand, the R{^, a)-matrix ( |2.17| ) is in 



dependent of A, but depends on the parameter a = hi], i.e on rj as well as on the 
quantum parameter h. The parameter a induced by the i?- matrix is also present in 
the quantum commutators of the operators involved in L„. The parameter rj repre- 
sents the coupling constant of the quantum models (e.g. anisotropy parameter in the 
XXZ spin chain) and defines at the same time the deforming parameter q = e*'' of 
the underlying quantum algebra. 

Therefore the classical limit h ^ affects the i?-matrix as 

R{X,a) = I + nr{X,r]) + 0{h) (2.22) 

yielding the classical r(A, r/) matrix with A = r/A , which determines the structure 
constants in the Poisson bracket relations among the corresponding Lax operators. 
However, the classical form of the Lax operator stays almost the same, except only 
for the obvious replacement of the field operators by field functions. 

On the other hand at A 0, the discrete Lax operators get affected by reducing 
to the continuum limit 

L„(A, r],A) = I + AC{x, A, r/) + 0(A), (2.23) 
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while the i?-matrix remains the same for both discrete as well as the field models. 
For example, at this limit ( |2.31| ) turns into the well known Lax operator of the sine- 
Gordon field model ( 2.25| ), but with the same trigonometric i?- matrix ( |!2.18 ). 

At — > 0, when the deformation parameter q = e^"^ l+irj one gets an altogether 
different class of models. Since both R{^,hri) and ?7)-matrices depend on this 
parameter, both of them transform at this limit with the spectral parameter reducing 
to ^ = e^^^ ^ 1 + ii]X. The trigonometric i?-matrix ( p. 18 ) reduces clearly to the 



rational solution ( |2.19| ), while the L operator ( |2.21| ) reducing to 



L(A) = A + S 



3 

a=l 



a 



s'^) , (2.24) 

with an irrelevant scaling of the spectral parameter. S*" s° are simply the genera- 
tors of the undeformed s/(2). Expressing through spin or bosonic operators, ( |2.24| ) 
can generate Lax operators of the isotropic XXX spin chain and the NLS model. A 
list of such models, associated with the same rational i?-matrix is given below. 



2.1.3 List of quantum integrable ultralocal models 

We present for ready references a list of important ultralocal models with defining 
Hamiltonian or the evolution equation, along with their representative Lax operators. 
The 2x2 discrete Lax operators Ln{X) and the continuum ones C{X), given here are 
expressed through matrices 

2^ ^1^00^' 2^ ' yo 1 J ' \0 ) V^O/ 

and = — e"^. The matrix elements of the discrete Lax operators are realized in 
Pauli matrices ct°, bosonic operators g„, ipn with commutation relations [qn,Pm] = 
bPn, V'm] = x'^nm or in g-osciUators of different types. The corresponding 
continuum Lax operators C involve bosonic fields u{x,t),ilj{x,t) with [u{x),p{y)] = 
h6{x — y) or [ip{x),ip'^ (y)] = h6{x — y). 

These ultralocal models are quantum integrable and associated with either 
trigonometric or the rational i?(A)-matrix and ,in principle, exactly solvable by the 
quantum inverse scattering method (QISM) 0, based on the R and L matrix for- 
malism and the Bethe ansatz method described in sect. 7. The underlying algebraic 
structures, generation of the Lax operators and the construction of Hamiltonians for 
these models will be discussed in the following sections. However it should be men- 
tioned, that for some models the explicit quantum Hamiltonian is difficult to find 
and therefore for them we provide only the corresponding classical analogs. 
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I. Models with trigonometric Rtrig{X)-mSLtri'x 

Listed below are some field as well as lattice models, all associated with the 
trigonometric i?-matrix ( p^.lSI ). We have used the notations p = Ut, ^ = e~*'^ and A 
is the lattice spacing, 
i) Field models: 



1. Sine- Gordon model (SG) |2 



2 

m 



u{x, t)tt -u{x, t).j;^ = — sin riu{x, t) , 

V 

j^sg = ip + m sin(A — r)u)a~^ + m sin(A + r)u)a^ (2.25) 
2. Liouville model (LM) [|l|] 

1 I V ^e^^ 

u{x,t)u - u{x,tU = -e2^"(^'*), A™ = W -pi- ^2.26) 



4. Derivative NLS model (DNLS) H 

iiplx, t)t - ip{x, t)xx + Ai{ilj\x.t)ilj{x, t))ip{x, t) = 0, 

Cdnis = % + fc-^V e') + + ii,a- + - Ki^H) e') • (2.27) 

5. Massive Thirring model (bosonic) (MTM) ||18[ 

L = I dxijil^d, - m)^ - ^gj'^j^ = ^j'^^, [V^a(x).^i(y)] = fi5ab5(x - y) 
ii) Lattice Models: 

1. Anisotropic XX Z spin chain (related to 6-vertex model) [^] 

N 
n 

L(r-)(0 = z[sin(A + r/aV^) + sinr^(aV,- + a-a,+)]. (2.29) 

2. Wu-McCoy spin chain (corresponds to asymmetric 6\^(1) model) (WM) PO], pT| 

^ i _ _ 

H = J2 cos ^(t^nO^n+i + ^n^n+i) + COS r/aj^d^+i + - sin ^(a+(T„+i - a„ a+^J , 

n ^ 
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3. Lattice sine-Gordon model (LSG) |^ 

+ a^mA sm{\ — rjUnjcr'^ + rri A sm{X + riUn)cr~), (2.31) 

1 

where = 1 + |m^A^ cos2ri{un + \] 

4. Lattice Liouville model (LLM) ||82| 



( ) where /(„„) = [1 + A^e-(-.«>li (2.32) 



5. Lattice version of DNLS model (LDNLS) ( g-oscillator model) 

(2.33) 

where A^, Aj^ are g-oscillators expressed through standard bosonic operators as 

6. Ablowitz-Ladik model (AIM) [|8l 



n ln(l + /l) 

4"^"^^ (0 = fr' + bla+ + + ^ e^) , (2.34) 



where bm, b\^ are the g-bosonic operators with commutators \bmi b\^ = h{l — 



7. Relativistic quantum Toda chain (RTC) (corresponds also to discrete-time Toda 
chain (DTTC)) 



H = ^2 (cosh 2riPi + rj^ coshr7(pj + Pi+i)e 



niPi+l-Pi) + {<li-<li+l) 



L([*^)(0 = (^(^e^^P" - + r/e''"a+ - r]e^"P"-'i" a'^ . (2.35) 
8. Discrete-time quantum Toda chain like model (DTTCL) |4^ 

Lf'\\) = ^[A + iVi + w4ei + ^cosh(r/(iVi-A))e2 

+ (3if{Ni)pia+ + aiCOsr] qif{Ni)a-, (2.36) 
with = f{N) = 2^^^^ i[2N + uj- t], -[uj- th],) 
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II. Models with rational i^^at (A)-matrix 

Here we list ultralocal models associated with the rational i?-matrix (|2.19| ). 
i) Field models: 



1. Nonlinear Schrddinger equation (NLS) ET 



(2.37) 



ii) Lattice Models: 



1. Isotropic XXX spin chain ||2^ 



TV 

.3 3 N 



n 

)(A) = z(A/ + (aV^ + a+a- + a-a+)) (2.38) 
2. Lattice NLS model (LNLS) g 

2 

s = -. (2.39) 



A2^(2s - A^/'V)^ A-s + A^/^V /' 



3. Simple lattice NLS type model (SLNLS) 

H = ^^{(j){k + l)^{k-l)-{N{k)+N{k + l))(t){k + l)^{k) + {3KA'^)-^N{ky) 
A k 

LniO = ^(^{-'- + ^)e' + ^e'-zK^{n)a^ + ZKh{n)a-^ (2.40) 
with N{k) = 1 + KA'^(f){k)^lj{k). 

4. Tofia chain (TC) (nomelativistic, time-continuous) 

^ = E {Ip' + e^"-"^^^) > 4"nA) = ( ,^ 1 . (2.41) 



5. Quantum Toda chain like model (TCL) |47| , p3 | 

Li (A) = (A + -fPiQi + iUi) + p,pia+ + aiq^a' + ^ (2.42) 

7 
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with constant parameters 7, aj, uii 

6. Tamm-Dancoff q-bosonic model (TDqB) ||52| 



An integrable g-bosonic model with q = e*^-deformed rational i?-matrix is given 
by the Lax operator 



L^(A) = ((1 - iXA + KN)f{N) (e^ + e^)-iy/Kca+ + i^ha-) , (2.43) 
with f{N) = involving Tamm-Dancoff type g-bosonic operators [|143 



Historically the above models and their Lax operators were discovered in isolated 
and seemingly unrelated way. However it will be shown in sect. 5 using their un- 
derlying algebraic structure, how these models are interrelated and how they can 
be generated systematically as descendants of some more general integrable ancestor 
models (see fig. 5.1). 



2.2 Integrability of quantum nonultralocal models: braided 
Yang-Baxter equations 

As already mentioned, there exists another large and important class of models, called 
nonultralocal models. Such models might have inherent nonultralocality in the sense 
that the basic fields may satisfy current like algebras with derivatives of (5-functions 
appearing in their commutators, or they may be like fermionic or anyonic fields 
noncommuting at different lattice points. Alternatively, the models with standard 
canonical fields might also acquire nonultralocality due to gauge transformation from 
ultralocal models. Note that , the gauge transformations in the quantum case unlike 
classical models, might change the ultralocality property leading to a nonequivalent 
system. 

Since the ultralocality condition is a key requirement in deriving the representa- 
tive QYBE, the above scheme must be significantly modified to make it suitable for 
nonultralocal systems. We introduce here a scheme with extended QYBE and other 
relations for describing nonultralocal quantum systems, serving to some extend as 
universal equations for a wide variety of such integrable models. We provide a list 
of important nonultralocal models at the end of this section, deferring the details to 
sect. 8. 



2.2.1 Braided structures and extensions of QYBE 

For understanding the algebraic structures in nonultralocal systems we return to 
the simple matrix relations (^]7,2^) and see immediately, that both of them fail 
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for nonultralocal quantum models, since now the Lax operators do not commute at 
the same as well as at different points. Therefore the algebra with somewhat trivial 
multiplication property (|2.8D needs to be generalized as 



{A(^B){C ®D) =ijBc{A{C ®B)D) (2.44) 

where the noncommutativity of B2, Ci should be taken into account. However the 
underlying coproduct structure, responsible for the globalization of QYBE must be 
preserved. Such an idea was realized in the braided algebra [^, ^ ^ and was 



implemented in formulating the integrability theory of nonultralocal models 
The basic point is to complement the commutation rule for the Lax operators at the 
same site e.g. the local QYBE with the braiding properties describing commutators 
at different lattice sites. The local QYBE must also be generalized to incorporate the 
braiding relations, such that the transition from local to the global QYBE as in the 
ultralocal case becomes possible. In general, since the braiding at different sites may 
differ widely with arbitrarily varying ranges, the picture might be too complicated for 
explicit description. Therefore we limit here to two types of braiding only, showing 
however that this limited choice is enough to describe most of the physical models. 



Referring to |]T9|, |2T| for more details we present the braided QYBE (BQYBE) as (see 
fig. 2.11) 

Ri2{\ - /i)Z2i^Li,(A)Z2iL2,(/i) = Z{2'L2j{iJ,)Zi2Lij{X)Ri2{\ - yu). (2.45) 

where the nearest neighbor braiding is described by matrix Z , while the braid- 
ing between other neighbors, taken to be the same for all pairs, are given by a 
single matrix Z. Note that these braiding matrices in general can also have spec- 
tral parameter dependence as Z21 = Z2i{\,fi), Z12 = Zi2(/i, A) and similarly 
Z21 = Z2i{X,fi), Z12 = Zi2{fi, A). We however suppress this explicit dependence in 
all the relevant equations for brevity. To describe nonultralocal systems, the BQYBE 
must be complemented by the braiding relations 

L2,+i(/i)Z2-/Li,(A) = Z2-/Li,(A)Z2iL2,+i(/x)4-/ (2.46) 

at nearest neighbor sites and by 

L2MZ^,'L^j{X) = Z:^^L^^{X)Z2iL2MZ2^ (2.47) 

with k > j + 1 answering for all nonnearest neighbors (see fig. 2.12). Note that 
along with the conventional quantum i?i2-matrix as in ultralocal systems, additional 
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-2^12, ^12 matrices appear here, which may or may not depend on spectral parameters. 
When Z matrices are spectral independent, the above BQYBE and braiding equations 
can be obtained through Yang-Baxterization PT| from the more fundamental spectral 
parameterfree algebraic relations like 



-R^l 2 -^2 1^ % -^2 1 ^ fj 
-^2i+1^21 ^ 



(2.48) 



etc. Note that the relations like ( |2.48|) are also obtained from the BYBE ( |2.45|) etc. 
at the limit A — ±oo. 



The braiding relations ( 2.46 - 2.47 ) are the generalizations of the ultralocality con- 
dition and bear the true essence of the nonultralocal systems, where nonnearest and 
nearest neighbor braiding can be treated differently. In analogy with the derivation 
of YBE (|2.13| ) as the consistency condition of the QYBE, one can derive also the 
corresponding relations for nonultralocal models consistent with ( |2.45| ). Such YBE 
type equations involving R and Z matrices are 



i?12(A— /i) = Zi3 -Ri2(A— /i) , Z13 Z12 -R23 (/^^ ^) = -R23(Ai— i^) ^12 ^13, 

(2.49) 

with similar equations relating i?, and Z . Moreover the Z, Z and R matrices also 
satisfy individually the YBE ( p^.l3 ). Strictly speaking, the validity of YBE ( |2.13| ) for 
the Z-matrix is a sufficient but not necessary condition. One gets rather the YBE for 
matrix -Ri2(A) = Zi2-Ri2(A)Z2i^. This typical situation can be seen in the example of 
nonultralocal quantum mapping listed below and discussed in sect. 8. 

The i?-matrices associated with the nonultralocal models are the same trigono- 
metric ( p.l8| ) and rational ( |2.19| ) solutions (elliptic solutions are not considered here) 



found earlier for ultralocal models. 



The set of relations ( |2.45| - p747D represent the universal equations for the integrable 
nonultralocal models within a certain class of braiding. We shall show in sect. 8 that 
much in parallel to the established ultralocal case, particular choices for i?, L, Z, Z 
can derive systematically from these relations the basic equations for concrete phys- 
ical models. It is readily seen that the trivial choice Z = Z = I reduces the above 
set into the standard QYBE (|2.9| ) together with the ultralocality condition, while 
the nontrivial Z's would lead to different types of nonultralocal models. For exam- 
ple, particular homogeneous braiding with Z = Z ^ I correspond to SUSY and 
anyonic models, while the choice Z = I with Z ^ I describes nonultralocal models 
like WZWN, quantum mKdV etc. with 5' function appearing in their fundamental 
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commutation relations. The case like Z = I , Z ^ coincides formally with the 



reflection equation |63 



For establishing the role of the extended QYBE along with the braiding relations 
for describing the integrability of nonultralocal models, we have to show first that the 
important transition from local to the global QYBE is possible for the monodromy 
matrix 

T^-'l(A) = LakiX)La,k-liX) ■ --LajiX) (2.50) 

acting in the space 7i = \4(8>Vfc_i®. ■ .^Vj. Similar to our derivation for the ultralocal 
models we start from the local BQYBE ( |2.45| ) at point j + 1, multiply both sides by 



Z2iLij(\)Z2iL2j{iJ,) (see fig 2.13a, 2.14a) and using the same BQYBE at j in the 
rhs transfer -Ri2(A — fi) to the extreme right (fig 2.14a-b). 

We use next the braiding for nearest neighbors ( p.46| ) (fig. 2.13a-b) and its 1 ^ 2 



interchanged form (fig. 2.14 b-c) to drag Lj through Lj^i in both the sides. Some Z 
matrices gets canceled with Z~^ allowing us to arrange together the Lax operators 
acting in the same space 1 or 2 to give again the same form of BQYBE for T^"'"^-'(A) = 
Laj+i{X)Laj{X) (fig. 2.13c, 2.14d). We repeat further the same steps multiplying by 
Z^^Lij_i(A)Z2iL2j_i(/i) again from both sides and using the similar arguments. We 
have to remember only that together with (|2.46| ) for Lj and Lj^i we have to use also 



the braiding relation ( p.47|) for commuting Lj+i and Lj^i. Continuing for j — k + 1 
steps we get the global BQYBE is the same form as the local one: 

R^2{\ - fi)Z^M'''\X)Z2iTt'\fi) = Z^,^Tf^'\^i)Z,2Tf'^'\\)R,2{X - fi) (2.51) 

confirming the existence of the coproduct structure of the underlying braided algebra 
in nonultralocal systems. 

For showing the quantum integrability one has to find the commuting transfer 
matrices through trace factorization, which becomes a nontrivial task due to the 
appearance of Z matrices in the BQYBE (|2.51|) . We shall see in sect. 8 how to tackle 
this problem and modify the BQYBE for periodic models. 

Investigations of some nonultralocal systems from other angles can be found in 



[nn, m, 



2.2.2 List of quantum integrable nonultralocal models 

Nonultralocal models are mostly nonfundamental ones, defined in the Hilbert space. 
They may correspond to genuinely integrable systems with spectral parameter de- 
pendent i?(A)-matrix and the Lax operator or one may be interested only in their 
nonultralocal algebra ( p. 48 ) without invoking any spectral parameters. We collect 
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such quantum integrable models, which can be described fairly well by the present for- 
malism described by the relations ( p.45[ - p^47D and list them below mentioning briefly 



the choices for their braiding Z, Z apart from the i?-matrix and the Lax opera- 
tors. Details of these models together with the systematic derivation of their basic 
equations will be given in sect. 8. 

I. Systems with spectral parameterless i?-matrix 

Since the algebraic structures of the models are of concern here, we may take the 
spectral parameterfree equations like ( |2.48| ) involving i?+ matrix ( p^.lSD for describing 
the nonultralocal features. 



1. Current algebra in WZWN model |61 



The model involves the nonultralocal current algebra 

L2iy)} = L,ix) - UyMx - y) + ^C5'ix - y) (2.52) 

with Ci2 = 2Pi2 — 1, where L = |(Jo + Ji) with = d^gg~^, being the current. 
For describing the discretized algebra one needs the spectral parameter independent 
choice Ri2 = Z\2 = Rq2i with Z = 1. 

2. Coulomb gas picture of CFT (CG-CFT)|8| 
The associated nonultralocal Lax operator is 

C,g{x) = P{x) - S+ = J2E^, {P{x)®,P{y)} = 6\x-y)Y.H,®H, (2.53) 

a i 

Integrable description is made by the solution R = Rqi2 with the choice Z = 1 and 



3. Braided algebra (BA) ^ 

The braided algebras described in sect. 3 can be reproduced from the spectral- 
free BQYBE ( p.48[ ) as a particular case, by choosing homogeneous braiding Z12 = 



Z12 = R21 ■ As R12 one may take various solutions of ( p.l4| ) like (p.l5|) , nonstandard 
(|6.26|) , Jordanian solution etc. to reproduce results of 



II. Models with rational i?rat(A)-matrix 

Nonultralocal quantum models associated with the spectral-dependent rational 
i?(A)-matrix ( |2.19|) are listed below. 



4. Nonahelian Toda chain (NATO) |^ 
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The Lax operator of the model is given by 
^r^(A) = f ^ - - I , A, = g,g,\ B, = g,^,g,^ (2.54) 

along with the rational -R(A) = ihP vr — A and the braiding Z = 1 and 

Zi2 = 1 + ih{e22 ® 612) ® TT 



5. Nonultralocal quantum mapping p7|| 

The system is associated with A dependent Lax operator L^^^^ = V2„V2„-i, with 

K = A„(A) I 1 + E ^^MEiA . (2.55) 

with A„ = XnEjsf^i + J2i^=i^ = A2„+2 = A , involving nonultralocal algebra 

[Vi,j{n),Vk,l{m)] = h{6„„m+lSk,j+lSi,NSl,l - Sm,n+lSi,l+lSk,NSjA) (2.56) 



For integrable description one needs the rational i?(A)-matrix with Z = 1 and the 

aN- 



spectral parameter dependent choice Zi2{u2) = I + ^ J2a ^ ^Na e, 



III. Models with trigonometric i?trig(A)-matrix 

Nonultralocal integrable quantum models associated with the spectral dependent 
trigonometric i?(A)-matrix ( |2.18| ) are listed below. 

6. Quantum mKdV model (QMKDV) [§§ 
The nonultralocal Lax operator 

CmkdvixX) = '-{-v-{x)a^ + C(r^), [v{x),v{y)] = 6'{x - y) (2.57) 
describes this integrable model with classical equation 

±vt + Vxxx - Qv^Vx = 0. (2.58) 

The quantum model corresponds to the trigonometric i?-matrix and Z = 1 , 

Zj\2 — Zj2\ — q ^ , 

IV. Models with homogeneous braiding 

Nonultralocal models with homogeneous braiding have no distinction between 
braiding at different sites. 
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7. Integrable model on moduli space (IMMS) ||126 | 

This integrable model, motivated by the Chern-Simons theory, may be introduced 
by the Lax operator 

Ln{X) = + XI, (2.59) 

where matrices are related to generating quantum algebra. This abstract 
nonultralocal model may be described by the trigonometric Rtrig{X) matrix along 
with the braiding Z12 = Z21 = Ri2q 

8. Supersymmetric models (SUSY) 

SUSY models can be defined from the nonultralocal point by choosing Z = Z = 
HVapQaPi where rjap = Caa ® 6/3^ and g = (— l)"'^ with supersymmetric grading a. 
i?-matrix is either rational or the nonstandard trigonometric type ( |6.26| ). 

9. Integrable Anyonic model (JAM) 

Generalizing the SUSY model one may choose Z = Z = J^Va/sgais, with anyonic 
phase 6 included in g^is = e^^""^ . Integrable anyons correspond to the rational R 
matrix , while their g-deformations are associated with the trigonometric case. 

V. Systems with specific braiding 

Here we list integrable nonultralocal models models, which correspond to braiding 
a bit different from Z, Z used in the present formulation. 

10. Kundu-Eckhaus equation (KEE) ||87|| 
Classically integrable equation 

#t + + vii'^^)^ + e\^^f^ + 2ie{^^),^ = 0, (2.60) 

corresponding to the Lax operator 

i{\-ei,H) \ ^2.61) 

r]2ijj — z(A — 6il)'^il)) j 

is nonultralocal as quantum model due to noncanonical field ip and can be de- 
scribed by the rational -R(A) matrix along with braiding expressed through Z = 
diag{e^^, 1, 1, e*^) . The success and difficulties with the model are discussed in sect. 
8. 

11. Reflection equation (RE) |6^ 
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The reflection equation may also be viewed as nonultralocal model given through 
a different type of braiding. However formally it fits with the present formulation for 
nonperiodic models with braiding Z = 1, Z = R{Xi + A2). 

We shall elaborate on all the above models in sect. 8 showing their systematic 
classification in fig. 8.1. 



fig.2.1 



Figure 2.1 Ultralocality condition on the Lax operators showing trivial comuu- 
tation of operators defined at points i ^ j ■ 



fig.2.2 



Figure 2.2 Local quantum Yang-Baxter equation ( ^.dj) . R, L matrices act at the 
intersecting lines with 1,2 representing auxiliary and 3 quantum spaces. 



fig.2.3 



Figure 2.3 Multiplication rule ( 2.1Q ) induced by the ultralocality condition of 
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the Lax operators. 



fig.2.4 



Figure 2.4 Flipping the R-matrix using QYBE (^.(\ ) at point i 



fig.2.5 



Figure 2.5 QYBE for the nonlocal objects Tt*^^*l(A) defined at two points i + 



fig.2.6 



Figure 2.6 a) Global QYBE for the monodromy matrix T(A). b) T(A) is 

a global object defined through the Lax operators Lj at points N to 1. 



fig.2.7 



Figure 2.7 Trivial commutation of transfer matrices r(A) = trT{X) at different 
A. This feature of integrability ensures involutiveness of conserved quantities. 



fig.2.8 



Figure 2.8 Yang-Baxter equation l \2.1^ ) showing equivalence of two scattering 
processes to reach from the arrangement of particles (1,2,3) to (3,2,1). 



fig.2.9 



Figure 2.9 Vertex model with periodic boundary condition. Due to repeated 
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summation indices ( H-space indices), the transfer matrix r depends only on the 
quantum indices . 



fig.2.10 



Figure 2.10 Boltzmann weights of the Q-vertex model with charge conservation 
and charge inversion symmetry, constituting the elements of the R-matrix 



fig.2.11 



Figure 2.11 Local braided QYBE ( ^.4W - Intersecting (wavy) lines 1,2 with 
crossing from below denote Z12 or Z12, while from above denote their inverses. 



fig.2.12 



Figure 2.12 Braiding relation l \2.4'V showing nontrivial commutation between 
Lax operators L2fc(/i) and Lij{X) at different points k > j + 1. For k = j + 1, i.e. for 
\2.4(^ , in the Iks of the figure Z should be replaced by Z . 



fig.2.13 



Figure 2.13 Successive steps of the BQYBE for transition from its local to the 
two point form for nonlocal T'-'+^'-'l (The figure shows the Ihs of the equation). 
a.)Local BQYBE for Lj^i{X) multiplied by Lax operators at site j. 
b) Use of braiding relation ( ^.4(\ ) to drag Lj through Lj^i{ use fig. 2.12 in the boxed 
portion of a)-b)). 
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c)BQYBE forT^^+^^^ = L,j+i{\)Lj{\). 



fig.2.14 



Figure 2.14 Successive steps of the BQYBE for transition from its local to the 
two point form for nonlocal T'-'"^^'-'! (The figure sliows the rhs of the equation), 
a) Local BQYBE at j + 1 multiplied by objects at j. b) Flipping of R-matrix using 
BQYBE at j (use fig. 2.11 in the dashed boxes in a)-b)). 

c) Use of braiding (^-jOj with 1^2 (apply fig 2.12 in line boxes in b)-c). 

d) Transition to BQYBE for T^^+^^\( Use of Z12Z12 = I unknots the lines simplifying 
the picture of b)) . 



3 Algebraic aspects 

In this more or less self-contained section we present some useful and elegant algebraic 
structures hidden behind the integrable theory and intimately connected with it. 
The application of these formal objects for constructing basic relations in quantum 
integrable systems, will be considered in subsequent sections. Though the subject 
matter is rather abstract, we try to maintain our presentation at simple physical 
level. 



3.1 Sklyanin algebra 



We focus on the local QYBE (p|) R{X, fi) Li{\) L2{n) = L2{fi) Li{\) i?(A,/i), 
the key relation in quantum integrable ultralocal systems and consider first the gen- 



eral i?-matrix solution of p.l3| ) given through the elliptic functions: 

R{\-fi) = 
UJ0 + UJ3 = a(A) 
UJ1+UJ2 = c(A) 



IJ,)a' (g) a\ a" = I 
sn(A + a), uo — UJ3 = b(A) = snA, 
sna, ui — UJ2 = d(A) = k sna snA sn(A + a) (3.62) 



The corresponding solution for the L(A)-operator may be chosen in analogy with 
(I^D as (12:211) : L(A) = Sf^o ^i(A)c^' ® S\ The algebraic properties of S\ 
which are yet unspecified operators, are to be determined from the QYBE. Inserting 
the above R{X — /i)-matrix and the L(A)-operator in QYBE (|2.9| ), we notice that 
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except operators 5* all other elements entering in this equation are known objects. 
Comparing therefore the terms with independent functions of A, /i, we arrive at a set 
of algebraic relations for S^, free from spectral parameters: 

[S', = tJfs-ylS'', S^]+, [^^ S^] = i[S', S^]+ (3.63) 

involving commutator [ , ] and anticommutator [ , ]+ along with the Casimir 
operators 

C, = J^UiS'Y, C, = T.l^MS'^f. (3.64) 

The triplet of indices (a,/5, 7) represents any cyclic permutation of (1,2,3) and 
the structure constants Jas = -i — I may be parametrized by J^n = -^—^ . and 
obey the constraint J12 + J23 + J31 + J12J23J31 = 0. Note that in contrast to a Lie 



algebra, ( |3.63D is a quadratic algebra, which depends on the parameters a, k. This 



quadratic algebra, discovered in |31] is known as the Sklyanin algebra 



3.2 Quantum algebra and its origin in integrable theory 

We will be interested mostly in the trigonometric limit: A; — of the Sklyanin algebra 
(p.63| , |3.6^) , since at this limit it acquires a clear Hopf algebra property and gets linked 



with the quantum algebra. 

Since at — > 0, snx — >■ sinx, cnx cosx and dnx — > 1 , the elements of the 
i?-matrix (|3.62| ) turn into trigonometric functions as 

a(A) = sin(A + a), b{\) = sin A, c(A) = sin a, d{X) = 0. 

The corresponding trigonometric limit of the Sklyanin algebra may be obtained from 
the limiting values of the parameters giving the structure constants 

J12 = 0, J23 = -J31 = tan — . (3.65) 



This clearly reduces ( p.63| ) to 



= Ttan2|[S^5^] + , =45°5l (3.66) 



The Casimir operators ( p.64| ) reduce to 



^ = ^l,{Sr, C,-C2-^C2 = cos' I {Sy + sin^ | {S')'. (3.67) 
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This trigonometric Sklyanin algebra (TSA) is also quadratic and depends on a single 
parameter a. 

Instead of obtaining the TSA as a limit of ( |3.6ij| ), we may get it directly by 
using the trigonometric solution of the i?-matrix and constructing in analogy the 
corresponding L operator. To complete this chain of relationships, we should connect 
the TSA with the well known quantum algebra (QA) and that would link the origin 
of the QA to the Yang-Baxter equations and thus to the quantum integrable theory. 
Let us notice first that the Casimir operator given in the form of an ellipse suggests 
a simple reahzation (putting C2 = \l) 



go 



cos[as ) 
2cosf ' 



sm[as j 
2sinf 



(3.68) 



through a single operator s^. This realization together with the same S^, reduces 
the quadratic algebra ( p.66| ) to 



[s^S^]=±S^, [S+,S-] = [2s%. 
The remaining Casimir operator in (|3.64]) reduces to 

mj + i], = [s\[s'-i],+s+s^. 

where the q — bracket is defined as 

— q~^ sin (ax) 



x\ 



q-q 



-1 



sm a 



This in fact is the well known quantum algebra Uq{sl{2)) 
our aim of finding the origin of QA in integrable systems. 
At the limit g — >• 1 (or equivalently at a — > ) , when 



(3.69) 



(3.70) 



30[], and that fulfills 



a 



2s-^ + —{2s-'y + ■ 



the deformed algebra ( |3.69| ), clearly goes to the familiar su{2) spin algebra 



±s~- 



[s 



+ S-] 



2s' 



(3.71) 



The related Casimir operator reduces to j{j + 1) = s^{s^ — 1) + s~^s~ . We should 
remember that the quantum algebra, sometimes called as the quantum group, is in 
fact neither quantum nor a group. Nevertheless, it has got some analogy with the 
quantization procedure (see for example [137] as well as with group properties. 
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Quantum or quantized algebras may be introduced in different ways. However the 
way it has been obtained here hnking it with the integrable systems, gives perhaps 
more physical insight of this abstract mathematical structure. The g-deformation 
may be viewed as g = e*^^ by relating the deforming parameter rj to the coupling 
constant of the physical models and h to the Plank constant like quantum character- 
istic. We will see that the QA is not given by just any deformation, but is constrained 
by important algebraic properties like the Hopf algebra and that makes such defor- 
mations a unique one in the simplest case. Note that even in our familiar example 
(p.71|) , the sum of two commuting spins yields again a spin solution, i.e. 

As^ = <^ I + I <^ s^, As^ = (g) / + J (g) (3.72) 



also satisfies (|3.71|) . In mathematical language this is the definition of the coproduct 



A, which takes an algebra A to its tensor product and represents the main feature of 
a Hopf algebra. The familiar spin algebra therefore is a Hopf algebra, though only a 
cocommuting Hopf algebra (characterized by the explicit symmetry of the first and 
the second spaces in relation ( p.72| )). However, the superposition rule ( |3.72| ), as can 
be easily checked by direct insertion, is no longer true for the quantum algebra (|3.69|) . 
Therefore the idea is to generalize this rule as 

AS^ = ^q'' + q-'' ^ S^, As^ = s^^I + I^s^ (3.73) 

to make it compatible with (|3.69| ). This generalization indeed is quite nontrivial. 



since it defines now a nonco commuting Hopf algebra (notice the asymmetry in the 
two spaces in the coproduct relation ( ^.73 )) and makes the quantum algebra Uq{sl{2)) 



unique. 

A simple but fascinating effect of this nontriviality may perhaps be given by 
the following example 0. In anisotropic spin chain H^xz = J '^nii^n^n+i + 
^n^n+i) + "^osh a;cr,^cr^^_]^], the rotational nonsymmetry is immediate. However, 
if we spoil the symmetry still further by adding another term to it as H = 
Hxxz + sinha(crj^ — c^), curiously the rotational symmetry is restored again. The 
fact is that, the Hamiltonian H is related to the Casimir operator of the correspond- 
ing quantum algebra and consequently, it is symmetric under Uq{SU{2)) and hence 
naturally under SU{2). This type of Hamiltonian has been found useful in solving 
open chain problems p9|. 



3.2.1 Quantum algebra as Hopf algebra 

In a more formal approach quantum algebra Uq{g) can be defined as a deformed (or 
quantum) universal enveloping algebra of some Lie algebra g. Example is Uq{sl{2)) 
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as given above. With the following properties it turns also into a noncommuting and 
nonco commuting Hopf algebra A. 

i. Coproduct A : A ^ A® A, example (|3.73|) . 

ii. . Multiplication m : A ® A ^ A, example ( p.Sl) . 

iii. Counit e : ^ ^ K, example e(S'^) = 0; e(g^^^) = 1. 

V. Antipode (~ inverse) s : A ^ A, example s{S^) = —q^^S^, s(g^'^^) = g'^'^^. 



3.2.2 Quantum algebra through noncommuting geometry 



This is a group like approach starting from the noncommutative coordinate space 
having the relation xy = q~^yx, which can also be expressed as R~^{x^x) = (x^x) in 

X 



the matrix form. Here x 



y 



and = PR~^ is the braid group representation, 



while is the upper-triangular solution ( 2.15 ). If we demand now the invariance 
of the above defining relation for x under the transformation x' = Tx (recall that 
the vector transformation x' = gx, g G (rotation group) in usual commuting space 
preserves the length of vector x), we obtain the relation on T as 



R+ T1T2 = T2T1 R-^ 



(3.74) 



which is a group like relation of quantum algebra A* in the matrix form (recall the 
property g^g = 1 for a rotation group). The relation (|3.74|) written in the elemen- 

a b 



twise form through noncommuting matrix elements of T = with quantum 

\ c d j 

determinant detqT = ad — qbc usually taken as the definition of the quantum group: 

ah = qba , ac = qca , be = cb , bd = q^^bd , 

dc = q~^cd , ad — da = (q — q^^)bc. (3.75) 

Let us stress again that, though this is not a group, the property that the coproduct 
AT = T ■ T', for [T, T'] = satisfies again the relation (|3.74|) bears similarity with 
the standard group multiplication rule. 

In recent years success in establishing conventional exponential relation between 



the quantum algebra and group has been achieved [B8 



3.2.3 Faddeev-Reshetikhin-Takhtajan (FRT) formulation 

The FRT formulation |^ , closely related to the theory of quantum integrable sys- 
tems, bridges between above two approaches by noting that the algebras A and 
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A* are dual to each other. Coproduct A in A* induces multiplication m in A: 
(/i/2,a) = hkia) = (/i(g)/2)(A(a)) , where hj2 e A, a e A*. Setting e M„(^) 
the duality condition links between them as ^L^^^T^ = R^. The main relations 
of the FRT algebra may be given by 

lS±) = L?) lS±) R+ , R^ Lf^ L?) = L^p if^ R^, (3.76) 

and an equation similar to the first of ( p.76|) with R^ replaced by R~ . The R^ 
matrices are same as R^ in ( [^.15D , L"^ are upper/lower-triangular matrices, the op- 
erator elements of which are expressed through the generators (/, s^, S^) of quantum 
algebra (|3.69|) as 



Thus (|3.76|) with ( p.77| , p.l5|) reproduces the quantum algebra ( p.69|) , though given 



in a matrix form similar to the defining relation for the quantum group ( p.74|) . 

Another important feature of the FRT algebra is that, it derives elegantly the 
coproduct structures of the quantum algebra ( |3.73|) . Defining operators = / 
• • • ® X" ■ ■ ■ / we may evaluate 

, (+) , {+) _ / -(g - q-') (S' ® g<+i + g-< ® ^"^ J 

(3.78) 

Comparing with (|3.73|) it is immediate that 

4- 4:L. ^ ^ f - -.r'''^''-"-" ] - ^45.«. (3.79) 

Similarly we get LiJ. i^Li+i = ^-^L-ni+i and also at the same time 

Ll^l+i -^In] = ^-^ln!n,+i5 (3.80) 

where A = VAV is another possible permuted form of the coproduct with V{x ® 
y)V = y ® X. 

Note that the same i?+-matrix ( |2.15| ) appears in some or other form in the variety 
of relations concerning quantum algebra described above. We will see in sect. 4 that 
the i?^ matrices representing structure constants of the quantum algebra enters also 
as the building blocks for the i?(A)-matrix constructions in integrable theory. The 
associativity of the FRT algebra gives the spectral parameter less YBE ( p. 14 ). 
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3.2.4 Multiparameter deformation of general reductive Lie algebras 

Apart from the simplest one parameter deformation of s/(2) discussed above, mul- 



tiparameter deformations of general Lie algebras can be formulated [^]. We give 
below a brief account of it. 

Let g he a reductive Lie algebra of rank A^. Namely g is the direct sum of, say, 
M simple Lie algebras plus an abelian center. Hi, H2, . . . , Hj^j is the basis of the 
Cartan algebra of which Hi{l < i < Ni) span the semisimple part and the remaining 
N — Ni number of Ha, {Ni < a < N) belong to the center of the algebra. Let a^ be 
vectors with components a^, such that a^k = for all k and (A^i < a < N) , while 
Qij = 2{ai ■ aj)/{ai ■ a^) with (1 < i,j < Ni) are the entries of the Cartan matrix 
related to the semisimple part. Let Xf be generators associated to the simple roots 
ttj. The universal deformation tlq{g) can be defined so that each simple component 
remains the same as that of the standard one parameter quantization with relations 

m 

[Xt, XT] = = 5,,m,^^ (3.81) 

22 

and f{ii)Xf = Xff(ii ± a^) plus the Serre relations 

0<k<n \ / an 



(3.82) 



with n = 1 — ttij and the notation qa = e^''(')^"''"'\ where (1 < p{i) < M) counts 
the number of simple components. Since p is constant on each simple component, we 
may use also the notation qp = e^p'-^'i . The remaining deformation parameters on the 
other hand can be relegated to the coalgebra structure defining the corresponding 
coproducts as 

A(Xf) = ®Af + (Af)-i®Xf 

A{Hi) = Hi®I + I®Hi, (3.83) 
where Af , containing the parameters Via, tit {tik = —tki), has the form 

The antipode {S) and co-unit (e) are given by 

S{Xt) = -qfMt SiH,) = -H, and e(X±) = e{H,) = 0. 
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This defines the deformation of general Lie algebra as a Hopf algebra and we will 
show below the existence of universal 7^-matrix in this case. 



Sometimes it is convenient to use q-deformed Cartan-Weyl basis for related 
constructions. Such generators satisfy the commutation relations 



75 ' 



(3.84) 

(Ip - %^ 

In the above expression, if the index 7 corresponds to the non-simple root = Yl,oii 
for certain simple roots a^, then = Y^hi and the relations between two basises are 

1/2 

(3.85) 



hi 



Hi, 



e±i 



ai ■ ai 



9? 



3.3 Yangian algebra 



The Yangian Y{g) |^ associated to a Lie algebra g is another deformation of the 
universal enveloping algebra of polynomial maps 5'[A]: C ^ g. The general represen- 
tation theory of the Yangian is discussed in |Q, though we will be interested only in 
its simplest realizations related to the rational i?-matrix solution of the YBE given 
by ( |2.19|) : -R(A) = A — rjP , due to its direct consequence in the integrable theory. 

If we consider the global QYBE ( p.ll| ) with the rational i?-matrix, the monodromy 
matrix T(A) satisfies the quadratic relation 

(A - fi)[T''\X),r\fi)] = 7]{T^{X)T'"^{fi) - T''\iJi)T''\\)) (3.86) 

and the quantum determinant defined as detgT{X) = T^^(A — ri)T^^{X) — T^^(A — 
ri)T^^{X) commutes with all matrix elements of T(A) . Normalizing detqT{X) = 1 
and expanding 



(3.87) 



fc=0 



we get a realization of Yangian Y{sl{2)) generated by the elements t"^-,. The algebraic 
relations given by (|3.86| ) are equivalent to 



[tf 



"(0)' ^(fc)J — " " ''(fc); 

\j.ab +cd] r ±ab A.cd ] ±cb ±ad ±cb ±ad (n qq\ 

Note that these are like recurrence relations and due to the constant quantum deter- 
minant all generators of the algebra can be constructed starting only from its first 
two components t^Q^ and ff^y Clearly the Yangian algebra (|3.88|) is a deformation of 
the loop algebra. We shall see in sect. 5 some useful finite dimensional representation 
of this algebra realized in integrable models. 



cd ] 
{k)\ 



^cb ^ad 



^ad ^cb 
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3.4 Universal 7?.- matrix 

Notice that for deriving the algebraic relations like L^^. L')^_ = L')^_ L^^. Ri2~^ 
in the FRT formulation (|3.76|) we use the i?i2-matrix acting only in the auxiliary 



spaces, while the operators Lam defined in both the auxiliary space ( a = 1 
or 2) and the quantum space Ui. For obtaining the coproducts we need to multiply 
L operators in the same auxiliary space, but in different quantum spaces ?T,j,nj+i. 
However if we allow certain abstraction and suppose the existence of a T^niWi+i-niatrix, 
such that 'R-n,n,+^ ^Sn? ^i^+i = ^i2+i ^n.n.+i, then we can have, as seen 



from dpg), the relation 7^„^„,+l A (lJ^^^^^^)) = A 7^„^n,+l• In 



operator form such 7?.-matrices satisfy therefore the intertwining relation 

7^ A(a) = A(a)7^. (3.89) 

between coproducts A and A. This abstraction of the quantum i?-matrix satisfying 
(p.89|) is called the universal 7?.-matrix. A Hopf algebra endowed with a universal TZ- 



matrix is called a quasitriangular Hopf algebra. As a consequence of the coproduct 
property we then have 

(A ® J)7^ = 7^l37^23, (/® A)7^ = 7^137^12, (3.90) 

which would yield 

7^(A ® /)7^ = 7^l27^l37^23, , (3.9i) 

where we have denoted the quantum space indices by bold face numbers as in sect. 
2, for convenience. Notice the ordering of indices in the two relations in ( p.90|) , using 
which together with the intertwining property (p.89|) we get 



7^( A ® /) (7^) = ( A ® /) (7^)7^l2 = 7^237^l37^l2 (3.92) 

Equating ( p.91| , p.92D we get finally 

7^l27^l37^23 = 7^237^l37^l2 (3.93) 

as abstract generalization of all YBE's. Therefore this universal equation and the 
universal 7^-matrix ,as we will see below, may be taken as the basic objects for con- 
structing other important structures like FRT algebra, L^, matrices etc. They in 
turn would generate the key elements in integrable theory, like the spectral parameter 
dependent Lax operators L{X) and quantum i?(A)-matrices as well as the QYBE and 
YBE. 
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Though in general it is difficult to find the solutions for the universal 7^-matrix 



is the explicit form, following the prescription of Drinfeld p5| such an object may 
be obtained exploiting basically the defining relation ( |3.8y| ) and using ( |3.73| ). The 
simplest of such solutions for Uq{sl{2)) is given by |^ 

n = qhi-'^-') expg-2((g - q-')S-q''' ® q-""' S+) (3.94) 

with the notations expg(a;) = Y.n>o (fp where {n)g := ^fj"- 

3.4.1 Twisting transformation and general 7^- matrix solution 

For obtaining a general universal 7^-matrix solution corresponding to all deformed 
reductive Lie algebras discussed above, let us first look into the twisting transfor- 



mation introduced in Under such transformation one gets a twisted 7^-matrix 
starting from the original one: TZq — > 7^i2 = ^i2'^i2^v2 provided the twisting 
operator T satisfies the conditions 

7^12^ 13^ 23 = ^ 23^ 13^12) ^ 12^ 13^ 23 = ^ 23^ 13^ 12 (3.95) 

along with the antisymmetry property J^^h = ^ha- The corresponding Hopf algebra 
also gets twisted with a transformed coproduct A = J^AqJ^'^, such that the defining 
relations ( p.8£ ) and (|3.93|) are valid again for the twisted TZ matrix. 



We find first the TZq related to the semisimple Lie algebras. Implementing then 
suitable twisting Q we get the general solution for the reductive case with multipa- 
rameter deformations as 

'R-{qu,tij,Via) = G~^(tij,Via)'JZo{qii)Q~^(tij,Via). (3.96) 

We may construct the universal T^o-matrix for the semisimple Lie algebras in the 
form ||141| , pO| TZq = RK, where K is expressed only in terms of the Cartan 



generators as 

K = exp (^Y: Ki^^^^^^^d^' H,®H^ > (3.97) 

with d^^ = {d~^)ij, dij = (ctj ■ aj) representing symmetrized Cartan matrix. On the 
other hand R is expressed in the factorized form 

M 

R = UiR^'^)^ with = n (^i"^), ^i'^ = exp,- (S'(9P - ® e_,)) 

(3.98) 
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where A^J_ is the set of all positive roots belonging to the p-th simple component with 
the prescribed normal ordering |9^. For finding now the general solution (|3.96| ) we 
may choose the twisting operator as 

= exp{^( Eja v,^a^\Hkmc.-H^mk) + Y.^''''tkrar'H,mj)] (3.99) 

ijkr 

with a*-' = {a~^)ij. It can be checked easily that Q exhibits the required antisymmetry 
property Q21 = Q12 and satisfy the required conditions ( p.95| ). It is exciting to 
check, that the general solution IZ (|3.96|) with the insertion of (|3.97| , |3.98| , p.99| ) is 



the intertwining operator (|3.89|) for the deformed general reductive Lie algebra with 
coproducts ( p.83| ). 



Note that specializing to a particular Lie algebra, which would be reflected in the 
corresponding choice of root systems, Cartan matrix aij and the normal ordering in 
A+, one can derive from the general universal 7?.-matrix obtained here, the solutions 
for all the classical algebras of A,B,C,D etc. types together with their semisimple 
as well as reductive generalizations. For example in the case of quantum algebra 
Uq{gl{2)), by applying the twisting transformation JF = g*f(s^®^-^®«^) on (|3.94|) , one 
derives the universal TZ matrix 

^ ^ ^i2(s3«s3+e(s3«z-z«s3)) y g (q^'-^S+r ® iq-'-''^^^S-r, (3.100) 

where Z is central to the algebra. 

In the next section we take up important applications to integrable systems using 
different representations in the case Uq{gl{n + 1)). 



3.5 Hecke algebra, Birman-Wenzl-Murakami algebra 

Such algebraic structures inherent to matrices play important role in integrable 
systems. Hecke and Birman-Wenzl-Murakami algebras may be represented by the 
braid group along with certain extra conditions. The braid group in turn is given by 
the relations 

hi hi+i hi = hi+i hi hi+i, hi hj = hj hi, \ i- j \>2, (3.101) 

where hi, i = 1,2, ■ ■ ■ ,n — 1 are braiding operators deflning a group, the elements 
of which are n-braid. The braid group can be used for obtaining link polynomials, 
which have importance in the knot theory for classifying links and knots [Q. We are 
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interested however in the integrable theory and observe that, defining R = PRg 



the spectral-free equation ( |2.14| ) rewrites as 



Ri2 R23 R12 



R23 R12 -^23) 



[R12, R. 



34J 



0. 



(3.102) 



Representing Ru+i = h the equivalence of (|3.102|) and (|3.101|) becomes immediate. 

The Hecke algebra may be defined by the braid group relations ( |3.101D along with 
an extra Hecke condition 

h^i + {l-t)h, = t (3.103) 

with e = const. Notice that, if we demand the i?^-matrices to satisfy in addition to 
( p.l4| ) the condition 



R-^ - R- = {q- q-^)P, R' = P{R^)-^P, 



(3.104) 



then one finds easily that R — R ^ = {q — q or {qRY + {l — q^){qR)=q'^I 



which is exactly the Hecke condition ( p.l03| ). Though the condition ( |3.104D may 
appear to be rather restricted, it is interesting to check, that not only the matrices 
( p.l5| ), but also the whole class of i?^-matrices related to Uq{gl{n + 1)) satisfy this 
condition. In the next section we will be concerned basically with such i?^-matrices, 
since they through Yang-Baxterization allow an easy and systematic construction of 
Lax operators and quantum i?-matrices representing integrable systems. 



Birman-Wenzl-Murakami (BWM) algebra ||138|| is of more general structure and 
may be defined through the braid group relations ( p.l01| ) along with the condition 

i 



{Ri - ci){Ri - C2) ■ ■ • {Ri - Cn) = 0, 



(3.105) 



where Cr,r = 1,2, are constant eigenvalues. Note that in the restricted case 

of n = 2, condition ( p. 1051 ) reduces to ( |3.103| ), though otherwise it is more general 
and applicable to all other i?^-matrices not belonging to the Hecke class. Some 
prominent examples of such i?-matrices are those of the higher spin models as well 
as of the models related to algebras other than gl{N). For the Yang-Baxterization 
and other applications of the BWM algebra, the readers are referred to E^, H, |91||. 



3.6 Significance of Hopf algebra in integrable systems 

To clarify the significance of the Hopf algebra in integrable theory let us return to 
the coproduct and multiplication properties of the algebra and the transition from 
the local to the global QYBE, explained in sect. 2. 
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We will discuss first the ultralocal case. Let 5*" = / (g) / (g) / (g) ■ • ■ (g) S"" ®/ ■ • • (g) / 

"1 V ^ , 

n\ 

and similarly 5*^^ be two copies of generators of the quantum algebra, which trivially 
commute. Due to the Hopf algebra structure the coproduct AS'"^^^^) acting nontriv- 
ially on both the spaces ni and satisfies again the same algebra. Continuing such 
tensor products one obtains the coproduct acting globally on all the spaces (see fig 
3.1 a,b,c). Note now, that S^^^^) and Lgl^fn^). SJ^^)) satisfy the FRT re- 

lations L^^^ L^^. = L*^] Li^l for i = 1,2 due to the related S"^. satisfying 
the algebra ( p.69| ) or (|3.71|) . The ultralocality condition [Li^], = 0, holds 

due to the vahdity of [S^.S^J = . Finally, the coproduct ALj^^ma) = ^iS^iS 
defined on both the lattice sites (nin2) is again in upper (lower)-triangular form and 
due to the Hopf structure of the elements AS'"^^^^), becomes a solution of the same 
FRT algebra 

(A^Stn.)) (A^tn.)) = (A^tn.)) (A^Stn.)) R^- (3.106) 

Thus the coproduct A : A ^ A of the Hopf algebra induces nonlocal operators, 
continuing which as ^ ^ ^(g^(g- • -(g^, we get the global object AL^^ ^-^(AS'^" 
defined on the whole interval [1, A^] (see fig. 3.2 a-c). Such a globalization of the FRT 
algebra caused by the coproduct of the Hopf algebra (fig. 3.1 a-c), in fact represents 
in the Yang-Baxterized form the construction of the global QYBE from the local one, 
described in sect. 2 (fig. 2.2-2.6). Let us now focus on another important feature of 
a Hopf algebra, namely the multiplication rule m : A^ A ^ A and look into the 
ultralocality condition we have crucially used. Though its action is easy to define 
for the quantum algebra: m : (a (g 6) ^ a ■ b, this operation for the higher tensor 
products, is given in a bit nontrivial way involving the the permutation operator 
t: bxc^cxb m® = [m x m) : {id x t x id). Thus we actually arrive at our 
familiar relation (|2.8| ) as 

(a(g6) ■ (c(g(i) = (m X m)(a(gr(6(gc) (gd) = m(a(gc) x m{b®d) = ac®bd (3.107) 

where we have assumed the trivial permutation property r(6 (g c) = c (g 6, which is 
also true for the quantum algebra. Note that transition from local ( ^.76[ ) to global 
(p.l06|) and similarly for the QYBE's related to the ultralocal models mimic the 
multiplication property ( |3.107| ). 

It is remarkable that even for nonultralocal models, as we will see below, the Hopf 
algebra structure of the underlying algebra induces again the transition from the local 
to the global braided QYBE (fig. 2.13, 2.14 ) crucial for the quantum integrability. 
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3.7 Braided and quantized braided algebra 



In general the permutation property of r in ( p. 1071 ) may become nontrivial :t23 
and unlike quantum algebras can modify the multiplication rule as 

mBtg)B = {mx m) ■ {I ® 1^23 ® I) 



23 



(3.108) 



which is the main essence of the braided algebra |3^. Generalizing the notion of 
superalgebra, the braided algebra B{R) may be defined through the commutation 
relation 

R21T1R12T2 = T2R21T1R12 (3.109) 



for braided matrix T with constant solution R12 of YBE (p.l4|) . If T' is another 
independent braided matrix having the braided statistics 



(3.110) 



then A(T) = TT' is also a braided matrix satisfying the same algebra B{R). This 
shows that the braided algebra, though have complicated multiplication rule, has 
the coproduct structure, an important ingredient of the Hopf algebra. The matrix 
elements a, b, c, d of the 2 x 2-matrix T, T' generates different braided algebras for 
different constant i?-matrix solutions like standard, nonstandard, Jordanian forms 
etc. Il35|. 



Generalizing further to the quantized braided group QB{R,Z) |^ |3^, along 
with i?-matrix another Z-matrix were introduced with R, Z being the solutions of 
constant YBE's [RRR], [ZZZ], [RZZ] and [ZZR]. Here we have used the obvious 
notations for writing constant YBE's like ( ^^.141 ). The defining algebraic relations 
for QB{R, Z) may be obtained from equations like ( p.48|) as a particular case when 
Z = Z. This type of relations are responsible for the homogeneous braidings for all 
copies of the algebra. 

For incorporating different braidings for nearest and further neighbors, Z and Z 
matrices were included along with the R matrix in the generalization of the quan- 
tized braided algebra in Yang-Baxterizing the braided FRT relations given by 
(|2.48| ), one may arrive at the spectral parameter dependent BQYBE's ( |2.45| - |2747| ) etc. 
suitable for describing integrable nonultralocal models. Therefore, in one hand such 
equations serve for describing physical models including their conserved quantities 
and on the other hand at the spectralfree limit ( |2.48| ) reproduce the abstract math- 
ematical structures like quantized braided algebra fS^] ioi Z = Z and the braided 
algebra for Z12 = Z12 



R 



'21- 
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The Hopf algebra structure of the braided algebra plays the key role through 
its coproduct in constructing the BQYBE for the monodromy matrix leading to the 
quantum integrability of the nonultralocal models similar to the ultralocal case. 



fig.3.1 



Figure 3.1 Coproduct property of quantum algebra. 
a)Generators S^. act nontrivially of the Ui-th space. h)The coproducts A^^^^^^^-j act 
nontrivially on both rii and 712 spaces. c) The coproducts AS'|\...^-) act globally on 
all the spaces. 



fig.3.2 



Figure 3.2 a-c) FRT algebra and its tensor product induced by the coproduct 
structure. 

4 Construction of models from universal 7^-matrix 

In sect. 2 we have given a list of ultralocal models along with their Lax operator 
L(A) and the associated quantum it!(A)-matrix. Here we focus on their origin and 
construct them for a large class of models. This scheme 

i) reveals the algebraic aspects of the ultralocal quantum integrable theory and 

ii) using underlying algebraic structures constructs spectral parameter dependent 
representative Lax operators in a systematic way. 

It also highlights the following basic questions pertinent to the integrable systems 

iii) explaining why a wide range of integrable models like sine-Gordon , derivative 
nonlinear Schrodinger equation, massive Thirring model etc. share the same trigono- 
metric i?(A)-matrix, while another large class including NLS model, Toda chain , 
isotropic spin-| chain etc. are all associated with the same rational -R(A) matrix (see 
the list in sect. 2). 

iv) Thus it focuses on the commonness between such diverse models and tries to find 
the guiding principle for classifying different classes of integrable models. 

v) It identifies an algebraic criterion for defining the integrable nonlinearity, i.e. the 
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nonlinearity that makes a system integrable. 



4.1 Useful representations of universal 7?-- matrix 

The aim is to start from abstract algebraic structures and build systematically the 
objects of direct relevance in integrable ultralocal models. We look first into the 
universal 7?.-matrix solutions like ( |3.94| , p.96| ) of the universal Yang-Baxter type 



equation ( p. 931 ) and denote them by TZ~^ and the related equation by (+ + +), since 
their finite-dimensional representations as derived for example from (|3.94| ), are of 
upper-triangular form. It is easy to check that, if we define 71^2 = (^^)2i) it would 
correspond to the representation in lower-triangular form and be a solution of the 
universal ( ) equation. More generally one would get the universal equations 

^12 ^13 ^23 ~ ^23 ^13^12 ) (4-1) 

of six different forms with el, e2, e3 occuring in the sequence + + +, H , -|- + 

— , , h, — + + . The validity of these equations can be proved starting from 

the original one + + +, taking the inverse of TZ and renaming the quantum spaces 
1, 2, 3 in a suitable order. At the same time one concludes also about the nonvalidity 
of (± =F ±) type equations, which can hold only for triangular Hopf algebras with 

i?21 = Rl2- 

It is important now to show that, algebraic objects related to the quantum alge- 
bras can be generated from the above universal equations and the 7?.-matrices through 
direct reductions. For finite-dimensional fundamental representations vr acting in the 
space 1 and 2, while 3 kept as the quantum space, one gets 

(tt (g) 7^)7^l2 = Rt2, (tt ® / ® I)nf^ = Rf^ = L[f (4.2) 

which reduce (|4.1| ) clearly to the FRT algebra relations (|3.76|) . Consequently, taking 
the universal solution ( |3.94| ) and using the map vr one may get the i?^ and L^"^^ 
matrices in explicit form and hence derive the quantum algebra ( p.69|) from the 
quadratic relations among the elements of L^^)a3. At the same time why FRT relation 

like (H h) does not exist, is also answered promptly from the nonvalidity of such 

universal equation. 

Note that the complementary action of tt on 3 ^ 3 with 1,2 kept as quantum 
spaces, derives from the same universal relations ( [4.1|) the intertwining relation of the 
universal 7^-matrix ( p.89|) and from the solution ( p.94|) the coproduct structures of 



the quantized algebra (|3.79| ). The action of vr on all the spaces reduces (|4.1| ) to the 



spectral parameterfree YBE ( 2.14 ). 
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4.2 Yang— Baxterization and construction of i?-matrix and 
Lax operator 



It is crucial to observe that the FRT algebra relations (|3.76|) are quite similar to 



the QYBE ( p.9| ) and the algebra related objects and resemble closely to the 
Lax operator L{X) and the quantum i?(A, yu)-matrix associated with the integrable 
systems. However the main distinguishing feature of the former algebraic entries is 
their spectral independence. The spectral parameter on the other hand is an essential 
item in the integrable theory. Therefore, for application to integrable systems we need 
to introduce spectral parameters in the formal algebraic structures or in other words 
Yang-Baxterize the FRT algebra. 

To frame the scheme for generating integrable models we adopt sequentially the 
following steps. 

i) Start from the universal 7^-matrix and reduce it to i?^ and L^"^^ involved in the 

FRT algebra. 

ii) From the reduced universal YBE determine the underlying quantum algebra and 

the related coproducts. 

iii) Yang-Baxterize the FRT algebra to construct the Lax operator L{X) and the 

trigonometric Rtrig{X — yu)-matrix of a generalized and exact integrable model 
on a lattice, which may serve as an ancestor model. 

iv) Find different realizations of the related quantum algebra and reduce conse- 

quently the ancestor model to generate explicit Lax operators for variety of 
integrable lattice models. Since the Yang-Baxterized ancestor Lax operator is 
a solution of the QYBE, all the models generated from it must also be quantum 
integrable and share naturally the same quantum i?-matrix. 

v) Continuum models are obtained from their lattice versions at the continuum limit 

A ^0. 

vi) For generating the class of models associated with the rational Rrati^ — f^)- 
matrix, we may take q —>■ 1 limit in the above scheme. As a result the underlying 
quantum algebra becomes a realization of the Yangian algebra and one recovers 
consistently the relevant objects for the rational class of models. 

This gives a systematic scheme for constructing the Lax operators and the quan- 
tum i?-matrices of different families of exactly integrable quantum models of ultralo- 
cal class, both at the lattice and the continuum level. However, as we will see below. 



50 



one encounters certain practical shortcomings of this scheme and has to prescribe 
some working remedies. 

For demonstration we take the simplest f/g(s/(2)) case and as itemized above, 
look first for the universal 7^-matrix solution (|3.94| ). Noting that, since 7f{S^) = 
0"^, 7r(s^) = maps into Pauli matrices, we obtain using ( [4.2| ) explicit forms of L*^^) 
and reproducing the FRT algebra with ( p.l5|) and (|3.771 ), which yield naturally 



the well known quantum algebra (|3.69 ) as the underlying algebra. 



The Yang-Baxterization of the FRT algebra , i.e. introduction of spectral param- 
eters C, = e*'''*', C = e*''^, is in general a difficult task. However in the present case it 
takes the simple form [03 



R{X - /x) = - ^R- , (4.3) 

for the quantum /^-matrix and similarly ^ 

L(A) = + ^L(-) , (4.4) 

for constructing the Lax operator. For showing that R{X — fi) and ^^(A) thus obtained 
indeed represent integrable systems, one has to check them as solutions of the QYBE. 
Therefore inserting the above Yang-Baxterized forms in ( p.9| ) and comparing the 
coefficients of various powers of spectral parameters, we arrive at a set of seven 
algebraic relations involving i?^ and L^"^^ . Six of them are same as the FRT 
relations ( p.76| ) and hold automatically, while the remaining one is 

- = R~L^^^Li-^ - L^f^L^^^R' (4.5) 

As we see, it is a combination of (± =F ±) type relations and clearly does not belong 
to the FRT algebra as explained above. 

However the key observation is that, if we demand the i?^-matrix to satisfy the 
Hecke algebra (|3.104|) , we may replace in ([4.5|) R^ through each other, using the 
Hecke condition R^ = R~ + c P, which results 

i?-LS-)4+) - 4+)LS-)i?-(cP)LS-)4+) - L^^^l[-\cP) 

= i?+LS+)4-) - 4-)^;+)/?+ + (-cp)lS+)4-) - Li~^L^^\-cP). 

The above relations hold obviously due to the permutation property of P and the FRT 
relations ( |3.76D . We see thus that the FRT relations can be Yang-Baxterized as 



jj) to give the QYBE for the class of i?-matrices belonging to the Hecke algebra. 
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At the same time one generates also the spectral parameter dependent genuine -R(A)- 
matrices and the Lax operators of quantum integrable systems. As mentioned, all 
i?-matrices related to the general Uq{gl{n + 1)) exhibit the Hecke property. 

It is illustrative to check that, the i?(A)-matrix built through Yang-Baxterization 
( [4.3| ) from the basic solutions ( |2.15| ) has the form 

Rll = Rll = sm{\ + a), Rll = Rll = sm{X), Rll = smae'^, Rll = smae-'^ 

(4.6) 

where a = hrj and A is scaled by rj. Notice that this is the same trigonometric 
solution ( |2.17| , |2.1g| ) obtained as a direct solution of YBE ( ^.13| ). Similarly the spectral 
dependent Lax operator, Baxterized as (|4.4| ) from L^"^^ (|3.771 ): 



sin(A + ?7s^), sinrjS e*'^ 
sinr^S'^e"*'^, sin(A — ?7S'^ 



L(A) = I rrol^-a 3, 1 (4-7) 



coincides with the expression ( |2.21| ), when the Sklyanin algebra is realized through 
the generators of the quantum algebra as (|3.68 ). It should be remarked, that the 



above Yang-Baxterized R{X~ ^) and L{X) matrices are in fact slightly different from 
their earlier forms, while the exact coincidence is achieved by a simple and allowed 
similarity transformation 

RuiX) ^ Au{X)Ri2{X)Ai,\X), and L^^iX) ^ Ai2{X)L,2{X)A^^{X) 

with Ai2(A) = e'W^^-^^^'l 

Now we come to the points v) and vi), i.e. to the practical applications of the 
scheme for constructing explicit Lax operators of integrable lattice as well as field 
models. For this we find different realizations of the underlying quantum algebra 
Uq{sl{2)) ( p.69| ) generated by S'^,s^. Inserting them subsequently in the Baxterized 
L(A) ([4.7|) , which serves as an ancestor model, we generate easily the representative 
Lax operators of other models. Since the ancestor model is quantum integrable with 
trigonometric i?tr.ig(A)-matrix ([4.6|), all its different descendants obtained through 
particular reahzations would evidently be quantum integrable and related to the 
same Rtrig- 

4.3 Models with trigonometric i?trzg (A)-matrix 
4.3.1 XXZ spin-i chain 

Simplest realization of the generators 5*^,5^ through Pauli matrices gives 

0"^, = and since e^*'''^^ = cos?7 ± ia^ smt] one gets [(J^]q = cr^, which reduces 
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the quantum algebra (|3.69|) evidently to the spin algebra (|3.71| ). The Baxterized L{X) 
operator reduces correspondingly to the form 



sin(A + rja^), sin rja 
sin 770"^, sin(A — r^a" 

Note that this is the Lax operator of the quantum integrable anisotropic XX Z spin-| 
chain listed in (|2.29|) . The Hamiltonian of the model is derived in sect. 6 and the 



eigenvalue problem is solved in sect. 7. 

4.3.2 Sine-Gordon: lattice and field models 

An important realization of the Uq{sl{2)) quantum algebra in bosonic operators Un,Pn 
is 

e^"^' = V2^e''T',S- = g{u)e'^P,S+ = {S-)\ g{u) = [1 + 2s^ cos2r]{u + ^)]'^ , (4.9) 

where s = Using the commutation relation [un,Pm\ = -^Snm, it is easy to check 
that the mapping ( |4.9| ) satisfy the g-deformed algebra ( |3.69| ). Therefore we may 
insert ( |4.9|) in the general L{X) (|4.71 ) to obtain 



r(\^=l "^^sin(A + r/u„) ^(u„) e*^"^ 

with g{un) given as in ( [4.9|) . The lattice Lax operator generated by cr^L„(A) = L^'*^-' 
represents exactly integrable lattice version of the sine-Gordon model listed in (|2.31|) . 

At the continuum limit A — > 0, as explained in sect. 2, we obtain Lll^^\\) = 
I — iACgg. Using the limits up to the order A giving m„ u{x),Pn — * p{x), g{u„) 
1, e^*^"^ ^ I±iAp{x) one derives finally from ([4.1CI|) the well known Lax operator 
of the sine-Gordon field model as Csg{X) = p{x)(T^ + kQ smr]u{x)a'^ + ki cosr]u{x)a'^, 
with the momentum vector ko = m cos A, ki = im sin A satisfying the relativistic 
condition k^ — kf = . The energy spectrum and the physical significance of the 
SG field model are discussed in sect. 7. 

Note that both the lattice and the field SG models are associated with the same 
trigonometric R{X — /i) (|2.18|) obtained through Yang-Baxterization. It should how- 



ever be remembered that, while the lattice version ([4.10|) is an exact solution of the 
QYBE, the field Lax operator Csg satisfies QYBE only up to first order in A, which 
is enough for the field limit. Similarly, the quantum algebra realized exactly in the 
lattice SG trivializes also at the continuum limit, though with some mysterious affine 
quantum group symmetry appearing in the conserved quantities [^, which goes 
beyond the present scope. 
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4.3.3 g-oscillator model 

g-oscillators or g-bosons may be defined through the commutation relations 

[A, N] = A, [A\ N] = -A\ AA^ - q^A^A = q^^ (4.11) 



An interesting realization of the quantum algebra ( p.69|) in g-oscillators may be given 



through the g-Holstein-Primakov transformation ||146 



5+ = ([2s-n],)U, S- = A\[2s - n]g)^ , = s - n, A^A=[n]q (4.12) 

while another simpler one through 

= N+^, 5+ = nAq^, S- = Kq-^A\ k = {q - g"^)i (4.13) 

Using such realizations it is possible to construct from ([4.7|) quantum integrable 
lattice models involving g-oscillators. 

4.4 Multicomponent models 

For constructing N (g) N matrix Lax operators and ® /2(A)-matrices associated 
with multicomponent integrable models, we may start with the universal T^o-^iatrix 
G Ug{sl{n + 1)) (g) Ug{sl{n + 1)) given by (|3.97] , ^l9^ ), where N = n + 1. For finding 



finite-dimensional fundamental representations ([4.2| ) of TZq, we may use the mapping 

niH,) = Eu-E,+u+i, 7r(g±^0 = l + ^n(#-l) + ^mmkf-l) 
7r(X+) = Eu+i, 7r{Xr)=E,+u, 7r{X±) = Ei,+,, 7r(X,^) = (4.14) 



where {Eij)ki = SikSji, along with the relations (p.85|) . Since the i?^-matrices obtained 



are again of the Hecke type, we may follow the steps (5]^,0) to construct the spectral 
parameter dependent operators 

n+l n+1 

R{X—fi) = sin(A— /i+a) ^ Ei,k^Ekk + ^{sin{X—fi)Ejj^Eii+sinaEij^Eji) (4.15) 

k=l iy^j 

and 

n+l 

LiN.N)^^^ = Y.i^uJk + ^'uJ,')Ekk + smaY.^E.,,_,,uJ,e,,-r'EJ+l^e^,uJi' (4.16) 

k=l i<j 

where 

ujk = q^^^=^ ~'^^='^~>-, co-n+i = g~ ^'=1 ~ , e±i = q~'^Xfq'^~ (4.17) 
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for /c = 1, ■ ■ ■ , n and the operators Cij correspond to the nonsimple roots I]s=j (^s- For 
generating muhicomponent integrable lattice models from ([4.16| ) associated with the 
quantum i?-matrix ( 4.151 ), we have to consider different realizations of the underlying 
quantum algebra in more physical variables like bosonic, g-bosonic or spin variables. 



4.4.1 Generalized spin model 

For fundamental models the Lax operator ( |4.16| ) coincides with ( ^.15| ) due to the 
mapping (^4.14 ) and represents quantum integrable model with interacting generalized 
spin operator Eij ||9^. Due to the regularity condition R{\ = 0) = P the Hamiltonian 
of the model can be constructed following sect. 6 and the eigenvalue problem can be 
solved using nested Bethe ansatz described in sect. 7. 



4.4.2 Multicomponent g-bosonic model 

Consider g-bosonic realization of the higher rank quantum algebra Uq{sl{n + 1)): 



n+l 



Hi = iVi-iVi+i, N^ = {n + l)s-J2Nkior 1 = 1,2. ..,n 

k=2 

Xt = 4^fc+i, X+ = ([iVi],)U2, X,- = (X+)t (4.18) 



with k = 2,3 ... ,n. Here s is an arbitrary parameter and Al, are n number of 
g-bosonic operators with deformed commutation relations 

COS 2 

With such realization and relation (|4.17| ), ( [4.1 6| ) represents the Lax operator 
-Z^n(A)^'^"'**~'''' of a quantum integrable multimode g-bosonic model EH. 



4.4.3 Toda field model 

Let us consider a realization of the above g-bosons in canonical variables [u'^jP^] = 

A, = q^P^ [iV4, 4 = [N,]g q-^P^ Nk = akU + s (4.20) 

with i = 2, 3, . . . , n+l and k = 1,2, . . . , n+l, where a, and pi are the simple roots and 
fundamental weights of s/(n + 1), respectively with the relation (oj ■ pj) = 6ij, i,j = 
2, 3, . . . , n + 1. Here we have assumed di = — Y^^^2 Pi — 0- curious 

to note that this realization of g-boson reduces the Lax operator L„(A)''''""'**~''^ or 
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equivalently, the original L-operator ( [4.16| ), with the use of ( [1.18| ) to the well known 
lattice regularized version of the Toda field model L„(A)(*°'^''^*"'"') g. Obviously for 
n = 1 we recover the lattice SG model. 

By rescaling p Ap and adjusting the parameters one can introduce the lattice 
constant A and go to the continuum limit at A 0, to generate the Lax operator 
of the Toda field models. 



4.5 Models with twisted i?^^j^(A)-matrix 

For constructing more general integrable models we may consider TZ G Uq{gl{2) 
Uq{gl{2)) constructed in ( |3.10CI| ). However instead of repeating the Yang- 
Baxterization [0, we can obtain the spectral parameter dependent i?-matrix and 
the Lax operator by twisting the previous Yang-Baxterized forms (|4.6| , |4.7|) with 
suitable representations of the operator JF = g«^'(s^®^-^®s^)_ This gives through 
Rtigi^ - = F{e)Rtrig{X - /i)F(^) with Fuie) = e^^(-'®^-^®-') the twisted quan- 
tum i?f^jg(A)-matrix 

= Rll = sin(A + a), Rll = sin(A)e~^^ RH = sin(A)e*^ R^ = R^ = sin a 

(4.21) 

and similarly the twisted Lax operator as 

L,„(A, e, Z„) = F,^{e, Z^)L,^{X)Fan{e, Z„), r^nie, Z^) = diag (e^^(^"-^^), e~^^(^"+^ 

(4.22) 

where Lan{X) is given by ([4.7|) . We have chosen here n{Za) = I, a = l,2to exclude 
the color parameters, which will be considered below. 



4.5.1 Generalized Wu-McCoy model 

As examples of concrete models we may consider simply the fundamental represen- 
tation of the Lax operator ([4.22|) coinciding with ([4.21|) , which would correspond to 



the spin Hamiltonian 

H = Y.l cosaalai^, + e^^^^+^^+^^V+a^+i + e-^(^-+^-+^)V^a++i. (4.23) 

m 

Note that for Z^ = 1 it reduces to the Wu-McCoy model ( |2.30| ) of ferroelectric in 
a constant external electric field 6, the Hamiltonian construction of which is dealt 
in ( |6.18 ). Therefore, (|4.23| ) may be considered as the Wu-McCoy or the asymmetric 



6- vertex model ||144|| placed in an inhomogeneous electric field Z^"^^9 and solvable 
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through the Bethe ansatz. Similarly, one may choose more general twisting with color 
representations 7i{Za) = Xa, relating to the spectral parameters as Ai = cA, A2 = c/i 
and construct the corresponding Lax operator from ( |4.22| ). In the simplest case (see 
(|6.17|) ), such color parameters induce interaction with an external magnetic field. 



4.5.2 Generalized lattice sine-Gordon model 

The Lattice sine-Gordon model obtained above can be generalized (^^LSG) to include 
nontrivial parameters 6 and Z„, retaining its quantum integrability. The correspond- 
ing Lax operator can be constructed by a 6' twist of L^'*^^(A) ( |4.10| ) related to the 
lattice SG model as 

(4.24) 



4.6 Multiparameter and color generalizations 

We can construct integrable multiparameter models by introducing parameters in the 
multicomponent Lax operator ( [4.16|) through twisting transformation. The concerned 
universal 7?.-matrix belongs to G Uq{gl{n + 1)) ® Uq{gl{n + 1)) and may be obtained 
using the twisting operator ( p.99| ) as 

g{(i)i,, V,) = exp j -Uy, v^a^\Hk®Z - Zmk) - 2 E <t^i3Hi®H^ j (4.25) 

by specializing M = l.qu = ,Ho, = Z, aij = 26ij - 6ij+i - Sij^i. 

Here the (f)ij,Vj are external parameters along with a central operator Z respon- 
sible for the color parameters. Note that, nontrivial (pij can be introduced starting 
from n > 2 only, i.e. at least from gl{3), while the parameter vi = 6 can appear 
even in gl{2), as considered above. Remarkably for reductive type algebras some 
additional nontrivial parameters might arise due to the presence of central opera- 
tors Z, since one can consider Z to have different eigenvalues in different spaces, i.e. 
{tt^ ® I)Z = \I ® I, while (/ ® vr^)^ = Jil ® I . Parameters A, /i may be taken as the 
color parameters. 

Therefore, we may construct the multiparameter Lax operators by twisting ( 4.16| ) 
through the operator 

n+l 

(ttx ® / ® /)^?i3 = Gi3(A, Vj) = J2 nicf)ij)WkiX, Vj)Ekk (4.26) 

k=l 
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where 

(4.27) 

for p = 1, ■ ■ ■ , n with = = 0, while iyn+i(A) = V ^^-^ . 

The Lax operator of multicomponent multiparameter quantum integrable models 
thus can be generated from the N x N Lax operator ( |4.16| ) using ( |4.26| ) 

as 

Lis{^,Kvj,(j)ij) = Gis{X,(f)ij,Vj)L['^'"^\^)Gis{X,(j)ij,Vj) (4.28) 

where 4>ij,Vj are additional parameters apart from the color parameter A. As men- 
tioned, the contribution from Tk{4>ij) introducing parameters (pij starts for the higher 
rank algebras with n > 2, while parameters Vj through Wp can enter even for 



4.6.1 Multiparameter multimode g-bosonic model 

Multiparameter generalization of the multimode g-bosonic Lax operator 
constructed above, can be obtained by transformation like ( 4.28 ) with twisting oper- 
ators obtained from (0|g]23) by using the g-bosonic realization ([4.18|). 



4.6.2 Multiparameter Toda field model 

Performing twisting (|4.28| ) on the discretized Toda field Lax operator iji^odafteid) (^^\^ 



1^] obtained above, it can be generalized to include parameters X,(j)ij,Vj, i,j = 
1,2, ... ,n along with the constant inhomogeneous field Z„. For this one has to make 
only the replacement through bosonic operators Hi = Ui in the expressions ([4.27|) . 
Note that n = 1 case recovers the generalized lattice SG model. 

It should however be remarked that, though such multiparameter quantum models 
are exactly integrable at the discrete level, whether they lead to any sensible field 
models at the continuum limit is not yet clear. 



4.6.3 Multiparameter trigonometric _R(A)-matrix 

The multiparameter models are associated with the multiparameter generaliza- 
tion of the trigonometric R{\ — /i)-matrix, which can be obtained by twisting of 
( [4.15| ): i?i2 = Gi2-Ri2(A — yu)G'i2 wherc G12 = GuiX, fi, (pij, ,Vj)) is obtained as 
(/ ® 7!')Gi2i\, (pij,Vj) = Gi2{\, fi, 4>ij,Vj) from G12 ( |4.26| ). In the explicit form the 
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i?-matrix is given by 

n+l n+1 

R (A - /i, A, /2, (f)ij, Vj) = sin(A - fi + a)Y^ glk^kk ® Ekk + sin(A - /j,) J^idlj^kj 

k=l k>j 

Ejj ® Ekk + g%^kjEkk ® Ejj) + sin a ^ gijgjiEij ® Eji (4.29) 
with the notation 

^^^■^/|(A,K})' /'(^'K}) = ^^ ^ (4-30) 



Note that the ^-deformed R ( [4.21|) can be recovered from the general ( [4.29|) as a 
particular case. 

Interpreting the color parameters as spectral parameters: X = cX, fi = cfi, one 
can get from ( 4.29| ) a i?-matrix with two sets of spectral parameters with nonadditive 



dependence. Recall that in standard integrable models the known quantum R{X — fi)- 
matrices (see (|2.17| )) depend additively on spectral parameters. 

For the special colorfree choice /a;(A) = 1 and /^(/i) = Xk, ( [4.29| ) reduces to 

n+1 n+1 

RiX) = sin(A + a)Y^ X^Ekk ® Ekk + sin A {^kjXlEjj ® Ekk 

k=l k>j 

+ ^kjX]Ekk ® %) + sin a 5^ XiXj{E,j ® E^, + Ej, ® Eij), (4.31) 

i<j 



recovering the well known Perk-Schultz model ||9^ . Thus it shows a remarkable con- 



nection between an integrable statistical model constructed in a traditional way and 
the representation of the universal 7?.-matrix constructed through twisting transfor- 
mation. The additional parameters = ±1 appearing in the diagonal elements of 
the Perk-Schultz model can be generated if one starts from a universal 7^-matrix 
related to noncompact groups (nonstandard i?-matrix [|142|] ). 

4.7 Models with rational i?rai(A)-matrix 

As explained in sect. 2, at the limit g — 1 both R, L operators are reduced to their 
rational limits. The rational i?-matrices can be obtained from the trigonometric 
( [4.15| ) or ( |4.29| ) by taking the limit 77— i>0inA = ?7A,a = r]h, which amounts to the 
formal replacement of sin(A — fi) X — fi, sin(A — fi + a) X — fi + h etc. As a 
result ( [4.15|) gives 



AT 



Rrat{X - /i) = (A - fi)I + nP, P = E Eij ® Eji (4.32) 
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with P being the permutation operator. 

On the other hand at the rational hmit the quantum algebra generators reduce 
to spin operators, which in turn transforms the L operator ( [4.7D to its rational form 




L(A) = I ^' , , (4.33) 



with s^,s^ being the generators of the standard s/(2) algebra ( |3.71| ). Different real- 
izations of the spin operators generate different integrable models from (|4.33| ) and all 



of them are associated with the same rational i?-matrix (|4.32|) for N = 2, which also 
coincide with ( p.l9| ) obtained earlier as a direct solution of YBE. 



4.7.1 XXX spin-| chain 

Spin- 1 representation through Pauli matrices, rewrites ( [4.33 ) into 



4-)(A) = + (a'al + a+a" + a'a^) (4.34) 

generating the Lax operator of the XXX spin chain listed in ( |2.38| ) and analyzed in 
sect. 6 and 7. Note that L^^^^^ can also be obtained from (|4.^ ) at the rational limit. 

Similarly, the multiparameter isotropic generalized spin models can be obtained 
at the rational limit of (^J^) or (|]2|). 

4.7.2 Nonlinear Schrodinger equation: lattice and field models 

Through bosonic realization of spin operators called Holstein-Primakov transforma- 
tion 

s+ = (A(2s- A^V))^, s" = 7/^^(A(2s- A^V))i = s - Atp^ ijj , (4.35) 

with [ijjn, ^In] = h^^, one can construct from ( [4.33| ) the Lax operator of the exact 
lattice version of the NLS model as 

'^)(A) = X + {s-A^i^r.y+ yA(2s-A^t^„)^ta+ + VAV^„(2s-A^t^„)^^" 

(4.36) 

as listed in ( |2.39| ). We return to this model in sect. 6 for the Hamiltonian constructing 
and in sect. 7 for its exact Bethe ansatz solution. Considering the spin parameter 
s = /tA~^ and taking the continuum limit A — 0, we obtain s'^ \/^^\ s~ = 
v^7/;t, As^ ^ 1 with V'n ^ iip{x), leading to ^a^L^,'"'") ^ I + -^A£„i,(A) and 
thus generating the Lax operator of the NLS field model 

£„;^(A) = iXa^ + K^(^t(3,y+ + iIj[x)cx-. (4.37) 



60 



The vector NLS model can be constructed from the rational limit of ( [4 .161) with 
proper bosonization. Another possibility of the vector NLS will be considered in the 
next section. 



4.8 Models with twisted rational i?rai(A)-matrix 

Note that the twisting operator F{9, Zn) involved in transforming the i?-matrix and 
the Lax operator remains the same at the rational limit, since it does not contain the 
deformation parameter rj. The rational i?-matrix ( [4.32| ) therefore twisted into 

R^^^tiX- IJ,e) = {X- lj)T + hP, T = diag{l,e-'^,e'\l) (4.38) 

while the Lax operator ([4.33|) turns into 



It is clear that the spin-i representation will lead to the Lax operator of a ^ deformed 
XXX chain given by ( [4.38| ) and corresponds to the Hamiltonian 



Hex.x = E ^^m^m+i + e'V+a„+i + e (4.40) 



m 



obtained also from (|4.23|) at the isotropic limit. 



On the other hand the Holstein-Primakov transformation (|4.35|) would lead ( 4.39| ) 



to a ^^-deformed lattice NLS model (6'-LNLS). Such 6'-LNLS would represent an exact 
quantum integrable discrete bosonic model associated with the twisted -R^^^-matrix 
( [4.38| ). Twisting transformations on multicomponent systems likewise would generate 
multiparameter models. However it is not clear, whether proper continuum limit 
exists for such discrete systems leading to any new field models. 

Note that the systematic constructions described here have partially answered the 
points raised in the beginning of this section. We shall elaborate on this issue after 
the constructions in the next section. 

5 Construction of models from extended trigono- 
metric Sklyanin algebra 

A systematic scheme for constructing integrable ultralocal models starting from the 
universal 7^-matrix is formulated in the previous section. We detect however its 
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shortcomings by noticing that, following this approach we could discover not all 
the ultralocal models listed in sect. 2 , but only a few of them, which constitute 
certain symmetric class. Such models include spin models, sine-Gordon, Toda field 
models and the nonlinear Schrodinger equation. The reason is simply the fact that, 
since we are confined to the universal 7^-matrix belonging to Uq{gl{n + 1)), we could 
generate only those models which can be obtained as realizations of the quantum 
algebra, i.e. the models with Lax operators having certain symmetric form. On the 
other hand, as seen clearly from the list in sect. 2, there are models with diverse 
nonsymmetric structures. This demands therefore that the present scheme should be 
modified suitably to incorporate also the generation of these models. It seems that 



one can do so by proposing a more general type of deformed algebra but at the 
cost of loosing some generality. 

It is found that that at the representation ( |4.2| ), when the universal YBE ( |4.1| ) 
reduces to the relations ( p. 761) , it acquires some extra freedom allowing more general 



choice for Lj^'' than ( p.77| ) obtained as a direct realization of 7^-matrix related to the 
known quantum algebra. To identify such general algebraic structures let us choose 
in upper/lower-triangular form 



L^-' = I 1 , L^-^ ={ " _ 1 , (5.1) 





with yet unspecified operators r. Note however that = are unchanged from 
(p.l5| ) and keeps track of the TZ from which it has been originated through repre- 



sentation ( [4.21 ). Interestingly, Using L^"^^ operators as ( p.lj ) along with ( ^.15] ) one 
extracts from the FRT relations ( p.76|) a new type of quadratic algebra 



[ri2, Tai] = -{q-q ^) (ri+Tg - 

^t^ij = l^^^ij^i^ . ^t^ji = l^^'TjiT^ , (5.2) 
for i,j = (1, 2); q = e*". The related Casimir operators may be given by 



= ) ^2 = T^T^ , D3 = T^T^ , L'4 = 2 COSa (t^T2 + Ti r2+) - [ri2, T2l]^ 

(5.3) 

To study the properties of this algebra one should find the coproduct A, antipode S 
and the counit e related to it. Using the formulas like (|3.79|) it results 

= ® A(r2i) = r+ ® T21 + T21 t+ , A(ri2) = T2 ® m + m ® r+, (5.4) 

^(^21) = -(r+)-V2i(r+)-\ S{n2) = -(r2-)-Vi2(rr)-\ Sirt) = (rt)-' (5.5) 
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and e(rj^) = 1, e(rjj) = 0. This establishes the essential Hopf algebra properties of 
(|5.2|) , though its quasi-triangularity i.e. the existence of universal 7^-matrix is still 
an open problem. 

For relating ( |5.2|) with known algebras, one may check that it maps into (|3.66| ) 



for the symmetric reduction 



+ 7^— r — ^, ri2 = 2z sin aS*-, r2i = 2i sin a;S'4 



2 cos f 2i sin ^ 

along with its Casimir operators. Therefore the quadratic algebra ( |5.2| ) is an exten- 
sion of the trigonometric Sklyanin algebra (ETSA) It is related also through 
symmetric reduction ( |5.6| ) and the mapping ( |3.68| ) to the standard quantum algebra 
(p.69|) with generators (/, s^, S*^) as 



= -rt = ±^'^q^'\ , m = rai = nS~ , n = {q - q'^) (5.7) 

with q = e* 



5.1 Ancestor model with trigonometric Rtrig{X)-ina.triyi 

Note that compared to the previous procedure, the L^"^^ operators have been 
changed here, while the matrices remain untouched. Therefore the above Yang- 
Baxterization procedure, valid for i?-matrices with Hecke algebra must also go 
through in the present scheme. As a result one arrives at the same Rtrig{\ — /i) 
matrix ([4.6|) , while the related L{\) operator takes the form 



with the upper and lower triangular parts coming from L"*" and L~ are combined 
through the spectral parameter ^ = e^'^. Since by construction this L{X) operator 
along with the i?-matrix solution ([4.6|) satisfy the QYBE ( p.9| ), one can take (|5.8| ) as 
the Lax operator of some generalized ultralocal integrable model. Such a quantum 
model would involve abstract operators {r} satisfying the g-deformed algebra (|5.2| ) 
and would be associated with the trigonometric Rtrig{X — ij) matrix. For physical 
applications, (|5^ ) can be looked as an ancestor model, which under different reduc- 
tions and suitable realizations through bosonic, g-bosonic or spin operators can lead 



63 



to various ultralocal descendant models. All such descendant models would naturally 
be quantum integrable, satisfy the QYBE's (P .91^.111) and share the same quantum 



Rtrig{\)-Taia.tYix inherited from the ancestor model. 



Among different reductions of the ETSA (|5.2| ), there is a special symmetric one 

(5.9) 



ri2 



^2 



which corresponds to (|5.7|) and links it to the standard Uq{sl{2)) . Such realization 
in quantum algebra generators reproduces clearly from (|5.8| ) the symmetric ancestor 
model (|4.7| ) and generates the important class of integrable models with the underly- 
ing Uq{sl{2)) symmetry. Recall that such models constructed in the previous section 
from ( |4.7| ) and its rational limit ( [4.33| ) include XX Z and XXX spin-| chains, SG 
and NLS models etc. 



5.2 Integrable reductions of ancestor model 

Note that though ETSA can be mapped into the quantum algebra as (|5.7| ) for invert- 
ible nonsingular generators, it allows in general richer possible realizations through 
nonsymmetric reductions. Consequently, much wider class of models can be gener- 
ated, which are not reachable from the quantum algebra without undergoing limiting 
procedures. Some of such interesting cases are listed below, though in no way they 
exhaust all possibilities. 

Let us consider a particular nonsymmetric reduction 



2 



0, = = c( r^^ ) \ = (r2i)"'', c = const. (5.10) 



which corresponds to the quadratic algebra 
[ri2,T-2i] = -2isina(r+)2 , t+Tu = e*"ri2r+ , r+r2i = e-*"r2ir+ , (5.11) 



reducible from ( |5.2| ). Through bosonic and g-bosonic realizations of this algebra we 
construct below nontrivial models. 

5.2.1 Liouville lattice and field models 



Consider realization of (|5.11|) with = in canonical operators having [u,p\ = ^ 
in the form 



r+ = A e*"", T2i = e-'^Pf{u), = f{u)e'^P, f{u) = (l + e^^^^^+i^^)^ 

(5.12) 
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This results from 

4"™)(0=| fM e^-^ \ where /K) = [1 + AV(")]^ 

(5.13) 

giving Lax operator of the exact lattice version of the quantum Liouville model (LLM) 
listed in ( p.32| ). At the continuum limit: A 0, when {un,Pn) — ^ (u^x) , p{x)) with 
[u{x),p{y)] = h6{x — y) and e^*^^" I + ±iAp{x), f{un) I, the lattice version 
(|5.13|) would lead to L^"™^(^)cr^ — / + A£/m(0 yielding the Lax operator 

CimiO = -P^' + (^^+ + e^")e''"" (5.14) 



of the well known Liouville field model (LM) listed in ( |2.26| ). Both the lattice as well 



as the field Liouville models have the same Rtrig{X) as the ancestor model ( |5.8| ). 

Note that for the choice c = — A^ in reduction ( |5.10|) one gets yet another simple 
and similar realization 

r+ = Ae-'^"^ , T2i = e'^'P , = f\u) e'^^^ (5.15) 



with f[u) being as defined in (|5.12| ). This leads to (after multiplying by from the 
left) 

representing the Lax operator of a nontrivial spectral parameter dependent lattice 



Liouville model, solved recently using algebraic Bethe ansatz in , where this model 
has been derived directly from the sine-Gordon model invoking some nontrivial lim- 
iting procedures. 

We consider next partially symmetric reduction of ( |5.2| ) 

^i" = -c (r,-)-' = (r+)-^ = r,, r^i = A, (5.17) 

giving the algebra 

[r2i,Tu] = -2ic sin a (r+2 + (r+)-2) . (5.18) 
5.2.2 g-oscillator model 



Note that though for c 7^ ( ^.17] ) can be mapped into the quantum algebra (|3.69|) , 
in general it gives more freedom of choice and for c = ^ through realization 

rt = ^r, =g-^ r+ = - ^rf = ru = -^A, r,, = kA^ (5.19) 
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where k = (| A sin 2a) 2 and A, are g-oscillators or g-bosons with the commutation 
relations we can construct an integrable lattice g-oscillator model 

In-^" - ^ oA^n+l ^ /It 



listed in ( ^^M 



Another such realization may be given as reduction ( |5.1(]|) by 

r+ = = AgiN), = 9iN)A^, (5.21) 

where g'^{N) = {q — q~^) ■ The Lax operators /.(^^'oson) ^j^^ corresponding 
quantum integrable lattice models involving g-bosons can be obtained by inserting 
( |5.21| ) in the ancestor operator (|5.8|). 



5.2.3 Derivative NLS: lattice and field models 



Curiously the g-oscillator model ( p.20| ) can be interpreted as the lattice DNLS model. 



To see this we express g-oscillators through bosonic operators with ['i/'„, = ^5nm 
as 

An = (f )^ ^"\/^' = J^i^n (5.22) 



and inserting in ( p.2C ) obtain the Lax operator L^'^^^^\/\^^) of an integrable bosonic 



model, which satisfies the QYBE in all orders of A with the trigonometric i?-matrix 
(|4.6|) . To identify L^''^"'''^(A, ^) as the exact lattice regularized version of the DNLS 
field model (LDNLS), we take the limit limA-.o ^^(^''^"'"^(A, = / + A£f' "(A) by 
using 5="=^" — > / ± rjiplipn, f^An — >■ /^ipn etc. with q = e^^ and demand it to satisfy 
the QYBE up to order A. However we face an initial difficulty due to the fact, that 
such naive limit does not satisfy the QYBE up to the required order and the coupling 
constants needs to be regularized as 



This gives 



1 + iK^h = e 1 — = e*", h=—sma. 



and defining the field as ^ipn = Ix"~^^ '^{^)(^^ with the commutator [ilj{x),ijj'^{y)] = 
— sin a6{x — y), we finally arrive from (|5.23|) at the Lax operator 

^dnis(^^^ = I (^-^^V + '^\x)ij{x)K + ^{ij\x)a+ + ij{x)a-)^ , K = diag{K_, -k+) 

(5.24) 
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of the DNLS field model listed in (|2.27|) . The model is solved exactly through Bethe 
ansatz in sect. 7. The explicit forms of the conserved quantities including the Hamil- 
tonian of the model if^*^^^ are 

r+ao r+oo 

Co = dx ^ Ci = Ai ip'^4'x dx , 



i^'"'^ = IC2= I ' V^Vx. - 22c(V^V)V^Vx ) dx , (5.25) 



etc. 

It is remarkable that, while the standard derivative NLS is nonultralocal |^ 
and could not be solved through QIST, the present DNLS model generated from the 
underlying quantized algebra ( ^.21) is ultralocal and Bethe ansatz solvable. 

5.2.4 Massive Thirring model 

Notice an important symmetry of the algebra ( [5. 21) , that the relations remain un- 
changed under the replacement 



"l"' ''"l ) ^ (''"2 5 ''~2 



(5.26) 



Using this property we can have another realization like ( p. 191 ) in an independent q- 



oscillator mode {A2, A^) giving another Lax operator similar to L^"'^^^^^ O). 

Since the product of mutually commuting L-operator solutions of QYBE is itself a 
solution, = Lv>sci(i) ]^qosci(2) j^^g^ solution of QYBE and consequently 

would represent the Lax operator of an integrable model with two independent q- 
oscillator modes ([0) 

Lr--(o = { ] (5.27) 

with the elements 

Ln - q -^q +—q - k A,A, 

U2 = I ( - ^ 4g^^^^) + < (^U-^^ - ^ g^^^^4) (5.28) 

and their complex conjugates. It is exciting to note that, when mapped through 
bosons as (|5.22|) , the Lax operator (|5.27| ) represents an exact lattice version of the 



massive Thirring model (LMTM) (bosonic), which satisfies QYBE in all orders of A 
with trigonometric i?(A)-matrix. 
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At the continuum limit the Lax operator of the corresponding field model reads 

(l-mtm = % (^^(^ - + «Pl - «P2 + (eV^I + + (eV^l + ^^2)a- j . (5.29) 

with pa = {ipltpa), ^ = diag{K^, —K-^-), and k = diag{Kj^, —k^). ( |5.29| ) coincides 
with the two component relativistic field model ||99| and also generates the well known 
massive Thirring model (bosonic) (MTM) with the Lagrangian 



where ip = (z/'i, ^"2), [V'a(x).V'|(?/)] = h6abS{x - y). 
Let us consider one more nonsymmetric reduction 



rf = r^- = 0, r+ = {t+)-\ [t^i, t^2] = (5.30) 

and its complementary one, obtained through the symmetry ( |5.26| ), both of which 
answer to the simple algebra 

[ri2,r2i] = 0, r+ri2 = e^'W , t+t^i = e-^"r2ir+ . (5.31) 

This allows generation of various other models through different realizations. For 
example, a lattice version of the light-cone sine-Gordon model (LCSG) studied in 
|100|| is obtained for 

T+ = e-'^P, ri2 = 4i = ^Ae'"". (5.32) 
On the other hand, the realization 

ri+ = -in)-' = -e''^, r2i = -{ruy'v' = (5.33) 
yields a discrete-time or relativistic Toda chain (RTC) with the Lax operator 

The Hamiltonian and momentum of the model are given by 

H = ^2 (cosh 2ripi + -rf cosh 77(^4 + pj+i)e'^'~'''+^ 

i 

P = Y1 (si^'^ + rf sinh ri{pi + pi+i)e^'"^"+i) . (5.35) 

i 

derivation of which from the Lax operator is shown in sect. 6. A family of such 
relativistic quantum Toda chains are constructed below. 
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5.2.5 Discrete-time quantum Toda like model 

We consider now another nonsymmetric realization of ( |5.2| ) through g-oscillators at 
lattice site i as 

= T-^e^^(^»+-«), ri = ^6^"^% = f3^At, n,, = a.cosa Ar, (5.36) 
z sm a 27 

where a = hT],h,'y,'r] are homogeneous and ai,(3i inhomogeneous constant param- 
eters, are g-oscillators, a bit different from and given by the commutation 
relations 

[iVfc, A^] = ±i5kih At A-^Ai - e^^'^'^AiAl = i^ufi e^'^^^Ar+o.-.;.) ^^^g^^ 

Through canonical variables these operators may be expressed as 

N = pq, A- = qf{N), A+ = f{N)p, (5.38) 

with P{N) = + UJ- zh], -[u- ^n],) and [x], = 

Inserting (|5.36| ) in ( p.S] ) and assuming C, = e^^ one derives the Lax operator 



ai cos hrj A~ cosh (77 (A^ 



'(A)= - :t .,z: ... , (5.39) 




of a g-oscillator model associated with Rtrig{X). Remarkably, this quantum integrable 
model with mapping ( ^.38| ) may be interpreted as the discrete-time as well as the 
quantum generalization of a Toda chain like model (DTTCL) proposed in and 
listed in ( p.36|) . At the time-continuous limit: 77 ^ 0, as seen from (|5.38|) , = 



pq, A —>■ q, v4+ =^ p , one recovers the Lax operator (|5.58|) , which is a quantum 



version of [p^] . 

We stress again that all the above models obtained as descendants from the same 
ancestor model ( [5.8[ ) are quantum integrable by construction, since they satisfy the 
QYBE as different exact realizations of ETSA (|5.2|) . Moreover they are associated 
with the same quantum i?tric/(A)-matrix ( ^4.6| ) of trigonometric type, inherited from 
the ancestor model. The family of integrable models is specified thus by the associated 
i?-matrix and the classification of integrable systems could therefore be based on the 
classification of the quantum i?-matrix assigned to them. Search for any new family 
consequently should start with a new solution of YBE ( ^131) . For this reason we 
switch next to the twisted i?-matrix solution related to gl{2) for constructing other 
kinds of models. 
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5.3 Twisted quadratic algebra and related models 

Twisted A) -matrices as ( [4.21D can be obtained from the original trigono- 

metric -Rjric/(A)-matrix (|4.6| ) by a simple twisting transformation Rf^^g{\) = 
F{9)RtrigiX)F{e), F{e) = e^^('^3^i-i^'^3)^ ^^icli takes a solution of the YBE to 
another solution and introduces an extra parameter 6. 

For constructing the associated Lax operator Lf^^g(,^,r) of the ancestor model, it 
may be taken in the same form ( |5.8| ), while the r operators involved satisfy now a 
6'-deformation of the quadratic algebra ETSA given as 

t ri2T2i - t^^T'2ih2 = - 2z sin a (fif^ - f{f:^^ , 

rtnj = t q^fijft , ^t^ji = Q^t-^^jif^ (5.40) 

where g = e*", t = e*^. Note that the coproduct structure for this algebra remains the 
same as ( |5.4| ). Redefining generators it is however possible to shift the ^-dependence 
from the algebra to its coalgebra as done in the general treatment of sect. 3. It is 
evident that at ^ = ±a, i.e. when the additional t-deformation coincides with the 
original g-deformation, the algebra ( |5.40 ) is much simplified. Interestingly, most of 



the physically significant models, as we see below, are generated at this particular 
resonating situation. 

We may consider as before different reductions of the 6'-deformed ETSA for finding 
relevant realizations. Notice that the symmetric reduction 

f+ = -f2- = -ze*("+^)^', ff = -f+ = ie-"^^-'^'' 

is related to the multiparameter quantum algebra known as Uq^p{gl{2))) we have 
already considered. 

Through bosonic realization of (|5.41|) we construct the ^-deformed lattice sine- 
Gordon model, derived previously in sect. 4 starting from the universal 7^-matrix 
related to deformed 5'/(2). Similarly realizing ( ^.41|) in Pauli matrices we get the 
Wu-McCoy model (p30|) . 



On the other hand, the nonsymmetric reductions ( ^.10| - ^730| ) may lead to the 9- 



deformation of all other integrable lattice models like LDNLS, LMTM, LCSG, LLM, 
g-oscillator models etc. obtained earlier having the same quantum i?j^j^(A)-matrix 
(pTD. 
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5.3.1 Ablowitz-Ladik Model 



Such a 6'-deformation of the g-oscillator reahzation ( p.21| ) gives 

^_,(e+«)7v+|^ ^21 = A~g{N), n2 = ~g{N)A\ (5.42) 



^1 = ^2 



where g'^{N) = ne**-""^^^ , = — q^^) ■ Observe that for the special value 
6 = —a, ( ^.42| ) is simplified to fi = fg" = 1, f2i = h, fi2 = 6^ , where h = Aq^il^ 
and = g'^f]^^^ is another form of g-oscillator with the algebraic relations 
q'^bb^ — b^b = — 1 . Remarkably, this deformed oscillator model represented by 

LniO = ' ( V' f) f^'"' " " ^"^"^ ^^-^^^ 

with h = 1 — is the Lax operator of the well known Ablowitz-Ladik model (ALM) 
1^], proposed as a discretization of the NLS equation years before the discovery of 



g-oscillators. 



5.3.2 Family of relativistic quantum Toda chains 



Exploiting 6'-deformed quadratic algebra (|5.4CI| ) one can construct a family of quantum 



integrable discrete-time or relativistic generalizations of the periodic Toda chain, 
apart from (|5.34|) obtained above. These models are canonically related ||101|| to 



the relativistic Toda chain of Ruijsenaars ||102|| and can be generated from a single 
ancestor model at different realizations of the underlying algebra [^] . 

Considering reduction complementary to (|5.30|) with = 0, one simplifies (|5.4CI| ) 

to 

ri2r2i = e'''%^f^2 , rfh, = e'^^'^^'^uf^ , ^i^^2i = e^^^'^-'h,,f^ (5.44) 

involving only Weyl type relations. This observation facilitates to find the realization 

f+ = e^'^+'^P, ff = -e-(''"^)P, fsi = r?e""P+'', fi2 = -r^e^^+^^^P"" (5.45) 

where a = rih, 6 = eh. We remark that certain operator ordering should be main- 
tained in all such expressions and consider the periodic boundary condition. ( |5.45| ) 



gives from (|5.8| ) the explicit form for the Lax operator 
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depending on the additional parameters e and c and associated with the twisted 
quantum i?(6', A)-matrix ( [4.21| ). The Hamiltonian H and the momentum P of the 
generahzed Toda system may be given by 



H = ^(cosh 2r]pi + if cosh rj{pi + Pi+i)e' 

i 

P = ^2 (sinh 2ripi + sinh r]{pi + pj+i)e' 



t{l+c)(p^+i-pi) + {qi-qi+i)\ 



-e(l+c)(pi-pi+i) + (qi-g,+i) 



(5.47) 



This integrable quantum model is able to generate the whole family of relativistic 
Toda chains for different choices of parameters c, Cq with e = CqT]. Moreover, since 
(|4.21 ) is independent of c, models with different values of c would share the same 
quantum A)-matrix ([4.21|) . 



Some members of this quantum integrable relativistic Toda chain family deserve 
special attention. The symmetric reduction with 6' = 0, corresponding to c = Cq = 
has already been found with the most simplification as ( |5.34 ). The related Hamilto- 
nian and momentum are expressed in symmetric form (|5.35|) . Interestingly, even for 
nontrivial 9 or e, one can achieve the same symmetric conserved quantities (|5.35| ) by 
choosing c = —1, as seen clearly from ( |5.47|) . However the corresponding i?-matrix 
will be different with a more involved Lax operator. 

For Co = 1 and c = 0, (|5.45|) reduces to 



-1 T21 



(5.48) 



with the Lax operator easily obtainable from ( ^.46| ). The conserved quantities are 
similarly reduced from (|5.47|) simply as 



E 



'^pili-Qi+i) 



(5.49) 



with H 



|(/+ + / ) and P = |(/^ — / ). This model corresponds to the quantum 



version of the discrete-time Toda chain | |101| | and the associated quantum Rtrig{o:, A)- 
matrix coincides with (|4.21|) for 6 = a. 

Another quantum generalization of ||101|| is achieved for cq = 1 but c 7^ 0, for 
which the Lax operator and the conserved quantities are obtained from the general 
forms ( |5.46| , p74^ ). This model though different from the above due to different c, 
shares the same quantum R{a, A)-matrix. 

Finally, one observes that the Hamiltonian and momentum (|5.47|) of the generat- 
ing model is transformed exactly to ( ^.49|) under canonical transformation {p, q) 
{P, Q) as 

p = P, q + {r]-t{l + c))p = Q (5.50) 



72 



Therefore, since the other models are obtained from the generating model at different 
parameter choices, all of them are canonically equivalent to (|5.49|) , which in turn 
has been shown ||1(J1[| to be equivalent to the Ruijsenaars model. Therefore the 



whole family of discrete-time or relativistic Toda chains discussed here are canonically 



equivalent to the relativistic generalizations of the Toda chain of ||102| . 

In sect. 7 it will be shown, how to achieve the separation of variables in the eigen- 
value problem for such quantum Toda chains by applying functional Bethe ansatz. 
Thus from the ancestor models through suitable realizations, one can generate sys- 
tematically the Lax operators of various quantum integrable ultralocal models of the 
trigonometric class. 

5.4 Models with rational i^^at (A)-matrix 

We have seen already that at q —>■ 1 or 1] ^ 0, the trigonometric Rtrig{)^)-'oia.thx 
reduces to its rational form (|4.32|) . We are interested to examine here how at this 
rational limit the underlying algebraic structures ( |5.2| ) are changed and how ( p.8| ) is 
reduced at this limit to yield the Lax operator for the ancestor model representing 
the rational class. 

It is not difficult to find the consistent limit of the operators, when the deforming 
parameter r/ — > as 

Using such operators along with the reduction of the spectral parameter ^ —>■ l + ir]X, 
the L operator ( |5.8| ) undergoes a smooth transition and simplifies to 

LiXU=\"°^"''i ^''-l. (5.51) 



The K operators satisfy the algebra, obtained at the g — > 1 limit of (|5.2|) : 

[Kl K^2] = e^KuKl, [K^ K^i] = -e^K^iK}, (5.52) 

with ei = 1,62 = ~1 and K\^K\ serving as Casimir operators. The algebra ( |5.52| ) 
corresponds to a representation of the Yangian, more general than Y{sl{2)) presented 
in sect. 3. 
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Since L[X) operator given in the general form ( ^.51|) satisfies QYBE and associated 
with the rational Rrat{X — fi) = X — n + hP matrix, it may be taken as the ancestor 
model for generating integrable models belonging to the rational class. 

Realization through spin operators: 



12 



K. 



21 



(5.53) 



recovers the sl{2) algebra (|3.71 ) corresponding to ( p.88 ). The ancestor Lax operator 
(|5.51| ) for this symmetric reduction goes to ( [4.33|) considered above. This produces 
the XXX chain ( [4.34|) for the spin-| representation in Pauli matrices and the NLS 
lattice and field models for the Holstein-Primakov transformation ( [4.35| ). 

The full algebra ( |5.52D on the other hand allows more freedom for constructing 
models with nonsymmetric reductions. For example if we consider K\ = i, K2 = 
Kl = 0, (|53^ ) simplifies to 

[Kl Ku] = tKu, [Kl K21] = -iK2i [Ku, K21] = 

and admits the bosonic realization 

This yields from (|5.51|) the Lax operator of the well known Toda chain (TC) 

Pn- X 



(5.54) 



(5.55) 



-qn 







(5.56) 



listed in ( p. 411 ). We return to this model in sect. 6 and 7 for highlighting other 
aspects. 

5.4.1 Quantum Toda chain like model 

Another nonsymmetric realization through canonical variables Pi,qi, [(li,Pj] = Sij : 



derives from ( ^.51|) the Lax operator 

U{X)=\ ^ + ^^^* + -^ 



7 



K12 = aq, K21 = /Sap (5.57) 



(3iPi 



7 



(5.58) 



with constant parameters 7, Pi, Ui. ( |5.58|) thus constructed represents the quantum 
generalization of a Toda like classical model (TCL) proposed in . The model can 
also be obtained directly from ( ^.36|) at 77 — > 0. 
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5.4.2 Simple lattice NLS model 

The lattice NLS generated in ( [4.36[ ) tlirough symmetric reduction from ( p.51| ) has 
rather involved structure and the corresponding Hamiltonian of the model, as will 
be shown in the next section, is also much complicated. Notably, the algebra ( 5.52| ) 
allows a nonsymmetric reduction realized as 

Kl = KA^(f)ij + 1, Kl = -A, K° = 1, = 0, 7^12 = lA^ij, K21 = -iAy^(f) 

(5.59) 

where the operators ilj,(f) obey the canonical commutation relation [ip{n),(f){m)] = 



■^Snm- Therefore, (|5.51| ) generates a AKNS type integrable model (SLNLS) given by 
the Lax operator 

iy/KAip{n) 1 



Through an irrelevant scaling L L = j^L (see ( p!.4(JD , it is convenient to con 



struct the explicit conserved quantities of the model in a much simpler way as 
will be shown in the next section. Though at the discrete level the operators ijj 
and are independent from each other, at the continuum limit : L„(A)''*'"'*^ 
1 -iA/:"'^(x,A) + 0(A2) , we get 

/:"'^(x, A) = ^(1 + (r^)\ + v^0a+ - v^t/;^- (5.61) 

Due to the restoration of the symmetry, one can put = ■j/'Mn the corresponding 
conserved quantities, which recovers the standard NLS field model ||52| . 

All the above models constructed as descendants of ( |5.51| ), share the same rational 
-Rrat (A)-matrix inherited from the ancestor model. 

5.4.3 Models with twisted rational i?^^^(A)-matrix 



We consider ^-deformed rational i?^„j(A)-matrix given by ( [4.38|) , which induces 
through QYBE a deformation of algebra ( p.52| ) to 



K^Ku = tKuiK^eiKl), K^K^i = r^K^iiK^, - e^K]), (5.62) 

with t = e*^ and i = 1,2 . Here operators K^,K} form a commuting set, though 
noticeably, K^s are no longer central elements of the algebra. ( ^.51|) with such K 
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operators would represent the Lax operator Lf^^(A) of another ancestor model gen- 
erating rational class of quantum integrable 9 deformed models, associated with the 
twisted R^^^{X). For example, symmetric reduction with realization in spin operators 

Kl = Kl = f ^ = -K^ = sV^ Ki2 = 7^21 = t''+^ s- (5.63) 

reproduces the Lax operator ( [4.39| ) and recovers the 6'-deformed XXX and the lattice 
NLS models constructed before. Through nonsymmetric reductions of ( ^.51| ) however 
we can generate wider class of twisted models. 

5.4.4 ^^-deformed Toda chain 

The bosonic realization 

= Kl = 0, Kl = ^e^^ = pe'^, = e-("+i^)+^^ K21 = e("+t^)+''P (5.64) 



satisfying ( ^.62| ) rewrites (|5.51| ) as the Lax operator Lf^ of a quantum integrable 



deformed Toda chain {6-TC). 

5.4.5 Tamm-Dancoff g-bosonic model 

By deforming the lattice NLS model Lji(A)^'''"''*'' found above, through the parameter 
q = e'^ one gets the Lax operator 

y i^b f{N) J 

of a quantum integrable system (TDqB involving Tamm-Dancoff type g-bosonic op- 
erator ||143|] satisfying [b, N] = b, [c, N] = —c, be — qcb = q^ . The g-bosonic 



models, as we have seen, are usually related to the quantum algebras and associated 
with trigonometric Rtrig{X)-v[va,tiix. However the present TDqB model is connected 
with the rational matrix (|4.21|) with 9 = a. Such g-bosons can be mapped easily to 
the NLS fields as b = f{N)ip, c = f{N)(j) with [ip,(p] = 1 leading to the rela- 
tion L'^{b,c,N) = / (X)L„ ('?/', 0)*^^'"''*^ Such 6'-deformed models though correspond 
to the rational class with g = 1 algebra, they exhibit properties resembling quantum 
algebra. 

The ultralocal quantum integrable models generated systematically following the 
above procedure, from the ancestor models with trigonometric and rational -R(A)- 
matrices as well as with their twisted forms are presented schematically in fig.5.1. 
The abbreviations used to denote the models are as in their list in sect. 2. 
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Thus the following results are achieved through the present scheme matching to 
some extend the original aims and answering the questions raised at the beginning 
of sect. 4. 

i) Ancestor Lax operators of exactly integrable quantum lattice models are con- 
structed in a rather systematic way. Through consistent realizations of these models 
one could generate different families of quantum integrable ultralocal systems repre- 
senting known as well as new lattice and field models. 

ii) The descendant models generated from the same ancestor model share the same 
inherited i?-matrix. And this classifies the models into two broad classes: trigonomet- 
ric and rational, explaining thus the commonness of diverse models pointed above. 

iii) The notion of integrable nonlinearity seems to be determined by concrete realiza- 
tions of the quantized algebras ( |5.2| , p.4q ) (or their g — >• 1 limit ( |5.52| ,thk-alg) , as given 
for example by ([4.9| , |5.12y5.55|) etc. Transition to continuum limit usually distorts such 
exact nonlinearity leaving only some remnants. 

Extension of the scheme presented here to higher dimensions and its application 
for constructing ( (A^ x A^)) Lax operators will be taken up in sect. 9. 



fig.5.1 



Fig. 5.1 Generation of quantum integrable models from respective ancestor 
models with Lax operators Lf, Lt related to the trigonometric Rf, Rt matrices 
and with Lax operators L^, Lr related to the rational R^, Rr matrices obtained at 
particular limits. The subscript t stands for the trigonometric and r for the rational 
class. 

6 Hamiltonian and other conserved quantities 
from Lax operator 

We have seen how the representative Lax operators belonging to an integrable class, 
determined by the associated i?-matrix can be constructed in a systematic way. How- 
ever, since the Hamiltonian of the model is an important physical object, we focus here 
on its construction along with other conserved quantities, starting from a more ab- 
stract element, namely the Lax operator. Recall that (see sect. 2) the conserved quan- 
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titles {Cn}, n = 1,2,... are generated from the transfer matrix r(A) as the expansion 
coefficients of Inr(A) = Z^^o^"-^" therefore C„ = const. ^ lnr(A)jA=AO, where 
the expansion point A*^ may be 0. The global operator r(A) is constructed from the 
local Lax operators Laj{X), j = 1, . . . N a.s 

r(A) =tr,(L,i(A)...L,^(A)), (6.1) 

where r(A) acts on the total quantum space 7i = ®n=i^ri- Therefore, from the knowl- 
edge of Lax operators it is possible to derive expressions for conserved quantities 
including the Hamiltonian. In general, they will be nonlocal objects and physically 
uninteresting. However some special properties of the Lax operators allow to over- 
come this difficulty and obtain Hamiltonians with localized interactions. The Lax 
operator Laj{X) acts on the product space Va hj, of the common auxiliary space 
Va and the quantum space hj at site j. The models with Va and all hj,j = 1, . . . N 
isomorphic are called fundamental, while the rest with inequivalent quantum and 
auxiliary spaces are known as the nonfundamental models. For fundamental models 
since Va ~ hj, the Lax operators Laj{X) formally coincide with their corresponding 
i?aj(A)-matrices. 

6.1 Fundamental models 

Lax operator Laj{X) of such models acts on V ® V and therefore the auxiliary and 
quantum spaces may be interchanged. To utilize this feature the regularity property 
of the Lax operator Laj{X = 0) = Paj, i.e. i^^^(O) = 5^5^ is demanded, where the 
permutation operator Paj can interchange the spaces as PajLak = Lj^Paj with the 
property = 1 and tra{Paj) = 1- Therefore, for regular models the transfer matrix 
O yields 

r(0) = tra {Pa.PaHl . ■ ■ Pa.^l) = {Pjj+1 ■ ■ ■ Pjj^l)tTa{Pa,) (6.2) 

for any j, using the freedom of cyclic rotation of matrices under the trace. Taking 
derivative with respect to A in (|6.1|) we similarly get 

N N 
t'(0) = tra E {PaA^A^) ■ ■ ■ Paj-l) = E^J+M ' ' ' ~l)tra{Pa,) , (6.3) 

where we have assumed the periodic boundary condition: LaN+j = Laj. Defining 
now H = cCi = c^lnr(A)|A=o = c r'(0)r"^(0), c = const. 
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we construct the Hamiltonian from (|6.2| ) and (|6.3|) as 

N 

H = cj:m,^MPn+i) (6.4) 

with only 2-neighbor interactions, due to cancelation of all other nonlocal factors and 
the property = P. 

Similarly taking higher derivatives of Inr(A) higher conserved quantities Cn,n = 
2,3,... may be generated. Note that the conserved operator Cn involves the interac- 
tion of n + 1 neighbors. 



6.1.1 Hamiltonian construction 

We take the Lax operator coinciding with the i?(A)-matrix in 4 x 4-matrix form 

\ / 1 \ 



I «+(A) 

6+(A) c+(A) 

c_(A) 6_(A) 
V a-(A) ) 



P 



1 

1 



V 



(6.5) 



with more general entries and demand the regularity condition Ljj+i(0) = cPjj+i. 
Multiplying L'^^j^^i^) with Pjj+i and inserting in ( |6.4D we get the explicit expression 
for the Hamiltonian 



N 



H = const.5:6;(0)a>-^i + 6'_(0)a-a+_,i + (/+- J+)a3a3^i 



n=l 



+ (/+ + J+)I + 2/-a^ + J-(a3-4^i) (6.6) 

which is general enough to cover all spin models considered below. Here the param- 
eters are = i(a^(0) ± a'_(0)) and = i(c'^(0) ± c'_(0)). It should be noted that, 
for the models with periodic boundary condition: (Tn+i = <?! , in which we are mostly 
interested, the last term in ( |6.6| ) representing a boundary term simply drops out. On 
the other hand, the Hamiltonian with open boundary condition takes the form 



N 



n 



n=l 



(6.7) 



which acquires an interesting Uq{su{2)) quantum group symmetry. Such models were 



investigated in connection with open chain problems |p.03| , |89 
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6.1.2 XXZ spin-i chain 

The Lax operator for this model is given by the trigonometric i?-matrix with 
regularity condition and corresponds to a^(0) = a'_(0) = cos 77, &'|_(0) = &'_(0) = 1 
and c^O) = c'_(0) = 0. This clearly reduces ( |6.4| ) or similarly ( |6.6| ) to 

N .1 



H^^^ = (^- cos 77(1 + (xjlaj^+i) + + a„ j , (6.8) 

which after neglecting the constant term reproduces exactly the anisotropic XXZ 
Hamiltonian (|2.29| ). 



6.1.3 XXX spin-i chain 

Since in this case the Lax operator may be given by the rational i?-matrix ( |4.32| ) 
having the regularity condition and with the values a',_(0) = a'_(0) = &'_,_(0) = &'_(0) = 
1 and c'_,_(0) = c'_(0) = 0, one gets from (|6.5| ) L^-j_^^(0) = I reducing ( |6.4|) simply to 



N 

H = cY.Pmi- (6-9) 

This is clearly equivalent to the known Hamiltonian of the isotropic spin-i chain 
model listed in ( p.38|) . It can also be obtained from as the degenerate case 

ri = 0: 

^1 

H^xx = 9^n^n+l + + (6-10) 

n ^ 

6.1.4 spin-1 chain models 

A simple way of getting an Integrable spin-1 model is to construct the Lax operator 
as a ^^/(S) related rational i?-matrix. This would lead to (|6.9|) with the permutation 
operator Pjj+\ = J2ab=i ^ib^ba'^^ which when expressed in spin-1 operators S takes 
the form 

TV 

Hspini = const. ((S„S„,+i) + e(S„S„+i)^) (6.11) 

n 

with e = 1. Another more interesting possibility is to use the Lax operator obtained 
through fusion of spin-i Lax operators ||104|| as 



i?ab(A) = (P+ ® P+) n n + ^^(2^ - ^ - ^■))(^a+ ® P^)92sW (6.12) 

j=l k=l 
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with spin s = 1 ||105|| , where is the symmetrizer in the fused spin-s space a and 
5^25 (A) is some multiphcative factor. 

The corresponding Hamiltonian is given in the form ( |ti.ll ) with e = — 1. 



6.1.5 Impurity model 

Simple integrable impurity models may be constructed by replacing the Lax operator 
Laj(A) of the regular model by an inhomogeneous Lam (A — i^) at a particular lattice 
site m, where u corresponds to the impurity at m. We see immediately, that the 
regularity condition no longer holds at m and consequently, the terms involving 



m — 1, m, m + 1 in ( |6.4| ) become nonregular. As a result, instead of ( |6.4| ) one gets the 
Hamiltonian of the impurity model as 

N 

jy^m—l.m 

+ Lm-lm{^)L'^-lm+li^)Pm-lm+lLj^_ij^{h')) (6.13) 

with next-to-nearest neighbor interaction. Using therefore the L-operator of the 
XXX or the XX Z chain one can derive from ( |6.13[ ) the Hamiltonian of the respective 
impurity model. 



6.1.6 Inhomogeneous model 

Inserting different impurities uj at sites j we can get an inhomogeneous spin model, 
where instead of ( |6.1|) one obtains 

r(A) = tra {Lai{X - ui) . . . LaN^X - i^n)) . (6.14) 

Note that at A = due to the loss of regularity at all points: Laj{—i'j) 7^ Paj-, the 
corresponding Hamiltonian becomes nonlocal and as such of no physical interest. 
However interestingly at A = Uj the regularity is restored and one can use the dual 
conditions on the L-operator at 77 = as 

Laj{X - Uj,r]) \\=uj= Paj Laki^j ~ l^kyV) \v=o= ^ (6.15) 

This allows to construct a new set of Hamiltonians Hj = ■^T{uj) \r]=o= Hkj^j Jij^i^j ~ 
i^k,v) \v=o ■ Notice that though this is a bit different definition of the conserved 
quantities, their commutativity also follows from that of r(A) at different A's. 

We can check easily the property (|6.15|) for the L-operators of the XXZ and 
XXX models. The long range Hamiltonian for such inhomogeneous XXX model 
therefore can be constructed as Hj = J2k=ij lyz^^j ' exhibiting similarity with 
the Haldane-Shastry model (see sect. 9). 
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6.1.7 Alternating spin chain 

The idea is to replace the spin-| Lax operators Lg-o-' at every alternate lattice sites 
by a Lax operator L^jg with the same auxiliary space a but different quantum space 
s, corresponding to the spin-s operators. Since at alternating sites the regularity 
condition does not hold, in analogy with the impurity model one gets a Hamiltonian 



with interaction between different spins ranging to next to nearest neighbors ||106|| . 



6.1.8 Twisted and colored models 

Twisted as well as color deformation of Lax operators have been considered in sect. 
4. It is interesting to see their effects in the Hamiltonian of the models. Note that, 
under twisting transformation Rab{X,9) = FabRab{^)Fba ^ the regularity condition is 
preserved as 

Rab{^) = FabPabFba^ = PabPhaF^,^ = Pab 

while the term -L'(O) in (|6.3|) is replaced by 

^..+i^;,+i(o)i^>i, + 2F;^^MFiAAmj+i- (6-16) 

We consider the realization of the twisting operator as 

Faj+i{e, A) = e3((^^+^^H+i-*^-') (6.17) 

with the choice of the color parameter through the spectral parameter A as = cA 
and all other parameters Aj = 0. 

If we consider first the twisting involving 9 parameter only and no A, then due to 
-Fjj+i(0) = we get from ( p.4|) and ( |6.16| ) the expression 



N 

H = Y.F,,+,{e){L',,^MPmi)FiAi{e). 

This leads to a twisted form of the XX Z chain given by the Wu-McCoy model ( |2.3CI| ) 

N -, 

i/-- = const. Y: e^V>-^, + e-V-a+^, + - cosT^a^a^^,. (6.18) 

n=l ^ 

One can derive also the same model directly from ( |6.6|) by considering twisting of the 
parameters 

6^(A) = e^^^6(A), a±(A)=a(A), c± = c. (6.19) 

Note that by transforming = e*"^cr^ the Hamiltonian ( |6.18|) can be reduced to 
the untwisted form (|6.8|) , however with the quasiperiodic boundary condition a^+i = 
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-iNe~ 



Oi- More general type of twisted boundary conditions in integrable spin chains 



have been studied in 107 



On the other hand, for twisting transformation ( |6.17] ) involving spectral pa- 
rameter cA induced through the color parameter as Fjj+i(A) = e^'^^"^^+^, we get 
2-Fjj_,_i(0)Fj^^i(0) = |ccr|^]^. This yields an extra term \cYl!j in the Hamil- 
tonian, signifying the interaction with an external magnetic field |c along the z 
direction. Note that this is consistent with the choice 

h^{X) = e^5^6(A), a±(A) = e^^^a{X), c± = c (6.20) 

which gives /_ = | and can derive directly from (|6.6| ) the Hamiltonian of the model 
as if = + § Eft^f. 

6.1.9 Free fermion models 

It is easy to check, that if we choose the parameters of the L operator in ( |6.5|) as 

a+(A) = a-(A) = cos A, 6+(A) = &-(A) = sin A, c+ = c_ = 1 (6-21) 

satisfying the free fermion condition + 6^ = c^, the regularity condition of the 
corresponding Lax operator /^^(A) is valid and the Hamiltonian ( |6.6|) reduces simply 
to the XY chain 

N 

H.y = cE + ^74+1- (6-22) 

j 

The explicit equivalence of the XY chain with the free fermion model can be shown 
by realizing the Pauli matrices though the anicommuting fermionic operators Cj,c] 
using the Jordan- Wigner transformation 

Ci = (^?---^|-i)^7 (6-23) 

This clearly leads to 

c]cj = rij = i((T| + 1), c]cj+i = cjc]+i = (6.24) 



etc. with CjCj^i = —Cj^iCj and rewrites (|6.22| ) as the free fermion model 



N 

Hjf = -c E Cj^+i + (6.25) 
i 

Another free fermion model can be constructed remarkably from the supersym- 
metry related nonstandard i?-matrix |JL4'2[| consistent with the parametrization 



a+(A) = sin(A + r]), a_(A) = sin(— A + r^), (6.26) 
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with other parameters as in the XXZ model. The regularity condition for this 
model follows due to a+(0) = a_(0) and we notice that, since a+(0) = — a'_(0) 
and c^(0) = c'_(0) = 0, the interacting term I]ct|c7|+i in ( |6.6| ) drops out from the 
Hamiltonian. This leads again to a XY chain placed in a magnetic field given by 
the term cosr]Y,f cr|- The transformation ( |6.23| ) reduces the Hamiltonian clearly to 
a free fermionic model. 



6.1.10 Hubbard model 

A nontrivial example of the regular models is given by the Lax operator of the one- 
dimensional integrable Hubbard model, where the Lax operator does not coincide 
with the i?-matrix and moreover the spectral dependence of the i?-matrix is both of 
the difference as well as the sum. The Lax operator corresponding to the model may 
be given as 

^r(^) = I A) ® I ^)) (6-27) 

where F^{a \ A) and /^^(r | A) are the Lax operators of two XY spin chains involving 
mutually commuting spin-| operators a and r and = j. 

We check first the regularity condition using the property of the XY Lax operator 
ps/(A);,=o = P as L^"''(0) = PS'((T I 0) ® ps'(r | 0) = P(a) ® P(r) = Pnub, where 
we have assumed h{0) = for the function h{X). Therefore following the above 
procedure , apart from the separate contributions of 1^^ related to a and r operators, 
we also get an interacting term h'{0)ajTj in the Hamiltonian coming from t'(0). This 
finally gives the Hubbard Hamiltonian in the form 

Hnub = H^y{a) + H,j.y{T) + Hinticrr) 

= E + ^i^t^i) + + + -^^M) ■ (6-28) 

We can transform Hamiltonian (|6.28| ) into its conventional form ( [r.96D by applying 



the Jordan- Wigner transformation (6.23-6.24) 



6.2 Nonfundamental models 

When the auxiliary space V and the quantum spaces hj = h are not isomorphic, 
the regularity of the Lax operator i^aj(A) does not hold and the construction of the 
Hamiltonian becomes difficult. However there are few methods, that can in principle 
derive the conserved quantities from the Lax operators of such models, which are 
mostly bosonic models. 
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6.2.1 Generalized spin model method 



The first method is based on the idea of considering bosonic models as the generalized 
spin models with arbitrary spin ||108|| . One has to take the fused -Rab(A) matrix (|6.12| ) 
with spin s along with possible gauge transformation and replace the auxiliary spaces 
by the quantum ones. We thus recover the regularity in the quantum space : i?ab(0) = 
Pah and may construct the Hamiltonian of such generalized spin models with nearest 
neighbor interactions. The isotropic and the anisotropic cases with arbitrary spins 
are determined by fusing i?-matrices corresponding to XXX or XX Z spin-| chains 
and belong to two different integrability classes: rational and trigonometric. In the 
first case realizing the spin s, su{2) generators through bosons we construct the 
Hamiltonian of the corresponding bosonic model. Similarly in the second case the 
bosonic realization of the Uq{su{2)) generators is required, which would lead to the 
nonfundamental model of the trigonometric class. 



6.2.2 Projector method 

In a more straightforward approach, known as the projector method [^, one 
starts directly with the Lax operator Laj{X) of the discrete bosonic model with- 
out demanding any regularity condition. However from the vanishing condition 
detqLj{X) = at points X = u, u* of the quantum determinant 

detgL = ^ ((LnL22 + ^22^11) - (L21L12 + 1^12^21)) (6.29) 

with L = L{X + m), we represent the Lax operator at these points both as a direct 
and an inverse projector. Therefore at A = z/, L°'^(z/) = ra{j)sj3{j) = Vp{j)ua{j) 
with r,s,u,v being two-dimensional vectors and similarly for A = u*. Therefore, the 
transfer matrix at these points becomes 

r{u) = tr {K{j - 2, J + 2 I u)L{j + 1 | iy)L{j \ u)L{j - 1 \ u)) 
= K^pU - 2, j + 2 I z/)r^(j + l)s,U + l)r^(j>5(j>5(j - - 1) (6.30) 

where K{j - 2, j + 2 | z/) = L(j - 2 | z/) . . . L(l | v)L{N \v)...L{]+2 \ v). We face 
with the problem, that in the quantum case due to the noncommutativity of rj,Sj 
at the same site, the last operator in (|6.30|) can not be dragged to the first place. 
However, using the inverse operator representation rs{i — l)sa{i — l) = Va{j — i)us{j — 
1) at the same site, we overcome this difficulty to get r(z/) = Qj^i {iy)M{j + — 1) 
where 
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M{j + l,j,j-l) = s{j + 1) ■ r{j) s{j) ■ u{j - 1) (6.31) 



with s{j + 1) ■ r(j) = si{j + l)ri(j) + S2{j + l)r2{j) etc. Similarly, we obtain r'(z/) 



J2f=i Qi-lii^) {s{j + i)L'{j I i')u{j — 1)) to yield finally the conserved quantity 



C, = c^lnr(A)|,=, = cr-i(z/)rV) 

TV 

= c + I, J, J - 1) + l)L'{j I z/)n(j - 1)) (6.32) 

i=i 

involving interactions of neighboring 3 sites, where M{j + — 1) is defined in 

( |6.31| ). Likewise one may derive also the conserved quantity from the expansion at 
z/*. 

Compared to the generalized spin method applicable to the symmetric models, 
the scope of the projector method seems to be more wide. 

6.2.3 Lattice nonlinear Schrodinger model 



The LNLS is a bosonic model with nonfundamental Lax operator ( [4. 361) , the Hamil- 
tonian of which can be constructed with both the above methods. 

To apply the generalized spin method to the lattice NLS, specialize to the spin 
s isotropic model and consider the Holstein-Primakov transformation ( |4.35| ) with 
spin parameter s = The Hamiltonian derived from RahW by fusing XXX 

spin-| models takes the form = /s(SjSj+i), where fs is a definite polynomial 
lU(j|| . Therefore from the equivalence of the spin model Lax operator with that of the 



lattice NLS: L'"''*(A) = j^a^L^^^ '^^^ obtain the corresponding Hamiltonian 

involving only the nearest neighbor interactions, which at the continuum limit A ^ 



reproduces the well known NLS field model ||108 . 

For constructing the lattice NLS Hamiltonian through the projector method, the 
Lax operator Ln{X) is taken in the form ( |2.39D with s = u = —-^ — -^-^ with the 



inclusion of a parity factor (—1)". Due to detgL{X) = |A^(A — z^i)(A — 1^2) = 0, with 
z/i = iuK, V2 = —i2n — vi the Lax operator L„(z/i) at A = z/i reduces to the direct 
projector at even n sites, while to the inverse projector at the odd n. Similarly, 
the reduction of -L„(z/2) is in the reverse order. Thus the necessary requirement of 
representing the Lax operator at the same lattice point as both direct and inverse 
projectors for the same spectral parameter value fails here. However, this difficulty 
is overcome by defining a new Lax operator as product of two Lax operators at even 
and odd sites: C{j \ A) = L{2j \ \)L{2j — 1 | A). At new points u = — ^ + in and u* 
the necessary condition is fulfilled and the Hamiltonian for the lattice NLS model is 
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obtained finally by adding the conserved quantities Ci coming from the expansions 
at A = z/ and X = u* along with a normalization factor involving A and an additional 
term proportional to the number operator N = Z^jLi'^j^iH' This discrete 
Hamiltonian H{A) is localized with interaction ranging over five neighboring lattice 
sites and has a complicated form, though at the continuum limit it reduces to the 
known field Hamiltonian ( |2.37|) . 



6.2.4 A simple lattice NLS type model 

To demonstrate that the projector method is applicable also to models with lesser 
symmetry, we look into the Lax operator 

MO^f^^'ff, (6.33) 

with N{k) = 1 + KA^(t){k)ip{k), listed in (|O0|) and constructed as the AKNS type 
integrable model ( |5.60| ) from the rational ancestor model in sect. 5. At the degenerate 
point ^ = 1^2 = when detgL = 0, the Lax operator ( |6.33|) becomes remarkably simple 



and turns into a field independent projector: L(0) = P"*". As a consequence, the 
construction becomes much easier resulting a local and simple Hamiltonian H along 
with other conserved quantities like the number operator N and the momentum P 



H = C3 = ^5](0(A; + l)^(A;-l)-(A(A;) + A(A; + l))0(A;+l)V^(A;) + (3fi:A2)"iA(A;)3) 
^ k 

P = C2 = Y.i<Pik + mk)-^N{kr), ^ = Ci = ^E^W- (6-34) 

It may be noted that the above conserved quantities of the discrete AKNS type model 
are not symmetric in and ip. The Symmetry however is restored at the continuum 
limit allowing the reduction (p = ip"^ and recovering thus the known conserved quan- 
tities of the NLS field model. Note that the expressions for the C„ (|6.34| ) including 



the discrete Hamiltonian C3 are much similar to the corresponding NLS field model. 



6.2.5 Lattice sine-Gordon model 

For applying the generalized spin method to LSG, we have to construct first the 
anisotropic spin s Hamiltonian by fusing XX Z chains. Subsequently, through bosonic 
realization of the quantum Uq{sl{2)) ( [4.9|) one can, in principle, construct the Hamil- 
tonian for the LSG model with nearest neighbor interaction [[108|] . However the 
construction in its explicit form is difficult to achieve. 
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The projector method can be apphed also to LSG model taking the discrete Lax 



operator as ( [4.10| ). However, the problem of representing the Lax operators both as 
the direct and the inverse projectors at the same degenerate point in the spectral 
parameter is difficult to solve and as a result the corresponding Hamiltonian H = 
^{Ci + C_i) and the momentum P = — C_i) obtained in this case turn out to 



be only quasilocal functions [50, 51 



6.2.6 Relativistic and nonrelativistic Toda chains 

In some models the construction of conserved quantities become much simpler due to 
specific structure of their Lax operators. Both relativistic and nonrelativistic Toda 
chains with the Lax operators constructed in sect. 5 are such examples. Though these 
are nonfundamental models, due to the triangular form of their Lax operators (see 
( f^.35| , |2.4l| )) the local conserved quantities can be obtained directly as the expansion 
coefficients of the transfer matrix 

N N 

r(0 = trTiO = trill ^n) = E ^nC + C-„r", 

n 

by multiplying number of Lax operators given in the explicit form (|2.35|) and 
then taking the trace. Thus we get for the relativistic Toda chain (RTC) conserved 
quantities like 

/+ = E (e'''''' (l + V^e^''-'~''^)) , I~ = J2 (e"'''''' (l + v'^e^'''"'^'^)) (6-35) 

i i 

where we have denoted = — C„_2(C'n)~^ and = — C_(„_2)(C'-n)~^- The Hamil- 
tonian and the momentum of this relativistic or discrete time model is constructed 
as if = |(/^ + and P = |(/^ ~ )• In the analogical way the conserved currents 
of all other models from the family of RTC found in sect. 5 can also be derived. 

Since rj —>■ corresponds to the nonrelativistic limit, we obtain the corresponding 
conserved quantities from the limiting values of in ( |6.35| ), resulting clearly the 
following expressions for the momentum and Hamiltonian of the standard Toda chain 



ra (ICTD 



P = ^^^ = Ep^ (6-36) 
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7 Exact solutions of eigenvalue problem, Bethe 
ansatz 



One of the most significant achievements of quantum integrable systems is that, these 
are some of the few interacting models, for which the eigenvalue problem in principle 
can be solved exactly. Though the energy spectrum of the Hamiltonian is the most 
important physical object, the integrable models allow to solve the eigenvalue problem 
for all its conserved quantities. The most efficient method for achieving this is the 
algebraic Bethe ansatz (ABA), though its coordinate formulation (called coordinate 
Bethe ansatz (CBA)) is usually more effective for finding the energy spectrum of 
concrete models. We describe both these useful methods and compare them to clarify 
their interrelations. 



7.1 Algebraic Bethe ansatz 

ABA is an universal and to a great extent a model-independent approach for solving 
the eigenvalue problem of quantum integrable systems. This method exploits signifi- 
cantly the integrable structures like the QYBE and is designed to find the eigenvalues 
of all conserved quantities including the Hamiltonian. 



We describe basic features of this method [0 in a general way to emphasize on its 
universal applicability, though for clarity we orient mostly to the simplest examples 
of XX Z and XXX spin-| chains. As we know, the monodromy matrix in ultralocal 
integrable systems : T{\)a = Hi' Lai{\) satisfies a representative equation given by 
the global QYBE ( p.ll[ ), which in component form reads 

E^}a (Tl^T^ = E {t;!t;i) r^^^. (7.1) 
{j} {j} 

Note that though T(A) given through explicit Lax operators Lj(A) carries detailed 
information about the concrete model, at the global QYBE level the information 
needed is only about the associated R matrix. Consequently, the R matrix (|4.6| ) with 
a multiplicative factor &(A)^^ for convenience: 

oil _ o22 ^ X « _ sin(^ + V) Di2 _ p2i _ T d12 _ d21 ^ _ c _ smr/ 
^11 - ^22 - / - ^ - ^.^^ , ^12 - ^21 - i> ^21 - ^12 - /i - ^ - ^.^^ 

(7.2) 

describes the trigonometric class of models, while its g — > 1 limit 

/ = ^. = ! (7.3) 
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answers for the models with rational R ( see sect. 2). We have witnessed in sect. 5, 
that the models may have common algebraic origin and the same ancestor models. We 
notice further in the ABA scheme, that the difference between the descendant models 
of the same class almost disappears, showing an elegant universality of the scheme. 
Therefore for the choice ( [7.2| ) the same ABA formulation with minor modifications 
can describe the whole trigonometric class of models , e.g. XXZ spin chain, sine- 
Gordon, DNLS model, Liouville model, relativistic Toda chain etc. Similarly, for 
( |7.3| ) almost the same scheme is applicable for the rational class, e.g. XXX chain, 
NLS model, Toda chain etc. 

ABA scheme aims to solve some general eigenvalue problem 

r(A) I >m= AmW I >m (7.4) 

for the spectral parameter dependent transfer matrix r(A) = traTa{X) = J2a Ta{X) = 
A{X) + D{X), where we have introduced the short hand notation 

T{\)n = A{X), T{\)22 = Di\), T(A)i2 = 5(A), r(A)2i = C(A). 

The operators T{X)ij act on the whole Hilbert space n = ®fhj, which is a product 
of N local Hilbert spaces hj at site j. Recalling that the expansion coefficients of 
Inr(A) = J2n CnX"', generate the conserved quantities C„, n = 1, 2 ... , one can solve 
the eigenvalue problem for all C„ simultaneously simply by expanding A (A) as 

Ci I >m= A:„(0)A-1(0) I >„, C2 I >m= {AL{0)A-\0)y I >m (7.5) 

etc. with commuting set of operators Cn sharing the same eigenvectors. The model 
dependence comes when a particular Cj is chosen as the appropriate Hamiltonian of 
the model. 

The basic trick in ABA is in making an ansatz for the m-particle state | >m as 

I A1A2 . . . A„ >= B{\i)B{\2) . . . BiXm) I >, (7.6) 

where | >= | >», lying in the whole Hilbert space is a pseudovacuum rep- 
resenting a reference state. For the XXZ and XXX spin-| chains for example, 

this is a state with all spins up defined by | >i= ^ g ^ ' which coincides with 

the ferromagnetic ground state of the models. To see the justification for the ansatz 
(|7.6|) , we have to convince ourselves that B{\i) acts like a creation operator on the 
pseudovacuum | >, while the operator C(Aj) is its conjugate representing an an- 
nihilation operator. For this we look into the example of XXZ chain, though one 
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may consider other examples as well. Starting from its local Lax operator Lj(X) 
and acting by the local creation and annihilation operators on local pseudovacuum 
as a~ I >j=| 1 >i, cr^ I >j= along with \ >j= | | >i, one obtains the 
upper-triangular form 

i^.(A)|0>,= I h (7-7) 




where a,/3,7 depend on the models. Taking now the product of such local states 
we obtain the action of the monodromy operator on the global pseudovacuum T(A) | 
>= nf ^i(A) I >i, which would again be in the upper-triangular form. We 
conclude therefore, that C(A) | >= and -B(A) | >=| A >, where | A > = 
fi'j(A)((®fc<i I >k)^ I 1 >i (®/>i I >;), mthgi = 7a^-'/5* is the superposition 
of single particle excitations at different lattice points and represents a 1-particle state 
with rapidity A. For diagonal entries we obtain in general 

A(A) I >= a^(A) I >, D(A) I >= /5^(A) I > . (7.8) 

The model dependence comes at this stage through a^, since they are determined 
by the structure of local Lj{X) operators as in ( [7.7| ). 

We observe that, it is possible to extract now the eigenvalues of r(A) = A + D, ii 
it can be dragged through the strings of B's in (|7.6| ), without spoiling their structure. 
This can be done by using commutation relations like AB = fBA, DB = f'BD and 
hitting the vacuum | > by A and D. The property (|7.8| ) then can rebuild clearly 
the same eigenvector \ (p >m and solves the eigenvalue problem. With this strategy 
in mind we write explicitly the QYBE ( [7.1| ) in component form 

AB = jBA-hBA, (7.9) 
DB = fBD-f,BD, (7.10) 

and [B,B] = 0, where A = A{X), B = B{X), D = D{X) and those with replaced 
arguments are denoted by / = f{fi — A), /i = /i(/i — X),A = A{fi), B = B{fi), 
etc. Here the functions / = f , /i = | and similarly are contributions from the 
i?-matrix ( [7.3|) or ( [7.2|) , while the operators A, B and A, B etc. are the components 
of T(A) and T(/i) appearing in the QYBE. We notice that ( [7.9[ - [7.10| ) are indeed 



the required type of commutation relations, but for the second terms in both the 
rhs, where the argument of B has changed: B ^ B, spoiling the structure of the 
eigenvector. However, if we put the sum of all such unwanted terms= 0, we should 
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be able to achieve our goal. Therefore ignoring the second terms for the time being 
and making use of the first terms only, as argued above we obtain 

m 

A(A) I A1A2 . . . A,, > = X{f{\,~\)A{\)\Q> (7.11) 

m 

D{\)\\^\2...X^> = X{f{\-\,)D{\)\Q> . (7.12) 

i=i 

Considering further the action ( |7.8| ) on the pseudovacuum, we finally solve the eigen- 
value problem to get 

r(A) I A1A2 . . . A„ >= A(A, {A,}r) | A1A2 . . . A„ > (7.13) 

where 

m m 

A(A, {Ajr) = n /(A. - A)«''(A) + n /(A - A,)/9^(A). (7.14) 
Note that equations ( [7. 1 l[ - [frT2D are universal within a class of models specified by 



the form of /(A — /x) and hence determined by the choice of the i?-matrix. In the 
explicit eigenvalue expression ( [7.14] ) the contribution from ,(3^ is the only model- 
dependent part and in case of XXZ spin chain for example, as seen from (|7.71 ), is 
given by 

«^ = sin^(A + |), /3^ = sin^(A-|). (7.15) 

Thus we have solved the eigenvalue problem for the whole set of conserved currents. 

However we have left open on the way two questions. Firstly, what constraint 
one gets when the unwanted terms are equated to zero, as we have done above. 
Secondly, what are the equations for determining the unknown parameters {Aj}. 
Curiously, the second question is answered by the first one itself , i.e. the equations 
for Aj's can be obtained by collecting meticulously all the unwanted terms following 
the prescription of and equating them to zero. It is however encouraging that 
there is an alternative and easy way to derive the same equation. Clearly, the general 
eigenvalue A(A, {Aj}) must be analytic in A, since it should generate the eigenvalues 
of all Cn through the expansion ([7.5|). Consequently, its residue at all singular points 
A = Aj,j = l,...,m, should vanish, since 6(0) = for both the rational as well 
as trigonometric i?-matrix. Therefore using the residue rule res\=\^ Hj /(Aj ~ A) = 
Ylj^k f i^j ~ ■^k) along with the symmetry &(— A) = —b{X), one gets from ( |7. 14| ) 
the equation 

=n^^r^= n ";f^"';' . *=i.2,...,™ (t.w) 
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The set of relations ( [7. 161) called the Bethe equation, are the determining equations 



for the parameters Xj,j = 1, . . . ,m and represents at the same time the condition 
for vanishing of the unwanted terms. Thus we get the basic equations of ABA as 
(|7.13[ - [7rT4D giving the eigenvalues through parameters Xj and the equations (|7. 16|) for 



determining such parameters. 
I. Lattice models 

The eigenvalue problem of various lattice models can be solved readily from the 
equations ( [7.13[ - [n^ ) and ( [7.16|) by knowing the functions /(A — /i) and calculating 
only a{\),P{X) for concrete models. Note that the function / given by ( [7.2| ) or ( |7.3| ) 
comes from the corresponding i?-matrix of the model and therefore is the same for 
all models belonging to one integrable class. Functions a, (3 on the other hand are 
related to the Lax operator of individual models and usually should be calculated 
separately in each case. Let us demonstrate this point again on the examples of spin 
models. 

7.1.1 XX Z spin-i chain 

The anisotropic XXZ spin-^ chain belongs to the trigonometric class with Hamil- 
tonian and the Lax operators given by ( |2.29| ). We have seen in sect. 6 how the 



Hamiltonian of the model can be constructed from the transfer matrix at an opera- 
tor level as if = Ci = t'{X)t^^{X)x=o. Here we aim to find its energy spectrum. For 
this we use the i?- matrix input ( |7.2|) along with the vacuum eigenvalue ( [7.15 ) in the 
Bethe ansatz equations ( [7.13| J7.14| , P'-1(^ ) to obtain the eigenvalue 



/v., ^ -^^ sin(A,- - A + 2) sin(A - A,- + 3?) ^ 

by making a shift A ^ A + ^ and the Bethe equation 

/ sin(A, + f) \^_ - sin(A,-A, + r/) 

l,sin(A,-f); sin(A, - A, - r^) ^ " > 

for j = 1,2, . . . ,m. Remembering that the Hamiltonian for the XXZ chain is Hxxz = 
Ci, we get the energy spectrum H^xz \ >m= E^} \ >m as 

m -I 

42 =A(A)'A-'(A) U.o=sm„i: ^,„(,^_,),,„(,^^,) +iVcot,, (7,19) 

Depending on the values of the anisotropy parameter A = cosrj, the model ex- 
hibits qualitatively different behavior. For A > 1 the ground state of the chain is 
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ferromagnetic, while A < 1 corresponds to the antiferromagnetic ground state with 
total spin = 0, attained at half filling m = y. A = ±1 is related to the isotropic case 
discussed below with +1 describing the ferromagnetic, while —1 the antiferromagnetic 
chain. 



7.1.2 XXX spin- 1 chain 

Isotropic XXX Heisenberg spin-| chain is given by the Hamiltonian and the Lax 
operator as in ( p. 381) and belongs to the integrable class with the rational i?-matrix 



(|7.3|) . However, since g — >■ 1 limit of the trigonometric model corresponds to the 
rational case, the eigenvalue problem ( |7.4| ) for the XXX model can be solved cheaply 
by taking this limit from the anisotropic result ( |7.17| - fr!T^ ). This in effect amounts 
simply to the formal replacement of sin A — A, cos A 1, resulting 

A...(A) = (A + n 7 : I + n , (7.20) 

for the eigenvalue expression along with the Bethe equation 



for j = 1,2, ... ,m. The energy spectrum for the m-spin excited state can also be 
similarly obtained from ( [7.191 ) as 

j=i 'I 



Note that for the XXX spin chain model (|2.38|) the parameter is to be set as 77 = ± 



■2 



in all the above formulas, which correspond to the isotropic case A = ±1 and is linked 
to the ferromagnetic or the antiferromagnetic ground states, respectively. Physical 
significance of the ABA results will become clear after we solve the same problem 
through the CBA and show their equivalence. 

7.1.3 Spin chains with impurity 

For describing integrability of the impurity models considered in sect. 6, the Lax 
operator L{X) at site s in the monodromy matrix may be replaced by L(A — u), 
keeping all other operators the same. As a result, for the eigenvalue solution of spin 
chains with impurity, the model dependent part in the ABA relations ( |7.14|) and 
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i\nM is changed as a^(A) -> a^"i(A)a(A - u) and P^{\) P^~\X)p{\ - u). For 



XXZ spin-| chain with impurities , for example, the expressions sin^(A±|) in ( |7.15| ) 
should be replaced by sin^~^(A ± |) sin(A — ± 2) influencing the result ( [7.14[ - [7?TE| ) 
accordingly. Therefore we may conclude, that the impurity parameter v acts like an 
extra (but known) rapidity parameter A^+i in addition to the set of such unknown 
parameters. The corresponding energy of the impurity model is 

4:i = 41 + cot(|-z.)-cot| (7.23) 

-,(m) 



where E)^l^ is given by the homogeneous model ( [7.19| ) without impurity. 



7.1.4 Inhomogeneous models 

In inhomogeneous models, extending the concept of impurity to every site the argu- 
ments of all Lax operators are shifted by Vj as (|6.14|) . Therefore in the ABA equations 
one should replace now 

TV N 

a^(A) ^\{a{\- z/,), /5^(A) v^). (7.24) 

The energy spectrum for the inhomogeneous XXZ model therefore results to 

N 

Ei:L = Etl + n cot(| - u,) - iVcot(|) (7.25) 

7.1.5 Other examples of regular models 

Let us start with the more general regular L{X) operator of the form (|6.5|) and extend 
eigenvalue relation ( [7.14|) for the spin models to 



and the Bethe equation (|7.16| ) to 



i) Wu-McCoy model 

This is a twisted XXZ model constructed in (|6.18| ), where the parameters are 
given by ( |6.19|) as &±(A) = e^^^sinA, a±(A) = sin(A + rf). Therefore the Bethe 



95 



equation ( [7.18|) of the original XXZ chain ,as seen from (|7.27|) is changed by the 
appearance of an extra factor e*^^ in its Ihs, while the energy spectrum remains the 
same as ( [7.19[ ) due to cancelation of 6'-factors. 

ii) XXZ in a magnetic field 

We have seen in sect. 6 that the color twisting of the XXZ chain amounts to 
putting it in a constant magnetic field. Therefore, using the related parameters (|6.2CI| ) 



one finds that due to = and ""ajlx)^ = ^ip^' Bethe equation remains 

the same, while an extra constant term ^cN appears in the energy spectrum coming 
from a^(A), which answers for the effect of the constant magnetic field |. 

iii) Free fermion models 

It is interesting to check how the freedom of the free fermionic models are reflected 
in the Bethe ansatz equations. From the parametrization (|6.21| ) of such a simple 



model we find that ^^liX) ~ ^"cL ~ ^' 'which readily gives the rhs of the Bethe 



equation (|7.27|) = —1, signifying trivialization of the m-particle scattering matrix and 
confirming thus the noninteracting nature of the model. 

Similarly for the free fermionic model generated from the nonstandard i?-matrix 
with ( |6.26D , we have = = 1 yielding again the trivial scattering matrix. 

7.1.6 String solution and thermodynamic limit 

i) String Solution 

Though the general solution of the Bethe equations (|7.18| , [7.21| ) is difficult to find. 



a particular form of complex solution, called the string solution of length m = 2M + 1 
can be given as 

Xn = x + in, n = -M, -M + 1, . . . , M - 1, M. (7.28) 

The corresponding excitation energy can be calculated by inserting the string solution 
in the expressions like (|7.19| , |7.22|) . For example, the XXX spin model gives 



ii) Thermodynamic limit 

The thermodynamic limit or the large system limit at — > cxd with finite lattice 
spacing A is important from the physical point, when the sum in the logarithm of 
the Bethe equations turn into an integral: jj: Z^jLi fi Ixo dxf{x) and the equations 
become easier to handle in many cases. For example, in antiferromagnetic XXX 



96 



chain with total spin zero for the ground state, we choose the half filling m = y and 
for calculating the momentum pj take the logarithm of the Bethe equations ( [7.21| ) to 
get 

N_ 

tan-i2A, = ^ + j^p^t^n-\X, - A,), (7.30) 

where the integers (or halfintegers) Qk take all successive values between —^^^ < 
Qk < -^^^ for the ground state. At the thermodynamic limit ^ — >• x, A^ — > A(x), 
and since for the ground state ^^^^ — "^^ S^t the ranges —\<x<^, — oo< 
A < cxD, which reduce the Bethe equation (|7.30|) to an integral equation 

1 / 4 1 

tan^ A(x) = TTx + dy tan^ {^{x) — ^{v))- (7-31) 

Differentiating ( [7.3 1| ) and introducing the density p(A) = (^^^^ \x=x{\)) of A in the 
interval d\ one gets 

The equation can be solved using the Fourier transform 

1 /-oo 1 

PW = 7r ^^T— ^e-^^= (7.33) 
2tt J-oo 1 + e 1^1 2cosh7rA 

Total energy of the antiferromagnetic ground state therefore is 

Eo-^ = N [ h{X)p{X)dX = -N r dX{—^)p{X) = -2iVln2. (7.34) 

Z J J—oo A H~ ^ 

The energy Eq of an anisotropic chain in the antiferromagnetic regime — 1 < A < 1 



is calculated in a similar way |53 



7.1.7 functional Bethe ansatz 

For some integrable models though the R matrix and L operators are known, due the 
specific structures of the Lax operator or the boundary conditions, sometimes it is 
difficult to construct the reference state | > forsuchmodels. As a consequence the 
standard ABA method is not applicable to them. Nevertheless, by the application of 
another method called functional Bethe ansatz (FBA) it is often possible 

to achieve factorization to reduce the many body problem to a single particle one. 
We briefly touch upon this method on the examples of Toda chains. 
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i) Nonrelativistic and Relativistic Toda chains 

We have seen the systematic construction of the i?-matrix and the Lax operators 
of the relativistic (or time discrete) and nonrelativistic (time continuous) Toda chains 
as well as the explicit derivation of the respective Hamiltonians. However due to the 
presence of exponential functions of fields in the Lax operators, the reference state 
with C(/i) I >= is difficult to find for them. We demonstrate in brief the 
application of FBA to the Toda chains following the argument of |59| . 



Using QYBE ( |7.1| ) we find first [C(^),C(C)] = 0, a crucial relation for the ap- 
plication of FBA, along with functional relations between C{E,)iA{C) analogous to 
(0): 

A(A)C(/i) = /(A - ^i)C{^^)A{\) - /i(A - M)C(A)A(/i), (7.35) 

and similar relation between C{^),D{C). Defining now commuting operators 9j,j = 
1, . . . , N — 1 as C{\ = 9j) = and conjugate operators A~ = A{\ — > 9j), A^ = 
D{\ — i> 9j) with the proper operator ordering prescription ||5^, we get from (|7.35| ) 
and from the relation between C, D the simple equations 

A7C(/i) = f{9, - /i)C(/x)AT, A+C(/x) = /(/i - e,)C(/i)A+ (7.36) 

Let us consider the case of the nonrelativistic Toda chain related to the rational R- 
matrix, for which we have f{9j — jj) = and interpolating the zeros at 9j we 

construct C(yu) = ]Vj=i{l^ — 9j), using the input of the trivial commutation of C(/i) 
at different /i. From ( |7.36| ) we conclude that the operators A^ shift the operators 
9j by =F?7 and since it is true for the symmetric polynomials in 9j, one gets for the 
symmetric function 0(6'i, . . . , 9n) the important relation 



t{9,)<P{9^, . . . , 9m-i) = {k-2K^ + K~-^K^)<P{9^, . . . , 9^_^) 
= Kh-''<f){9i, 77, ... , 9n^i) + K-h''<f){9^, ...,% + r/, ... , 9n-i){7.37) 

where i'^'^ is the eigenvalues of AJ. To achieve separation of variables for the — 1- 
particle eigenfunction we represent 0(6*1, . . . , 6'Ar_i) = Y\f~^ to get 

t{9j)(j)j{9j) = i-^K^^(f)j{9j -r])+ z^K-50j(% + T]) (7.38) 

The FBA method described above is applicable also to its relativistic (discrete 
time) generalization related to the trigonometric i?-matrix by expressing the 
spectral parameter as ^ = e^^ and considering the zeros of C (/z) at 9n with , = e*^" . 
The separation of variable equation takes in this case the form 

^(ei)0i(e,) = t-^'f^Hji^je-'^) + ^^«:-^0i(e,e^'') (7.39) 
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Possibility of such solutions is also related to the pair propagation through a 
vertex for the i?-matrix [^21, i.e. the existence of vectors /j^^C""*) such that 



R{^Yj^v^pK-^ = v^^h°^^. A recent application of FBA to solve XXZ spin chain with 



antiperiodic or more general boundary conditions can be found in |^8|, |109|| . 
Field models 

We focus here on the specializations of the ABA, that one has to make for appli- 
cation to the field theoretic models, over that for the lattice models described above. 
We have seen in sect. 4 that the integrable field models like NLS, SG , DNLS etc. 
can be lattice regularized with exact Lax operators exploiting their underlying alge- 
braic structures. Therefore such lattice versions of the field models defined in a finite 
interval [—/,/], I = NA with lattice spacing A satisfy the global QYBE ( [7.1| ) and 
can be solved through the same ABA scheme described above. 

Our aim however is to explore the continuum models with A — > and at infinite 
interval / ^ oo. Involvement of these two limiting procedures makes the related 
ABA a bit different from and slightly difficult than the ABA for the discrete models. 
However, we achieve in turn some simplifications, e.g. we completely get rid of the 
problem with unwanted terms and consequently of the Bethe equations. We may 
consider first the A ^ limit by defining 

2N 

TKA) = hm n ^n(A) (7.40) 

n=l 

which still satisfies relation (|7. 1|) . Another advantage is that, since we are interested 
in field models at A 0, in most of the cases it is not necessary to deal with much 
complicated forms of the exact lattice versions of the Lax operators (see the list in 
sect. 2), but is sufficient to use L„ correct up to order A, i.e. Ln = 1 + Af/„(A), 
where Un{)^) is the naive discretization of the Lax operator U{x, A) of the continuum 
model. 

Let us find now the asymptotic at n oo, when the field operator, say ipn 
vanishes in the weak sense | ipn 0, resulting 

L„(A)^V(A). and L„(A) ® L„(/i) ^ H^(A, /x), (7.41) 

which may differ from V^(A) V{fi) due to normal ordering of the operators. At the 
infinite interval limit the QYBE (|7.1|) is needed to be regularized, for which we define 
the scattering matrix as 

T(A) = lim \/(A)-^T,(A)y-^(A) = | 'J/^^ 7;;;^ I , e = T. (7.42) 
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Observe that at oo interval limit SU {2){SU{1, 1)) symmetry is usually restored for the 
monodromy matrix. The regularized QYBE using the equation R{X — fi)W{X,fi) = 
W{fi, X)R{X — fi) now looks like 



R+{X - /i)T(A) ® T{fi) = (/ ® T(/i))(/ ® T(A))i?_(A - fi) 



(7.43) 



with R± = S^^RS±, where S± = S±{X, fi) = lim^^oo W^^{X, fi)V^^{X) O V^^{fi) 
and S± = S±{n,X) carry information about the asymptotic behavior. Unlike the 
lattice models involving the R matrix only, the explicit form of S± and consequently 
R± matrices become model dependent. We consider below some examples which will 
clarify this point. R±-ma.tnx generally takes the form 



( f 



R^ 



1 fi± 
/it 1 



V 



f ) 



c 
a 



(7.44) 



with /i±(A,/i) = lim;_^oo /i(A, /u)e^*^'''^'^^5, where the function p depends on the con- 
crete model. The limits of such singular terms must be taken carefully in the principal 
value sense 



lim V 

Z— too p 



±.m5{p) 



and after properly normal ordering the operators. 



7.1.8 NLS model 



Let us elaborate the subsequent steps on the example of NLS field gl], |T6[ model 
involving the fields ipiip'^ with [%p{x),'ip'^{yy\ = 6{x — y), where 



'^n = -r ip{x)dx 



(7.45) 



and described by the Hamiltonian and the Lax operator as in ( |2.37| ). The discrete 
Lax operator for our purpose may be chosen simply by discretizing ( |2.37| ) as 



L„(A) = I - iA/:„(A) = I - 



A r]2ip^ 
■r]2^jj 



-X 



(7.46) 



Therefore the asymptotic expressions at give V{X) = I—iAXa^, W{X.p) = I — 
iA{Xa^®I + pI^a^ — irja' ^a^) which in turn construct S± and R± matrices. Using 
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the associated rational /2-matrix (|7.3|) we get R± in the form ( [7.44|) with /i±(A,/i) = 
hm;_,oo A^^*^'''^"^''^ ~ 2Tciri6{X — /x), which yields from ( |7.43|) the relations 

AI3 = fBA - 2mr]6{X - ^i)BA (7.47) 

with [A, y4] = and the notations same as in ( [7.9[ - [7.10] ). In handling such singularities 
it is always better to extract the limits only at the final stage after performing all 
other operations. 

Compared to (|7i9| - |7.10|) in place of unwanted terms there appears in (|7.47|) a 
term with 5{\ — fi), which disappears for A 7^ /i. Therefore we can repeat here the 
arguments of the above ABA scheme for constructing the reference state | > as 
well as the m-particle eigenstates as | >m=| A1A2 ■ ■ ■ Xm >= -B(Ai)-B(A2) . . . -B(Am) | 
> . The eigenvalue expression consequently is given again by ( |7. 14| ) considering the 
regularization ( fT^ and yields limjv^ooa^ | >= lim^v^oo e*^^'"'^Tz(A)e*^^'"'^ =| 
> However, due to the absence of any constraint like Bethe equation, parameters 
{Xj} are arbitrary in continuum models. This is an important and major difference 
between the ABA for the field and the lattice models. 

Note that the NLS model is intimately connected with the spin models through 
Holstein-Primakov transformation ( ^4.35| ) and belongs to the same rational class as 
the XXX chain. We can formally derive the eigenvalue of A{X) for the NLS from 
( fr:20|) as AnisiX) = n"=i considering = 1 and a shift X ^ X - q. 

Consequently the eigenvalues for the conserved quantities for the NLS system are 
generated by the expansion InAnisiX) = J2n C-nX~^ at A — 00 and = Z^jLo C'nA" 
at A — *• in correspondence to the respective conserved opeartors. Since Hnis = C*_3, 
one finds therefore the energy spectrum for the NLS model as Enis = H^=i Aj- The 
case T] = — I ?7 |< corresponds to the NLS model with attractive interaction and 
allows also bound or soliton states given by the string like solutions Xj = Pq — i \ 
?7 I ( j — ^^Y^), i = 1, . . . ,m. Due to such choice of Xj the poles of the individual 
factors in Anis{X) are canceled by the zeros of the neighboring factors. As a result we 

get AgoiitonW = Tf^, which though constructed from m-number of particles as 

a bound state, mimics clearly the structure of a single particle and hence is called 
soliton. Note the curious fact that the same structure with zeros and corresponding 
poles at discrete points in the complex A plane also appear in the classical soliton 
solution for the scattering matrix element a(A), as mentioned in the introduction. 
Thus for such m-particle bound state, the particle number N = C_i, the momentum 
P = C-2 and the energy spectrum Em = C_3 can be obtained as 

N = m,P = Po, Em = ^-^im^-m) (7.48) 

m 12 
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The other set of conserved operators C^, n >0 leads usually to nonlocal quantities. 

Interestingly, the Hamiltonian of the NLS field model is equivalent to the 6- 
function Bose gas as the interacting N body system with the Hamiltonian 

N fj 

^ = -EnZ2+2^E^(^^-^.) (7-49) 
j=i ^-^j i<j 



Therefore the energy eigenvalue problem H \ n >= En \ n > for the NLS model ET 



also coincides formally with that for the (5-function Bose gas solved in |110, 56 



7.1.9 DNLS model 



Unlike the well known nonultralocal derivative NLS pBI, the DNLS model given by 



(p.27|) is ultralocal with the field commutation relation [ip{x), ip'^iy)] = — sin h6{x — y) 



and, as shown in sect. 5, has a remarkable connection with the quantum algebras and 
g-oscillators. We show here how one can solve the eigenvalue problem of this quantum 
field model exactly through the ABA scheme [^. Though the exact lattice version 
for the Lax operator of the model is available as ( |5.20| ), it is sufficient to consider 
such operator correct up to order A : L„(A) = / + 2A£„(A) by using (|7.45|) , where 
Cn is a simple disctretization of the field model Cdnis given in ( p.27|) . The coupling 



constants K-t are nontrivially renormalized in the quantum 

Since the model belongs to the trigonometric class, the ABA should incorporate 
the case ( [7.2| ) and follow the formalism described above for the field models involving 
the limits N ^ oo, A —* 0. The structure of the Lax operator and the commutation 
relations derive the asymptotic matrices (|7.41|) as 

ViX) = (l-Y^Va 



4 

The R± matrix has the structure of (|7.44|) with 



/i±(e,C) = lim /i(A,/x)e^^(«'-^')i = 2<^sinr/5(e^ - C') 



I— too 



and the QYBE relations have the form ( [7.43[ ) with specializations for as mentioned 
yielding 

AB = fBA - 27tX'^ sin r]6{X'^ - fi'^)B A. (7.51) 

Repeating the steps described above and using further the conditions like A{X) \ 
>=| >, one gets the eigenvalue solution AdnisiX) = Djli sinii(A-v-il) - '^^^^ allows 
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also the string-like solution ( [7 .281) as in the NLS model, resulting the m-particle bound 
state with the spectrum 



Co = rjm, , Ci = 2psm{rim), C2 = SrjHm = sin(2?7m) 



(7.52) 



and the binding energy 

AE = H„ 



mHi 



(sin 2am — m sin 2a) < 0. 



7.1.10 Sine-Gordon model 

This relativistic field model has been studied intensively highlighting its various as- 



pects p4|, We have generated its Lax operator in sect. 4 as a direct realization 
of the well known quantum algebra and discussed the problem of constructing its 
quantum Hamiltonian in sect. 6. Here we apply the ABA scheme to the model for 
determining the exact energy spectra. 

i) Lattice regularization 



For solving the SG field model by lattice regularization ||2J], one finds that the 
naive disctretization does not work here and one has to start from the exact Lattice 
Lax operator given in ( [2.311 ). Since the SG model is associated with the trigonometric 



i?-matrix ( [7.2|) same as the XXZ spin chain, the discretized model follows the same 
ABA scheme giving relations like (7^, 7.1C ), ( [7.17 ), ([7.18| ) etc. with specialization 

Di{X) = Ai{Xy and Ai \ >= e"(^)' | >, Di \ >= e"*(^)' | >. However there are 



also some major differences. Firstly, the action of the discrete Lax operator ( [2.31 



on the pseudovacuum: L„ | > does not reduce to the triangular form, which should 
be tackled by defining a new Lax operator 



L„(A) = L„(A)L„+i(A) 



(7.53) 



as the product of two neighboring operators. 

Unlike the above discussed models, the physical vacuum Q is obtained here by 
filling the state | > by particles as = llpB{Xp + «f ) | > [H]. Therefore 
1-particle state is described as the excitation over the physical vacuum as | 1 >= 
-B(Ao) Y[pB{Xp + i^) I > . At the limit / 00, the second term in ( [7.14[ ) becomes 
negligible and the corresponding eigenvalue of | 1 > normalized by that of the vacuum 
state gives 



A(A,Ao) 



A(A,Ao 
A(A) 



/(Ao-A) 



np^fe/(Ap + 2f -A) 
^p/(Ap + ^f -A) 



(7.54) 
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Now making a shift X X — ^ with a scahng of A, Aq , using the i?- matrix input (|7.2| ) 



and finally taking the length Z — > oo in (|7.54|) , we get for the m-particle excitation 



. / , N , TT- tanh(A — A^ — i-y) , 

A{X) m >= TT -4- m > 7.55 

^ ^ ' tanh(A - Afc + ^7) ' ^ ^ 

where 7 = and | m >= B{Xi) . . . B{Xra)IlpB{Xp + i^) \ >. Note that 

each factor in ( |7.55| ) corresponding to a single particle excitation has two zeros and 
corresponding poles in the complex A plane, which describe a breather like quantum 
state representing bosonic particle of the theory. One can consider also a M-particle 
bound state by considering a string-like solution for the parameters Xk- Curiously, 
the more fundamental kink state with nontrivial topological charge is difficult to 
obtain in this model. Expanding lny4(A) from ( [7.55|) as I]„C„A" at A ^ and as 
J2n C-nX~"' at A — » 00, we can get the energy and the momentum eigenspectra as 

H = ^XC-^ - C,), P = ^{C^i + C,). 
8 I 7 I 8 I 7' I 

ii) Continuum approach 

Alternatively, one can solve the eigenvalue problem of the SG field model starting 
from the continuum Lax operator ( [^.25D itself by considering carefully the short 



distance behaviors ||112|| . This more direct approach leads also to the same particle 
spectra given through the breather state ( [7.55| ), but draws an important conclusion 
that such solutions are possible only for the parameter value | 77 |< vr. For getting 
the asymptotic solution lim^^oo Ln ^(A) in the diagonal form, one has to make a 
gauge transformation over the Lax operator ( |2.25| ) amounting to a cyclic replacement 
of 0"", a = 1, 2, 3. As a result the i?-matrix changes to a form 




X = [^^7(A,/x), Y=(^l-^^tanh'-^jg{X,f,) (7.56) 

where g{X, ^) = f{X,fi) (l — coth . Therefore in place of (|7.2|) one gets 

the function / = | as /(A,/i) = ' which after defining 7 = 2tan^-'^?7 

and making a shift A —* X + i^ coincides with the factors appearing in ( |7.55| ) . At the 
continuum we construct ( [7.44| ) using this modified i?-matrix and the QYBE (|7.43| ) 

to get 

A{k)B{k') = f{k,k')B{k')A{k) + 2'K^ko6{ki- k[)B{k)A{k) (7.57) 
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with ki = 2sinhA, ko = 2 cosh A, which consequently derives the same breather-hke 
spectrum ( |7.55|) for the eigenvalue expression. Note that the only difference appearing 
here is in the value of the renormalized coupling parameter 7. 

The SG model in the light-cone coordinates using ABA has also been solved in 



7.1.11 Liouville model 



This relativistic field model is given as in (|2.26|) . Lattice regularized version is usually 



represented by the exact discrete Lax operator (|2.32|) . However in (|5.16| ) another 
such Lax operator is derived, which appears to be more suitable for the Bethe ansatz 
solution Like the SG model one has to use the product of Lax operators (|7.53| ) 



in constructing the pseudovacuum. Since the model belongs to the trigonometric 
class, the ABA scheme goes parallel to the spin-| XX Z chain described above with 
eigenvalue (|7. 14| ) , Bethe equation ( [7.16|) etc. by incorporating (|7.2|) . The difference 
appears only in the model-dependent part ( |7.15| ) defining the action of A{X) and 



D{\) on the pseudovacuum. This yields the Ihs of the Bethe equation as 



Comparing with the Bethe equation for the spin-S* case one may conclude, that the 
lattice Liouville model corresponds formally to a spin (— |) XX Z chain with an extra 
phase factor e*''^ EBI. Recall that, similar phase factor arises under twisting through 



(|7.27| ). The momentum and energy eigenvalues of the discrete Liouville model is the 
same (with opposite sign) to those of the spin +|, XX Z chain. Since discretized 
model is defined in a finite volume 27r = AA^, the N 00 limit would lead to the 
Liouville field model yielding the corresponding eigenvalues at this limit. Bound state 



solutions, conformal properties and other details of the model can be found in . 



7.2 Nested Bethe ansatz 

The ABA described above is designed for the SU (2) related models with 2x2 Lax 
operators. We mention in brief how to generalize this technique for models with 
SU{N), N > 2 having N x N Lax operators. Examples of such models are the Perk- 
Schultz ||9^ type model associated with the trigonometric i?-matrix (|4.31| ) and the 



vector NLS model with rational R. The method, called nested Bethe ansatz (NBA), 
consists of a recurrent application of the ABA or their nesting for subsequent solutions 
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of the higher rank problems using those of the lower ones. We demonstrate it for the 
case N = 3, since the generalization to N > 3 then becomes straightforward [p^ . 



We start as in the standard ABA with the QYBE ( |7.lD , where the R matrix 
is taken now as a more general solution ( obtained also in (|4.15|) through Yang- 
Baxterization) 



i,k 



ki 



(7.59) 



with /, /i given by the same expressions as in (|7.2| , |7.3| ) and the monodromy matrix 
is in the 3x3 form 





( Too 




B2 


\ 


T(A) = 




Tn 








I C2 


T21 


T22 


J 



This gives the relations 



and 



TaaBp 
fBaBp 



fBaToo — flBaT{ 



00) 



(1)75 



(7.60) 
(7.61) 
(7.62) 



with all Greek indices taking values 1 and 2 and with the notation Tqq = Too (A), Too = 
Too(/^) etc. Note that the relation ( [7.60|) is similar to ( [7.9|) , while ( [7.61| , [7.62D are 
complicated and unlike the usual ABA case do not preserve the same set of operators 
we started with. Nevertheless, there are similarities in the structure and the 2x2 
matrix R^^^ appearing in the above relation is again a i?-matrix solution with one 
step lower rank. 

Observing in (|7.62|) that a particular string of 5's generates its other combina- 
tions, we guess the Bethe ansatz for the eigenvector in the form of superposition of 
all such strings of -B's as 

I >=| A1A2 . . . A„ >= E Ba,{X,)B,,{X2) . . . 5,„(A„) I > F^'-^-, (7.63) 

{a} 

where the operator coefficients F"'^'""'" would be determined later. Using the similar 
pseudovacuum reference state (|7.8|) one obtains 

Too(A) I >= a^(A) | >, Tn(A) | >= T22{\) \ >= /?^(A) | > . (7.64) 

For this model we have a = /(A) and (3 = 1. We assume for the time being, that 
the unwanted terms in the commutation relations (|7.60| - fr6ll) (those with changed 
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arguments like Ba, Tap) are simply absent. We will find later that as in the ABA, the 
Bethe equations for this system represent also the condition necessary for vanishing 
of such unwanted terms. 

Repeating the standard arguments described above we obtain the eigenvalue of 
^M) as 

Too(A) I n > = Aoo(A I Ai, . . . , A„ I n >, 



. N 



a 



(A)n/(^^-^) l^> > (7-65) 



1=1 



However, the relation ( [7. 6 ID can not solve the eigenvalue problem directly and instead 
yields 

n 

Y^Taaix) I n >= (n B,,{Xi)) r(i)(A | {Mt;::^ f^'-"" I o > . (r.ee) 



(=1 



It is curious to note that the t^^^ operator appearing here with t^^\\ \ = 
E«T«(A I {A,}) = tr(T«(A | {A,})), where 

I {Az}) = i?«7^^HA - AO ■ ■■R^^]2±{X - A„)T.,„(A) 

= (L,,(A-Ai)---L„„(A-A„)T(A)U (7-67) 

represents the transfer matrix of A^ + n site chain with n inhomogeneous sites (coming 
from R^^^ matrices). One can continue the arguments further by attempting to solve 
the eigenvalue problem of t^^\\ \ {A;})^^.'.'.^" pa-i- a„ | g ^ since we see immediately 
that ( |7.67| ) involves only 2x2 monodromy matrices and hence is the standard problem 



we have already solved. This nesting is the basic idea in NBA and we choose therefore 

fpayan _ I ^^(l) ^a;---a„ 

1 I /V]^ , . . . , ^ 

= (5«(Af ^) ■ ■ ■ 5«(Ai\)) I O(^) (7.68) 

with B^^\xf) and | 0^^) > living on a lattice of n sites and thus is an element of 
direct product of n Hilbert spaces. The part of r'-^^ involving Ta-y„ however acts on 
nonoverlapping A^ product Hilbert space and hence commutes with B^^\\^py?, and 
can hit directly the vacuum | > . 

Using therefore the eigenvalue result (|7.13|) , with trivial modification due to in- 
homogeneity in ^^^(A — \j) and its form ( [7. 591 ) affect the action on the vacuum as 

n 

^^'H^ I {\}) 1 0^'^ >= n /(A - 1 0^'^ >' ^^'H^ I {\}) 1 0^'^ >=i o« >, 
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along with T^a | > as in (|7.64|) . Thus we get the part of eigenvalues as 



ni n ni 

(r«(A I {A,})F) I >= /3^(A) | > | J] /(A^ - A) IT /(^ - AO + 11 /(^f - A) 1 F. 



(7.69) 



The final expression for the eigenvalues of the original transfer matrix r(A) 
trT{X) = j:LoTaa{X) is 



r(A) I n >= ^ A„J n >= A(A | {AJ, {Ap^}) | n > 



a=0 



can be obtained in the closed form by adding ( [7.65| ) and ( |7.66|) incorporating (|7.69|) . 

We observe that in the NBA additional sets of parameters appear, for which 
the determining equations called Bethe equations are obtained from the vanishing 
of unwanted terms (already assumed). One can derive such equations by carefully 
keeping track of all such terms throughout the nested steps ||113|| or an easy way out 
as before is to obtain the same Bethe equations from the vanishing conditions of the 
residues of A(A | {Ai},{Ap'*}) at A = {Xi},j = l,...,n, and A = {Xf^},l = 
1, . . . ,ni. Thus from the condition resA{X)x=\^ = res(Aoo(A) + Aii(A))a=a™ = 0, 
(since A22(A) gives no contribution) we obtain the 1st set of Bethe equations 



n f(x 



f{Kn — A;) 



A.„ 



ni 

n /(Aj 



(1) 



Am) 



(7.70) 



where /(A^) 



a(Am) 



m 



n. The 2nd set of Bethe equations similarly can be 



obtained from the condition resA(A) = res(Aii(A) + A22(A)),^,(i) = as 



n /(Ai^^ - Az) = n 



1=1 



(1) 



fiK" - Af ^) 



j^k f{Xj - Afc 



(7.71) 



with k = 1, . . . ,ni. 

Choosing the function /(A) as ( [7.2|) or as (|7.3|) we can apply the above result for 
solving concrete quantum models with internal degrees of freedom and belonging to 
the trigonometric or the rational class, respectively. Examples of such rational models 



are the vector NLS model |114|, matrix'NLS model, quantum N— wave system |57, 115 



etc., while the Perk- Schultz model ||116|| or the generalized six- vertex model ||117|| , the 
anisotropic n- vertex model etc. are the examples of higher rank models belonging 
to the trigonometric class. 
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7.3 Coordinate Bethe ansatz 



The coordinate formulation of the Bethe ansatz (CBA) [Q, ^ |118|| aims to solve the 



Hamiltonian eigenvalue problem directly, without assuming explicitly any integrabil- 
ity condition. However for integrable models allowing ABA treatment, the CBA is 
also found to work successfully. Let us demonstrate it on the example of anisotropic 
XXZ spin-i chain we have solved already through ABA. 

7.3.1 XXZ spin-i chain 



The Hamiltonian of the model ( |6.8| ) after addition of an irrelevant constant term is 
written as 

^ _ _ A 

H^xz = - E f^n+l + (^n ^^n+l + -^((^l^^l+l " 1), (7-72) 
n ^ 

where values of the parameter A = cosr] defining the anisotropy of the model de- 
termines also its anti-ferromagnetic or the ferromagnetic behavior. We work now in 
the coordinate formulation, i.e. explicitly on the lattice, where the spin states can be 

either wp or dovjn with basic vectors e+ = ( | or e~ = | 1 = cr^e+, and 

we consider the spin excited states | xi, . . . , >= o^^ . . . a^^ | >, with m down 
spins e~ at lattice sites xi, . . . , in the background of the rest — m up spins. 
The set of such basis vectors are orthonormal satisfying 
< x']^, . . . , x^ I Xi, . . . , Xm >= 0, for different arrangements of {x'}, {x} and 
= 1 for {x'} = {x}. The m = or the vacuum state | > state is the ferromagnetic 
ground state with all spins up. 

The ansatz for the eigenvector is given as the superposition of all such states: 

>= E ^(xi, . . . ,x„) I Xi, . . . ,Xm > (7.73) 

X\<X'2.<...<Xm, 

the explicit form of which will be specified later. Equations for determining the 
functions ^'(xi, . . . ,Xm) can be obtained from the eigenvalue equation H \ \E'm >= 
^m. I > by the Hamiltonian ( 7.72 ) acting on the vector (|7. 73|) . We start from the 



elementary knowledge cr^e^ = 0, cr^ej = ej, cr^e^ = ±ej, to obtain 
(c^n t^n+i)e^ ® e^+i = , ((Tja^+i)e^ ® ej+i = ® e^+i. 



along with the trivial actions of all the above operators on ® ^n+i- interesting 
conclusion follows now, that on local excitations involving the term . . . e"*" ® e~ ® 
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e+ . . . =1 . . . , 0, x„, 0, . . . > the XY part of the Hamiltonian h^y = Yin ^l^y^"^-, 
where = -(a+o-^+i + o-^cr++i) acts as 

-/li"^'' I ...,0,x,_i^,... >=| ...,0,0^,... > (7.74) 

shifting it towards right by one site, provided that site is not aheady occupied. In 
the analogous way 

-h'i^^ I ...,0^,0,... >=| ...,x^0,0,... > (7.75) 

acts with a similar left shift. On the other hand, the remaining part of the Hamilto- 
nian hz = Yn h^'^^^, h^'^^^ = — y(c"^o"^+i — 1) preserves its position yielding 

-/i^-^'' I ...,0^,0,... > = -A I ...,0^,0,... >, (7.76) 
...,0,Xi^,... > = -A I ...,0,Xi^,... > . (7.77) 

For better understanding of the effect of spin interactions on the eigenvectors, let us 
look into the m = 2 excitations with | >= J2xi<x2 "^{^1,^2) \ xi,X2 >, before 
going to the general case, where + sign signifies the particular ordering of Xi < X2- 
Since the Hamiltonian involves only nearest neighbor interactions, we consider first 
the situation xi + 1 < X2, when xi, X2 are non-nearest neighbors, resulting the action 
of the Hamiltonian on eigenvector | \E'^ > as the noninteracting one. Therefore 
implementing the above mentioned properties of h^y, we obtain the action of the 
total Hamiltonian in the eigenvalue equation. Considering further the orthonormality 
of the set of basic vectors | xi, X2 > we get an equation for the functions \E'(xi, X2) as 

-(^(xi-1, X2)+^(xi+l, X2)+^(xi, X2-l)+^(xi, X2+l))+4A^(xi, X2) = ^2*(a;i, X2] 

(7.78) 

Note that since Xi,X2 are not nearest neighbors a particular state | Xi,X2 > can 
be reached under the action of h^y, from both right and left sides as in (|7.74| , |7.75| ). 
This results symmetric expressions of \E'(-, ■) in ( |7.78| ) and due to the same reason the 
contributions of hz in (|7.78|) are from both (|7.76| ) and ( [7.77|) summed to 2A\E'(xi, X2), 



coming from each of the excitations at Xi and X2- 

For incorporating interactions we consider now Xi,X2 as nearest neighboring 
points with X2 = Xi + 1. It is crucial to note that any state | Xi,X2 >=| x,x > 
can only be reached in this case by a right shift of a;i — 1 as in ( [7.74| ) and a left 



shift of X2 + 1 as in ( [7.75| ), resulting the presence of only such functions ^'(■, ■) in 
its eigenvalue equation. Furthermore, since only the left neighboring site of Xi is 
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empty and correspondingly the right site for X2, the contributions under hz for the 
excitation at xi come from ( |7.76|) only, while that for X2 is only due to ( |7.77| ). As a 



result compared to ( |7.78| ) less terms appear now in the eigenvalue equation giving 

- (^(a; - 1, X + 1) + a; + 2)) + 2A^(a;, a; + 1) = ^2^(3;, x + 1). (7.79) 

However we demand that our eigenvector | \1'^ > should be the same for all values 
of Xi < X2, i.e. (|7.78|) by setting formally X2 = Xi + 1 should coincide with (|7.79|) . 
Consequently the condition for vanishing of their differences gives the constraint 
equation 

^(x + 1, X + 1) + ^(x, x) = 2A^(x, X + 1). (7.80) 
The difference eqns. ( [7.78|) and ( [7.801 ) being at our disposal, an ansatz for the solution 



is made in the form of superposition of plane waves 

^(xi, X2) = Ai2e^(Pi^^+P2^2) + yi2ie*(P2^i+Pi^2)^ (7 81) 

where the coefficients ^12,^21 or rather their ratio and the momenta Pi,pi are 
unknown parameters. They will be determined as shown below, by the equa- 
tion (|7.80| ) and the boundary condition. Generally speaking one should denote 



A12 = Ap-i^p^: ^21 = ^p2Pi' stress the ordering of Xi,X2 and Pi,P2- Since the 
eigenvalue E2 must not depend on the positions of the excitations, we may find it 
from ( [7.7^ ) inserting ansatz (|7.81|) to yield 

2 

E2 = 2^(A-cospj). (7.82) 

For determining the coefficients we similarly obtain from ( [7.80| ) 

_ei0{pi,P2)^ (7.83) 



A12 g{P2,Pi) _ je(puP2) 



Ml 9{PuP2) 
where 

g{p^^ P2) = 2Ae*^'^ - e'^P'+P^^ - 1, (7.84) 

9{pi,p2) = —0{p2,Pi) may be interpreted as the phase change suffered in the scatter- 
ing of two pseudoparticles (magnons) due to the exchange of momenta. 

For deriving now the equations for the parameters pi,P2, we invoke the boundary 
conditions considering them to be periodic. That is we assume the chain to be closed 
by identifying the states at sites X and x+A^ given by | '^ti<x2 ^~ I '^xi+n=x[>x2 ^' 
where | > denotes the state with reverse ordering of x^s. However it is easy to 
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see that | > does not give any new solution and can be obtained from ^+ > by 
just interchanging xi ^ X2- Therefore the periodic boundary condition is compatible 
with the relation 

^{xuX2) = '^{x2,xi + N), (7.85) 
which with the substitution ( [7.81| ) yields the system of equations 

defining pi,P2- 

Thus we have completely solved the eigenvalue problem for the m = 2 case by 
finding explicitly the energy spectrum and the eigenvector. Repeating similar ar- 
guments one can solve the general case for arbitrary m < N, yielding the energy 
value 

m 

= 2E(A-cosp,) (7.87) 



and the eigenvector ( [7.73|) with 



vl/(xi, ...,xj = E^Pd....,-) e^^-^ ''^'^^'^ (7-88) 
p 

where -P(l, . . . ,m) stands for the permutation of all m indices 1, . . . ,m involving 
ml number of terms in the expression and P{pj) corresponds to the the same set of 
permuted indices j in pj. The amplitudes are connected through the phase factors 6 
as 

^l,...,m 

where the summation of the ^-factors are extended over those pairs of j, k which are 
necessary to interchange to arrive from the arrangement of indices 1', . . . , m' to the 
initial arrangement 1, . . . , m. More physically, for each interchange of particles at xj 
and Xk the scattering amplitude picks up a phase factor e'^^Pj'Pk) and the total change 
in amplitude is given by the multiphcation of such phase factors due to two-body 
scattering. This fascinating fact of the factorization of many-body scattering into the 
two-body one is the essence of the integrability ||119|| and has been witnessed in all 



integrable models. The ± sign appearing in ( 7.89| ) is determined by the parity of the 



permutation. In particular case of m = 2 we have derived this formula with — sign 



in (|7.83|) . The periodic boundary condition, as explained already, gives the equations 



for determining the parameters pj, j = 1, . . . ,m generalizing ( [7. 861 ) as 

^ipaN ^ (7.90) 
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accounting for the total phase change incurred in scattering of particle at Xj through 
all the rest round the chain and returning to its original position. For considering 
the thermodynamic limit it is convenient to take first the logarithm of ( |7.9(JD to get 

m 

p,N = 2nQj + I 9{pj,pk). (7.91) 

k=l 

where Qj = half integers for m even and integers for m odd. The expression for 
9{pj,pk) is given as ( |7.83| , |7.8^ ) by replacing 1, 2 by j, k or equivalently as 



_^ I -Asini(pj -pk) 



e(pi,pk) = 2tan-i — ^ (7.92) 



7.3.2 XXX spin-| chain 

All the formulas for solving the isotropic spin model can be derived easily by just 
putting the anisotropy parameter A = 1 in the above equations. This simplifies 
considerably the corresponding expressions giving for example 

™ 1 / 1 1 \ 

E::^ = 2 5:(1 - cosp,), r^^(p„p.) = 2 tan-i - [cot -p, - cot (7.93) 



etc. with the same ansatz (|7.88| ) together with ( [7.891) . 



7.4 Interrelation between algebraic and coordinate Bethe 
ansatzes 

Though the ABA and the CBA differ considerably in their approaches, with ABA 
dealing with global objects having more universal and model-independent formula- 
tions, while CBA being concerned with concrete Hamiltonians and their local actions, 
there exist also intimate connections between them. Though the above results ob- 
tained for the same models using both the methods look different, we intend to show 
here their equivalence basically through a change in variables. Indeed the spectral 
parameter related rapidity variables \j involved in the ABA may be linked to the 
momentum parameters Pj of the CBA by the relation 

e«=/(A,).^= !i4^. (7.94) 

h[\j) sm(Aj - f ) 

for the trigonometric models like XX Z spin— ^ chain. Similarly for the rational 

1 A + — 

models like spin— | XXX chain, one gets the relation e*^-' = f{\j) = j-zi more 
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simply as Xj = | cot y- With mapping ( [7.94|) between the variables we derive from 
(|7.92| ) interesting relations for the scattering amplitudes as 

sm{Xj-\k-j) Xj-Xk-j 

for XXZ and XXX models, respectively, concluding in general e'^'-^^'^*-' = 
Note the remarkable fact, that the scattering amplitude becomes a function only of 
the differences in the rapidity variables Xj and is linked with the function /(A — fi) of 
the i?-matrix. Therefore function /, which is also relevant for the eigenvalue problem, 
is the most important element in the i?-matrix. 



Using the basic relationship ( [r.94| ) between the variables, we can easily verify 



the equivalence between the eigenvalues (|7.19|) and ( [7.87] ) obtained for the same 
XXZ model through ABA and CBA methods, respectively. Similarly one checks 
the equivalence between (|7.22 ) and ( 7.93 ) for the XXX model, which also gives a 



closed relation E!^^ = — I]jLi Using the relationship ( [7.94|) , (|7.95|) we establish 
the intriguing fact that the Bethe equation ( |7.16| ) imposed in the ABA approach for 
eliminating unwanted terms by using an involved argument of vanishing residues for 
the eigenvalue A(A) of the transfer matrix, coincides exactly with condition (|7.90|) , 
obtained in CBA as a consequence of the periodic boundary condition. 

The Nested Bethe ansatz applicable for systems with internal degrees has been 
described above in the ABA formulation. CBA approach can also be formulated 
similarly for such systems. We consider below such an example, where the Hubbard 
model is solved through CBA. The ABA treatment of this model remained unclear 
until very recently ||120|. 



7.5 Hubbard model 



One dimensional Hubbard model is an integrable system, for which i?-matrix , Yang- 
Baxter equation and even ABA formulation is now available ||122|, |121|, |120||. We 



have constructed the Hamiltonian (|6.28|) of the Hubbard model starting from its 
Lax operator in sect. 6 . Our aim here is to present in brief its eigenvalue solution 
through CBA ||123| , |53| mainly to demonstrate the applicability of this method also 
to the systems with internal degrees of freedom. 

The model is described conventionally by the Hamiltonian 



hub 



N 

E 



(7.96) 



<T=± 
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with fermionic creation and annihilation operators c)^j, c^i and the number operators 
Ufji = c^j, Cai with spin projections a = ±. The equivalence of ( |7.96| ) (scale it by t and 



put t = — 1) and the Hamiltonian ( |fci.28 ) constructed before can be shown through 



the Jordan- Wigner transformation (|6.23| ) 



Eigenfunctions for the Hubbard model may be chosen similar to the above CBA 
ansatz. However we notice, that in constructing the m = 2 case due to the internal 
indices a, the states | "^ui^xi < X2) > and | '^2i{xi > X2) > unlike the above 
description correspond now to two independent solutions. This is because together 
with coordinates the corresponding internal indices also get permuted. Therefore we 
have to keep track of both the corresponding functions in the eigenvalue equation. 
Considering first the nonnearest neighbor situation xi + 1 < X2, when the interacting 
U term does not play any role, we obtain 

- (^12(3:1 - 1, X2) + ^12(Xl + 1, X2) + ^uiXu X2-l) + ^12(Xl, X2 + 1)) = ^2^12(Xl, X2) 

(7.97) 

and a similar equation for ^^21- Note that the equation coincides with (|7.78|) for A = 0. 



This may be explained by the fact that the action of h^y from the spin Hamiltonian 
is the same as the Hh^h with U = (their operator equivalence can be shown through 
the Jordan- Wigner transformation ||122 . 



The case X2 = xi = x invokes the interacting term and mixes the \E'-functions to 
yield 

-(^i2(x-l, x)+^2i(a;+l, x)+^2i(a;, x-l)+^'i2(x, x+l))+U ^12(0;, x) = ^2*12(2;, x). 

(7.98) 

Both the functions and '^2i{x,y) arise here depending obn the situation x < y 
OT X > y and the consistency demands 

^'i2(x,x) = ^2i(a;,x). (7.99) 



Repeating the CBA argument described before, we put coinciding x's in ( |7.97| ) and 
take its difference from ( [7.98| ) to get 



-(*i2(x+l, x)+^i2(x, x-1 ))+'^2i{x+l,x)+'$2i{x,x—l)—U \I'i2(a;, x) = 0. (7.100) 

In analogy with ( [7.81| ) we make the ansatz 

^i2(xi, 0:2) = Ai2(piP2)e'(^^"^+^^"^) + Ai2(p2Pi)e'(P^"^+^^"^) (7.101) 

for \l'i2 and similarly for \E'2i, with the corresponding coefficients taken as A2i{p2Pi) 
and ^21(^1^2)- The energy spectrum can be derived from ( [7.97 ) as E!^^^ = 
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— 2 X)|=i(cospj), with obvious coincidence with (|7.82| ) for A = 0. The ansatz 
also extracts from (|7.99| ) the relation Ai2ipiP2) - ^12(^2^1) = ^2i(P2Pi) - ^21(^1^2) 
and using it derives from ( |7.1(J(JD the important interrelation between the coefficients 

A.,.AP2Pi) = E S:i^Jlipm)A^'yMP^) (7.102) 
where the scattering matrix may be written in the operator form 

5l2(A)|A=A0~A0 = + 4^12)|A=A0-A0 (7.103) 

with the notation A° = sinp^ and g{X) = (A + i^)^^ 

The periodic boundary condition 12(3^1, 3^2) = ^21(3^1 +N,X2) induces addition- 
ally the relations like e'P^^ A2i{piP2) = ^12(^2^1) = E^i^ S^'^}{piP2)A^y^{piP2), 
which should be diagonalized to determine the equations for p^'s. Note that the S- 
matrix is exactly like the rational i?-matrix (|7.3|) or the Lax operator Ljk{\) |;^^(;^o_;^o 
associated with the XXX spin-i chain. 

For the general m-particle excitations, since the scattering matrix is factorized 
into the two-particle matrices as described above, the periodic boundary condition 
gives the operator form e^'^^^^ A{I) = TjA{I) where Tj = Yik^^j LjkiX'j — A°) posing 
the problem of diagonalizing the transfer matrix = T(A = A^) of a inhomogeneous 
XXX spin- 1 chain, which however has already been solved. Therefore without re- 
peatation we present the final result for determining {pj},j = 1, . . . , n as 

e^^^ = n/(A.-A°), n/(A^-A^) = n , X x ^^-^q^) 

^=1 1=1 s^Wl^s-Atj 

where A^ = sinpj and the function / is in our standard notation and given here by 
/(A) = 



8 Nonultralocal models and quantum integrabil- 
ity 

As we have mentioned, ultralocal systems received more attention and the integra- 
bility theory is well established for them |TB], while parallel progress has not been 
achieved for nonultralocal models and the related theory is still in the stage of devel- 
opment. However based on recent results [jl9|, some progress in formulating such 
a scheme is reported in sect. 2 along with a list of nonultralocal models proposed in 
different contexts, for which the scheme can be successfully applied. 
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Recall that the basic equation representing integrability for nonultralocal models 
is the braided QYBE (|2.45| ) together with braiding relations (|2 .461 - ^^471) , which lead 
to the global equation ( [2.51| ). The underlying algebra is the braided or the quantized 
braided algebra p 



6 37 



Our aim here is to focus on the concrete models listed 
before and show how their basic equations can be derived from the general formulation 
in a systematic way. 



8.1 Braided Yang-Baxter equation for periodic models 

For closer contact with the physical models we incorporate the periodic boundary con- 
dition Laj+p{X) = Laj{X) and define the global monodromy matrix Ta{\) = T]^'^](A) 
for the closed chain with N sites. Notice now, that for deriving equation for the mon- 
odromy matrix Ta{X) along with (p.46| , p^47D one requires also the braiding relation 



like 

Lbj{fi)Z^^^Laj+i{\) = Zi;'^^Laj+i{\)Z~f^Lbj{fJ.)Z^^^, (8.1) 

since L at the extreme ends again become nearest neighbors due to the periodicity 
Li(A) = L7v+i(A). The compatibility of (|8.1| ) with (|2.46|) demands the constraint 
ZabZba = I putting some restriction on Z. Therefore one may continue the global- 
ization beyond equation (|2.51| ) repeating similar steps and using (|2.46| , |2l47| ). Finally 



at j = applying ( |8.1| ), the braided QYBE for the monodromy matrix for periodic 
models can be expressed as 

i?i2(A - ii)Z:,^^T,{\)Z^iTM = Z^^^T^{ii)Z^^'T,{\)R^^{\ - /i). (8.2) 
Note that the Z matrix has disappeared from the equation. 



8.2 Trace factorization and quantum integrability 

For ensuring quantum integrability we need a commuting set of conserved quantities. 
However by taking the trace of (|8.2| ) we arrive at 

tri2 (z^^^Tr{u)Z^^^T^{v)) = tru {z^iT^{v)Z^^^T,{u)) , (8.3) 

which due to the appearance of Zab does not readily allow the trace factorization . 
It is obvious that for Z = I corresponding to the ultralocal models, the above trace 
identity is trivially factorized proving the integrability. For tackling this problem for 
general Z we consider the following factorization procedures. 
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8.2.1 Factorization by k— trace 



The procedure is based on the generahzation of reflection-type algebras [p^]. 

Let T{X) satisfy the quadratic relation 

A2(A, ^)ri(A)5i2(A, fi)T2{fi) = TMC,2{X, /i)Ti(A)Di2(A, ^^)), (8.4) 

and commute with [Ti(A), K2{fi)] = 0, where A, B, C, D are numerical matrices. Let 
operator i^' be a realization of a similar algebra 

/lr2*(A,/i)i^lnA)fii2(A,/i)K*^(/i) = K*^(^)C'i2(A,^)irlnA)Dit(A,/i) (8.5) 

where 

Si2(A,/i) = {{B%{\, A-,^{\,^^) = A-,^-'-{\,^^), Cu{X, f^) = {{c{i{x, f,)y'f\ 

(8.6) 

Then it is not difficult to show by direct calculation, that the k-trace r(A) = 
tr(T(A)) = tr {K{X)T{X)) forms a commutative subalgebra 

[f(A),f(^)] = 0. 

We like to mention that the spectral parameter dependent braided extensions 
of QYBE's for nonultralocal models can be derived through Yang-Baxterization of 
algebraic relations like braided FRT algebra etc. for matrices satisfying Hecke 
condition The procedure is similar to that of the ultralocal case presented in 
sect. 4, only with the drawback that the spectral parameter dependence of the Z- 
matrices can not be achieved in this way. The R{X — /i)-matrices associated with the 



nonultralocal models is the same ([4.15|) as for the ultralocal models, which with the 
inclusion of the parameters = ±1 is given by 

RtrigW = (Si^ -^-^aa ® + sin(A + tar])Eaa ® E^a + Smr]Eaf3 ® E^^) ■ 

a,l3 

(8.7) 

All the models discussed below are associated with this trigonometric i?(A)-matrix 
or its rational form at at g — * 1. The standard form (with all e = +1) is the most 
used one, while the nonstandard solutions (with some e = — 1) are needed usually for 
the SUSY models. Such i?-matrices exhibit the symmetry 

[i?12(M),SiS2] = 0, S = ^eaa0a- (8.8) 

a 

Therefore whenever s commutes with Z, along with the commuting traces trT{X) 
the s-traces tr{sT{X)) also naturally represent a commuting family. This property 
is used for example, in SUSY models for constructing conserved quantities, due to 
physical reasons. 



118 



8.2.2 Simplifying ansatzes for Z 

For deriving commuting transfer matrices from ( |8.2| ), we notice by comparing with 
that 

Ai2(A,/i) = Zi2i?i2(A-^)Z2iS B = Z-\ C = Z^^\ D{X,fi)=R{X-fi). 

For finding the operator K{X), which commutes with T(A) and satisfies ( ^.5|) , the 
simplest possibility isi^'(A) = k{X)I, where k{X) is a numerical matrix with constraint 

h{X)Z2Mf^)Zi2Ri2iX - /i) = RuiX - fi)k2ifi)ZuhiX)Z2,, (8.9) 



where Z = {{{Z ^) \ It is true that equation (8^) determining the condition 



for factorizability of traces is rather difficult to solve for general Z. Nevertheless, 
for some special forms of the Z-matrix listed below, one gets Z = Z simplifying the 
equation. These forms for Z are 

a) the diagonal form , 

b) factorized form: Z = A ® B, 

c) specifically chosen form as Zu = 1 + J2i ® Bi, (Bi)"^ = 0, 

d) Z chosen as rhe i?-matrix: Zu = Rq2i, 

e) Z expressed as Z12 = e^Si"^»®^' with arbitrary invertible matrices Ai, B^. 
It may be checked easily that for all the above ansatzes the condition Z = Z holds 

and ( ^.9|) reduces to 

RuiX - fl)k2{fl)Zi2kl{X)Z2l = ki{X)Z2lk2{fl)ZuRl2{X - fi). (8.10) 



Consequently any solution of A;(A) obtained from (|8.10|) for concrete R, Z related to 



nonultralocal quantum models would lead to transfer matrices r(A) = Tr(k{X)T{X)) 
commuting for different A and that confirms the integrability of the system. 



8.2.3 Factorization by gauge transformation 

As stated, the nonultrlocality may be acquired through gauge transformation, which 
can be used for trace factorization. Let the transforming operator / be a matrix of 
operators acting in the Hilbert space and satisfy the relations 

[i?i2(A),/i/2] = 0, [/i,/2] = 0, Z^,'Ti = l2Til^'. (8.11) 

Starting from ( p?2| ) and using (p.ll| ) the factors Zf^^Ta, a,b = 1,2 appearing in the 
equation can be replaced by IbTJ'j^^. Multiplying from the right by I1I2 and from the 
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left by its inverse and using further from (|8.11) the commutation relations between 



la,R{\) and between /^'s themselves we arrive at 

RuiX - fi)i^'T,{\)hi^'T2ifi)i2 = i2'T2{fi)i2k'Ti{fi)hRi2{X - /^). (8.12) 

One detects the appearance of gauge transformed monodromy matrices T(A) = 
l^^T{X)l and taking trace of ( ^.12| ) finds the required commutation relation 



[trT{\),trT{fx)] = 0. 

8.3 Examples of nonultralocal models and their basic equa- 
tions 

We focus on the examples of nonultralocal models listed in sect. 2 and show that 
the basic equations representing their quantum integrability can be obtained system- 
atically from the set of relations (|2.45| - p^47|) , ( ^.2|) etc. derived before from general 



arguments. These examples involve either spectral parameter independent or the 
spectral dependent trigonometric matrix ( ^.7] ) or its rational limit. At the same time 
they correspond also to different braiding, i.e. different choices for Z, Z matrices. 
Curiously most of these choices for the existing models agree with the above forms 
of Z, proposed on the argument of simplicity. 

8.3.1 Supersymmetric models 

Quantum integrable supersymmetric (SUSY) theory |^9|, |113| , |124|| is possible to for- 



mulate in a convenient alternative way expressing the gradings in a matrix form |[147 
We observe that such a formulation is in agreement with the framework of nonul- 
tralocal models for the choice of homogeneous braiding Z = Z = rj with rj as the 
diagonal matrix 

1112 = ^gapeaa ® 6^/3 with g^fi = (-1)"'^. (8.13) 

a,l3 

Here a = 0(1) is the supersymmetric grading depending on the even (odd)-ness of 
the indices. Due to the obvious property ?7i2^72i = 1, we get from ( |8.10| ) k = 1 and 
integrability follows from the commuting traces. However supertraces, obtained from 
the s-traces ( p.8|) by choosing (pa = (—1)° are usually used for defining conserved 
quantities for such models. The super QYBE's for the SUSY models can be derived 



from the general relations ( p.45| , ^.2|) as 



^12 (A - n)r]i2Lij{\)r]i2L2jifi) = rii2L2j{^i)rii2Lij{\)Ri2{\ - /i) (8.14) 
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and 

Ri2{\ - fJ')r]i2Ti{X)r]i2T2{fi) = ?7i2T2(/i)r7i2Ti(A)i?i2(A - /i). (8.15) 

On the other hand, the nonultralocal braiding (|2.46| , p~4^ ) correspond to the super- 
commutation relations 

L2k (/i)?7l2^1j (A) = ?7l2^1i ( A)r7l2^2A; (/i)'7l2 (8.16) 

for k ^ j. For comparing with the conventional results on SUSY models, we express 
the QYBE (8.14) and the braiding relation (|8.16|) in matrix elements to give 



(8.17) 

and 

(LlliXML'^W), = (-l)('^^+^^)('^+^^H^aUA)),(^^^,(/i))., (8.18) 

respectively. Note that this reproduces exactly the results of |113|| and [p.24|| with 
i?-matrix as rational ( ^.21|) or the trigonometric solutions (|8.7|) of nonstandard type 
with proper choice of ea = ±1 depending on the specific gradings. 

8.3.2 Integrable anyonic models 

Quantum integrable anyonic models (JAM) can be introduced by generalizing ( |8.13| ) 
in the nonultralocal description of the SUSY models by choosing Qajj = e*^ 
where 6 is the arbitrary anyonic phase and d = 1 is the anyonic grading, while 
a = gives bosonic commutations. For such a Z = Z-matrix the braiding relation 
describing commuting relations between Lax operator at different lattice sites k > j 
is given by 

Lt:^kM)Ltl^j)W = e-^^(--^^)(^-^^)L-(^.)(A)L- ,)(/i), (8.19) 

with all matrix indices running from 1 to N = m + n. The algebraic relations at the 
same point / on the other hand are given by the corresponding braided QYBE 

{k} {k'} 

(8.20) 

The quantum i?-matrix may be chosen in the trigonometric form ( p.7| ) or as its 
rational limit 

i?i2(A) =ti2A + -Pi2, (8.21) 
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where P is the permutation matrix and t = I for the standard and ti2 = ?7i2 for the 
nonstandard solutions. The consistency equations for Z and R matrices obviously 
hold for all these solutions. A possible Lax operator of the anyonic model may be 
proposed as 

^^(0(A) = A5.,pr + |e^^(%S- (8.22) 
Using rational i?-matrix ( B.21|) one may derive from (^M ) and ( FTl) the anyonic su- 



per algebras by matching the coefficients of different powers of the spectral parameter 
For standard and nonstandard solutions of i?-matrix these algebras also differ. 



They are nontrivial generalization of the well known graded super algebra ||6^, where 
together with an anyonic parameter 9 a set of additional N operators appear. 
Nonstandard i?-matrix and the choice 6 = ir recovers the known super algebra. 

For constructing quantum integrable models involving m number of bosons fo^^-* 
and n number of anyons fj^''\ one may consider a realization of the generators as 

Pl? = e^^\ pL'I = f:!^''' bi'\ Pi'^=bi^''e\ v'S = n'''f?- (8-23) 

The bosons satisfy the standard commutation rules and commute with anyonic op- 
erators, while the anyonic operators satisfy the commutation relations 



for k > j and 



and 



\pT\ f's''] = 0, /i') - e^' ft''' = (8.25) 



pTfS'=e-'Ui;'pf\ [/«, f?] = 0, (8.26) 



at the same lattice points. Note that these relations are consistent with the known 
anyonic algebra fT^], though defined in 1-dimension. Even though explicit construc- 



tion of integrable anyonic model has not been achieved yet, Kundu-Eckhaus equation 
considered below is however an interesting example of a similar realization. For dis- 



cussion on various physical aspects of the anyon theory we refer to [125 



By choosing trigonometric i?-matrix solution and the same Z-matrix, one may 



obtain a g-deformed anyonic algebra |21|, which represents a generalization of the 



extended trigonometric Sklyanin algebra to include anyonic parameter 6 and 
nonultralocal braiding. 
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8.3.3 Integrable model on moduli space 

In considering Hamiltonian Chern-Simons theory, moduli space of flat connections 
on the Riemann surface with marked points appear. From monodromies of flat 
connections along fundamental cycles of the Riemann surface an integrable model 
may be constructed (listed as IMMS in sect. 2), with spectral parameter dependent 



Lax operator |126 



Ln{X) = K" + XL (8.27) 

The if" matrices are defined as 

= (Lr)-^ ■ • • {L-_,)-'K,,L-_, • • ■ Lr, (8.28) 

by setting Ki = {L^)^^Lf, where L"^ satisfy the FRT relations ( p.76|) . Using the 
definition ( |8.28| ) and the equations for the equations satisfied by the Lax operator 
{^2T\) can be shown [12|, [12^ to be 

S{\, fx)L,n{\)RL2n{f^) = L2nif-i)RL,n{\)S{\, fX, ) (8.29) 

where S{X,fi) and S{X,fi) = {R)~^S{X, fi)R are solutions of YBE ( p.l3| ), along with 

RL^U^)R-^L2n{fi) = L2n{li)RLim{X)R~\ Til > u (8.30) 

These equations induce for the monodromy matrix T(A) = Hi Li{\) the relation 

S{\ ^l)T^{\)RT2{^l) = T2{^x)RT^i\)Si\,^l) (8.31) 

giving commuting transfer matrices 

t{X) =trg{T{X)) = tr{MT{X)), M = diag{q-\q) (8.32) 

for different values of spectral parameters and thereby producing a family of conserved 
operators in involution. 

This model can be shown to be well described by the general nonultralocal theory 
presented here, for which we identify S{\,fi) as the trigonometric Rtrig{\ — jj) and 
choose Zi2 = Z21 = Ri2- This clearly rewrites the general BQYBE ( p.45| ) into ( |8.29| ) 
and the braiding relation ( p.47|) for different points into ( ^.301) . The braided QYBE 
(p.51|) for the monodromy matrix related to nonperiodic models coincide therefore 
with ( ^.31| ) recovering the earlier results. We discuss below examples of models, which 
unlike the homogeneous long range nonultralocality allows braiding only with nearest- 
neighbors. First we look into some models described by the spectral parameter 
independent i?-matrices with equations like ( |2.48D and related to the conformal field 
theory, without any demand on trace factorization. 
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8.3.4 Current algebra in WZWN model 

The WZWN model is described by the unitary unimodular matrix-valued field 
g{x, t) : M ^ su{2), M = x R^. Defining the chiral left current as L = \{Jo + Ji), 
where = d^gg~^, one gets the well known current algebra relation ( |2.52| ). Intro- 
ducing dxU = Lu one can write the monodromy as Ml = u{x)~^u{x + 27r) and 
the similar relations for the right current. With the aim of unraveling the quantum 
group structure in WZWN model a discrete quantum version of the current algebra 
was formulated in [^. We demonstrate that the basic equations found in can 
be reproduced by choosing Z12 = Rq2i with Z = I'm. the spectral parameterless limit 
( |2.48| ) of the braided equations as 

-R^12(-R^12) ^^Ij^tj = (-RJl2) ^-^2j-^1j(-R^12)' 

which considering the relations R^i2 = (R^u)'^ = Rq2i with q = q^^ and Z21 = R^u 
gives 

(8.33) 



where we have suppressed the + symbol in Laj. Similarly the spectral independent 
limit of the braiding relations for the present choice of Z, Z reduce to 



L2j+lR^l2'^lj — ^Ij^ 



li-^2i+l, 



L')hL 



2k-l^lj 



(8.34) 



for k > j + We observe that the nonultralocal relations (|8.33| - ^^34D coincide exactly 
with the current algebras obtained in . Identifying further the chiral components 
as Uj = T'-^'^I, we get from ( |2.51D a relation analogous to ( ^.331 ): 



-^512'"lj'"2j-Rgl2 — '^2jUlj 



(8.35) 



Since 



Ur 



Ln...Lj^iUj, multiplying (|8.35|) from the left successively by 



^2n 



^2j + l; J^2j+2^ ■ 

35]) the algebra at n > j as 



and using the braiding relations ( |8.34|) we can derive from 



UlnU2jRgi2 = U2jUin (8.36) 

Similar exchange relations can also be derived for Por monodromy matrix 

Ml on the other hand, we rewrite equation (|8.2|) suitable for the periodic models: 



M^{R^,2r'M^R^,2 = {Rt2)-'M^Rl2M^ (8.37) 
We notice again that the relations (|8.35 - p37 ) recover correctly the results of |61 . 
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It is checked easily, that the present Z\2 = Rq2i satisfies the conditions (|2.14|) and 



(|2.49| ) for Ri2 = Rqu and since it corresponds to the trace factorization case d), the 



model might serve as a candidate for constructing nonultralocal quantum integrable 
models. Moreover, if the long range braiding given by Z is chosen nontrivially, it 
might also lead to a generalization of the current algebra of the WZWN model. 

8.3.5 Braided algebra 

As described in sect. 3 the braided algebra (BA) may be given in the compact form 
by the matrix relations (|3.109|j3^110|) , where the elements of the matrices T and T' are 



two copies of the generators of BA. It is easily checked that by a homogeneous choice 
for the braiding as Z12 = Zyi = R21, one can generate from the same spectralfree 
BQYBE the braided algebra relations (|3.109| , |3.110|) . 



8.3.6 Quantum group structure in Coulomb gas picture of CFT 

The Coulomb gas picture of CFT (CG-CFT) is based on the Drinfeld-Sokolov linear 



system ||6^ dQ = L(x)Q, given by the Lax operator ( p.53| ) involving periodic field 



P{x) with nonultralocal bracket. In a lattice regularized description was given 
at the quantum level. We show again that the same result can be obtained from 
the general relations formulated here for the nonultralocal models at the spectralfree 
limit considered in the above example. 

We choose Z = 1 with nontrivial braiding Z12 = q^i^^®^\ = en — Cj+ij+i, in 
agreement with simplification e), which reduces ( p.9| ) to [R,C] = 0, C12 = ^^2^1 ^2- 
Noticing that, Hn + Hi2 commutes with the parameterfree i?^-matrix, we may 
extract a solution k = q^i^^, since C12 = 2'^2^i^2 = gSi^^'i"'"^'^)^^ Defining 
Laj = B~^Laj, B = we can eliminate the Z dependence from the parameterless 
BQYBE yielding 

-Ri2-^ij-^2j = -^2jLijR^2 (8.38) 
The braiding relations on the other hand give 

L2j+iAi2Lij = LijL2j+i, L2kLij = LijL2k, k > j + 1 (8.39) 

where A12 = Zi2. These equations recover the related findings of [^, while their 
algebras of monodromies can be found from ( p.51| ) and ( |8.2| ) at the parameterfree 



limit, by defining Qj = B ^T^^'^'^ and S = B as 

Rt2QijQ2j = Q2jQi,Rt2 (8-40) 
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and 

i?^2'S'l^l2'S'2 = S2A12S1R12. (8-41) 

The following examples are of genuine integrable models with trigonometric or 
rational i?(A)-matrix and spectral parameter dependent Lax operators. 

8.3.7 Nonabelian To da chain 

NATC is a one-dimensional evolution model on a periodic lattice of sites described 
by the nonabelian matrix- valued operator gk G GL{n). The model, representing a 
discrete analog of principal chiral field, was set into the Yang-Baxter formalism in 



with the nonultralocal Lax operator ( p.54|) . The associated i?- matrix is the rational 



solution -R(A) = ihP — A with P = 11 ® tt, where n is a 4 x 4 and vr a x 
permutation matrix. We show that the model can be described well by the general 
formalism with the choice Z = I and Z12 = 1 + ih{e22 ® 612) ® vr. This is an example 
of models with acquired nonultralocality through gauge transformation. Using ( |2.54| ) 
one shows that Z2iLij\[ = /2-^ijv^V^ with / = diag{I , g^q) and derives Z21T1 = l2Tii2^ 
and other conditions of ( p. 111 ). Therefore the quantum integrability of the model 
follows from trace identity ( ^.12| ). With this input for R^Z^Z our main formulas 
(p.45| , |8.2|) recover the QYBE's of while the braiding relations are given by 



L2,+i(/i)Z2-/Li,(A) =Li,(A)L2,+i(^), L2fc(/i)Li,(A) =Li,(A)L2fc(^) (8.42) 
for k > i + 1. 

8.3.8 Nonultralocal quantum mapping 



The quantum mapping '^j is related to the lattice version of the KdV type 

equation 

^2i-l = ^2j, V2j = V2j+1 , 2^ V2j = 2^ V2j+1 = C. (8.43) 

The generalization of integrable NQM is associated with the Lax operator L„ = 
V2nV2n-i, with K„(A) (|2.55|) involving hermitian nonultralocal operators Vij (|2.56|) . 
The i?(A)-matrix is given by the rational solution 

h P 

i?i2(Ai-A2) = / + — ^ (8.44) 

Ai — A2 
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with the braiding 



h 

Zl2{\2) = 5l2(A2) ^ = 1 + ^ ^ CTVa ® e«jv (8.45) 

^2 a 



and similarly for Z2i{Xi) along with the choice Z = 1. This is an interesting example 
for the nonultralocal scheme, since it is the only model where the braiding matrices 
are spectral parameter dependent. We can see that the choice coincides with the 
simplifying ansatz c) presented above and therefore trace factorization is possible with 
k = 1. Consequently trT{X) constitutes the commuting family along with tr{sT), 
where Sa = diag{l, 1 + ^) We should stress here on a subtlety indicated in sect. 
2, that the present Zi2{X2) does not satisfy the standard YBE (p.l3|) , while the same 
equation holds for the matrix i?i2(Ai, A2) = (5'i2(A2))^^-Ri2(Ai — A2)5'2i(Ai) . 

Putting Z = I'm the braided QYBE's ( |2.45D and introducing -R(Ai, A2), we derive 



therefore 

i?i2(Ai, A2)Li,(Ai)L2,(A2) = L2,(A2)Li,(Ai)i?i2(Ai - A2), (8.46) 
and from the braiding relations ( |2.46| - p747D the nonultralocal condition 

i^2,+i(A2)52i(Ai)Li,(Ai) = Li,(Ai)L2,+i(A2) (8.47) 

and the trivial commutation L2fc(A2)Ivij(Ai) = Li^ (Ai)L2j!c(A2), for k > j + 1. 

The equation for the monodromy matrix with the periodic boundary condition 
can similarly be obtained from (|8.2|) in the form 

i?i2(Ai - A2)Ti(Ai)Si2(A2)T2(A2) = r2(A2)S2i(Ai)ri(Ai)i?i2(Ai - A2) (8.48) 

Thus all the basic equations of NQM ([0) can be obtained systematically from the 
general nonultralocal scheme. 

8.3.9 Reflection equation 



The reflection equation (RE) due to can also be derived from the nonultralocal 



point of view, by considering a monodromy matrix arranged as T(A) = Jlj Lj{\) 
such that 

r V^^(A) forj>iV 
Lj{\) = \ K~{\) for j = N (8.49) 

forj<iV 

with Lj{X) being any ultralocal Lax operator. As a result the local QYBE becomes 
the standard one: -Ri2(A — fi)Lij{X)L2j{fi) = I/2j(/i)I/ij(A)-Ri2(A — /i) with the 



127 



braiding 

-^^2(2Af-j)(/i)-Rl2(A + yU)Zl(j + l)(A) = Li(j + i)(A)i?i2(A + n)L2(2N~j)if^)- (8.50) 

Due to the braiding relation ( |8.50| ) and the equation for K^{X): 

i?i2(A - /i)Kr(A)i?i2(A + fx)K2ifx) = K2{^i)Ri2{\ + fi)K^{\)Ri2{\ - ^i). (8.51) 
the global QYBE takes the form 

i?i2(A - ^)Ti(A)i?i2(A + ii)fM = T2(/i)i?i2(A + /i)Ti(A)/?i2(A - /i), (8.52) 



and recovers the reflection equation of |63[. It is interesting to note that, though it 
refers to a braiding different from that considered here, the choice Z = 1, Zyi = Z21 = 
-^12(^1 + W2), with symmetric i?- matrices derives formally the reflection equation 
(|32D from BQYBE ( p3ID . 



8.3.10 Kundu-Eckhaus equation 

Higher order nonlinear equation KEE (p.60|) can be derived from the NLS equation 



by gauge transformation |8^. However the quantum model becomes nonultralocal 
with the field ipn = An+iipnAn, An = e~*^^^»=i^'^»^'^'^ exhibiting anyon like algebra 

^n^m = e^''-^tpm'4'n, foT U > Til + 1, tpm+l^m = e''^ (8.53) 

along with the standard relation [ipn, i^U = at the same point. Note that the 
braiding relations are inhomogeneous and different from the pure anyonic relations. 
The lattice regularization of the Lax operator (|2.61| ) can be obtained as L„ = 



hn+iL^n^^^^hJ" by considering the gauge transformation = A'^ over the lattice 
NLS Lax operator L(,'"'^) ( glSOl) . 



Using the commutation relations of tpn the braiding relations can be derived as 

L2j+l{fJ')Lij{\) = Lij{\)Zi2L2j+l{fi)Z^2 y -^2fc(/^)^l2^-^lj(A) = ^12^-^li (A)2'i2//2fc(/^)^12^ 

(8.54) 

with k > j + 1 and the braiding matrix Z12 = Z21 = diag{e^^, 1, 1, e*^). For periodic 
BC on ipn, the KEE model does not show periodicity :L^^i = h]^_^^LihN and we 
have il^N^^i = e^'^^ipi'ipN' giving 

L2^(/i)Zf2'Ln(A) = Z^,^LniX)Z^,2L2Nifi)Z^2'- (8-55) 

The braided QYBE with rational i?(A)-matrix has the simple form 

i?i2(A - n)Lij{\)L2j{n) = L2j{fi)Lij{\)Ri2{\ - fx), (8.56) 
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since [-R(A), Z\ = 0. Observe that the braiding for this model, differs from the types 
assumed for the nonultralocal scheme. Since the braiding in general can be widely 
different, certain cases naturally might go beyond any particular class of choice and 
should be treated individually. The KEE model due to peculiar inhomogeneous braid- 
ing also gives trouble in achieving trace factorization. By defining the monodromy 
matrix in a special way T(A) = Ljv(A) . . . L2(A)I/i(A), where I/i(A) = Z^^Li{X)Z 
one may derive the global BQYBE 

R{X - i2)T{X)ZT{iJ,) = T{iJ,)ZT{X)R{X - (8.57) 

However, since T = T12 is defined now in both the auxiliary spaces, it is difficult to 
factorize the trace in tr{T{X)ZT{n)) = tr(T{^)ZT{X)), obtained from ( p. 571 ). 



8.3.11 Quantum mKdV model 

Modified KdV equation is a well known classical integrable system with wide 
range of applications ||128|| . The quantum mKdV (QMKD) model is important for 
its relation with the conformal field theory [|7^. However, due to its nonultralocal 
algebraic structure {v{x),v{y)} = ^26'{x — y) the quantum generalization of the 
model could not be done through standard QISM [0, but was achieved following 
the scheme for quantum nonultralocal models and presented here . 

A lattice regularized Lax operator involving a single field ||129|| had a problem 
of mismatching the number of fields, which was overcome by introducing the Lax 



operator []65 



(8.58) 



with two different quantum operators W^. 
tralocal commutators 

ft 



e 3 . 



The operators ivf satisfy nonul- 



and 



4, 



] = 'ilT^ih-u - 24,« + 4,/-i) 



which go at the continuum limit A ^ 0, to 



^.59) 



(8.60) 



±^72 (4(a; - y) - Sy{x - y)) 



±ijhS'{x — y) 



(8.61) 



with discrete operators turning to the quantum field Av (x) and dk-^ — >■ 6x{x—y). 
The reduction of operators ^ I + iAv'^{x) at the continuum limit yields from 
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(^381) Lk{v~,v+,^) ^ I + AC{v-,^) + 0{A'^{v+)), where Cjir, Q men by (|237|) 
is linked to the well known mkdv Lax operator |]13[ as cr^£(w~,,^) = U{v{x),^) . It 



is remarkable that, though the lattice Lax operator ( |8.58| ) contains both variables 
having nontrivial commutation relation ( ^.601 ), dependence of one of them drops out 
from the continuum Lax operator (in the present case it is v~^). 

Inserting the lattice Lax operator (|8.58|) in the braiding relations (|2 .461 , ^74^ ) 
and using commutators (|8.59| , |8.60|) , we derive the matrices = Z21 = 



q~'^ ^'^ and Z = I. Plugging the Z-matrices in ( p.45| ) one finds that the R{X)- 
matrix of the lattice mKdV model is the same trigonometric solution related to the 
XXZ-spin chain. Since this i?(A)-matrix commutes with the present Z, the braided 
QYBE's simplify to 



and 



RuiX - /i)Ly(A)L2j(/i) = L2j{fi)Lij{X)Ri2{X - fi) 



Ri2{\ - /i)Ti(A)Zr2'T2(/i) = T2{i^)Z{^'T^{X)Ru{\ - fi), 



(8.62) 



^.63) 



for the periodic case. 

For identifying conserved quantities from ( p.2|) , one finds easily that the k = 1 
is a trace factorization solution giving [trT{X),trT{fi)] = . Similarly due to 
[Ri2,SiS2] = 0, the s-trace Tg{X) = tr{sT{X)) = q-^ A{X) D{X) with s = q"^''', 
also generates the commuting family of conserved quantities. Exact solution of the 
eigenvalue problem of this model will be considered below. Recently some other forms 



of quantum lattice KdV model have also been investigated in a different way ||8^, pq] . 

Thus we have seen that nonultralocal quantum models proposed in different con- 
texts at different times and solved in individual ways, can be unified under one 
scheme. All their basic equations can be derived systematically starting from the 
same braided QYBE and the related braiding relations for different choices of Z, Z 
and /^-matrices (see fig. 8.1). Therefore, much in parallel to the QYBE for the 
ultralocal models, its braided extensions can be considered to be some representing 
universal equations for describing integrable nonultralocal models of certain class. 



8.4 Algebraic Bethe ansatz solution for the quantum mKdV 
model 

To demonstrate that the Bethe ansatz technique described in sect. 7 is applicable also 
to nonultralocal quantum models, we present here the ABA scheme for the quantum 
mKdV model [Q. For solving exactly the related eigenvalue problem ( |7.4| ) including 
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that of its energy spectrum, one has to follow the standard ABA formulation, keeping 
simultaneously the track of nonultralocal signatures. For calculating the eigenvalue 
of the transfer matrix ( |7.1ij| ), the eigenvectors are taken in the usual form ( |7.t)| ), while 
the commutation relations of A{X) and -D(A) with B{n) should be obtained now from 
the braided QYBE ( |8.63| ) using R, Z-matrices as 

A(A)i?(/i) = q-'f{f,-X)B{f,)A{X)-q-'f{f,-X)A{f,)B{X) 
D{X)B{fi) = q'fiX- fi)B{fi)D{X)-q'f{X- fi)D{fi)B{X), (8.64) 

where the factors coming from the i?-matrix are f{X — fi) = '^"^J^^^^^^ and /(A — 
n) = ^if^;z^- Note that the multiplicative factors g^^ appearing in ( |8.64| ) are the 
contribution of Z due to nonultralocality. 

It is evident from the form of ( |8.58| ), that like the case of SG or Liouville model 
(|7.53|) , instead of a single Lax operator one has to use = LkLk+i for constructing 
the pseudovacuum. This gives the pseudovacuum ||65[ as 



|0>= n |0>^' \0>k=S{W^_,,-qW^)S{W^^, + {W,-)-'). (8.65) 
(^=1) 

and from the action of \ >k, 

a{X) = 2ze'(^+^) sin(A + ^), = e'^ (8.66) 

and 

/3(A) = 2«e'(^-^) sin(A - ^—). (8.67) 

Similarly one derives the crucial triangular form for the Lax operator acting on the 
pseudovacuum due to Lfc | >k= 0. In analogy with the standard ABA, the eigen- 
value problem ( |7.13|) for the mKdV model may therefore be solved by following the 
relations ( [mi , [rT^ [TT^ ) and using (|]M|,^]T3|) as the s-trace 

A(A) = g-ta(A)^n(9"'/(A,-A))+gt/5(A)^n(?'/(A-A,) 

3 3 

= 2.e-[g--f sin-(A + ^) ft 

2 fj^ sm(A-Aj) 

+ g-(-f -^-) sin^(A - ^) ft ^^"(.^ - + :;^^ ]. (8.68) 

2 ^-^^ sin(A — Xj) 

The energy eigenvalue can be obtained by expanding ( ^.68|) in the spectral parameter 
and considering C_3 = H as the Hamiltonian of the model. The determining equation 
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for the rapidity parameters Xj may be obtained from the vanishing residue condition 
of the eigenvalue equation ( p.68| ) in analogy with the standard Bethe equation ( |7.16| ) 



as 

Remarkably, these results for the quantum mKdV model, apart from an extra q 
factor, coincide exactly with the XXZ spin-^ chain. It is true however that, unlike 
XXZ the Hamiltonian of the mKdV model is H = C_3 and moreover one has to 
consider also the continuum limit for the mKdV field model. 



fig.8.1 



Figure 8.1 Classification of the nonultralocal models of spectral parameter de- 
pendent or independent types belonging to trigonometric and rational classes along 
with the respective choices for Z and Z matrices. 

9 Recent progress and unsolved problems 

We briefly touch here on some recent progress in the fields related to the quantum 
integrable system or influenced by it, which can be followed from the material pre- 
sented in the review. We also indicate for each topic some unsolved problems, which 
are either long-standing problems or have been opened up along with the recent 
development of the subject. 

9.1 Higher rank generalization of the ETSA and construc- 
tion of multicomponent models 

We have seen in sect. 5 that through different realizations of ETSA wider class 
of models can be reached, than those constructed from the quantum algebra. The 
Uq{sl{n + 1)) algebra on the other hand has been shown to generate multicomponent 
discrete Toda field models. Therefore, it may be hoped that the higher rank general- 
ization of the ETSA would lead to multicomponent models with richer structure of 
their Lax operators. 
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Such a generalized ancestor Lax operator related to the trigonometric Rtrigi^)- 
matrix ([4.15|) was constructed in |Q in the form 



N 1 N 1 

= E(7 Etk + E 7 ^ik Eh + e m El,. (9.1) 

k >• k,l ^ 

, Tki are generators of the higher rank ETSA given by the relations 
TkTki = e TkiTf, , T;. r^fc = ^/fc^^fc , 

[th ,r,fc ] = 2isina ( t^t+ - T+Tf ) , (9.2) 

for all different k,l,m,n indices. The diagonally placed operators commute 
among themselves in addition to 

[ Tmk , Tki] = 2i e sin a Tk^^^Tmi , (9.3) 

where e = sign(/i; — + sign(/ — m) + sign(m — k) , r^^'* = for e = ±1 
and 



[ Tkl , Tmn] = 2i p slu « TmlTk„ , (9.4) 

with p = +l for I > n > k > m ( and all its cyclic inequalities ), while p = — 1 
for the reverse inequalities and p = , otherwise. 

It may be noted that the above relations are the generalization of the ETSA 
(|5.2|) and reduce to it at N = 2 . On the other hand for a symmetric reduction 
through generators of the quantum algebra, (|9.1|) reduces clearly to ( [4.16|) producing 
generalized spin model and Toda field model. Therefore as descendants of and 
its multiparameter generalization one can expect to generate models with asymmetric 
form of their Lax operators like multicomponent DNLS, MTM, RTC, ALM etc. 

At the limit q 1, ( |9.1D reduces to 

^ i\ ^ 

L = Y.{Kt + -Kr)Eu + E Ki,E,u (9.5) 
I ^ j^i 

and is associated with the rational (N^ x N^) i?-matrix ( [4.32| ). The operators K 
satisfy the algebra 

[Kr^k, Kki] = K^K^i, [Kki, Kik] = K+Kf - K^K+, 
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(9.6) 



along with [K^^Ki^n] = [KkuKkm] = Kki,Kmi] = [Km, Kmn] = 0, where Kj^ com- 



mute with all other generators and form an abelian subalgebra. Notice that (|9T 
is an higher rank generalization of (|5.52|) and recovers the known sl{N) at particular 
symmetric reduction. 

Different bosonic realizations of this algebra and its multiparameter extensions 
should be able to generate from (|9.5| ) quantum integrable models like multipar- 
ticle Toda chain or Toda chain like models, as well as the vector generalization 
of the SLNLS, generically not exhibiting unitary symmetry. Let us take up the 
last example where through a set of independent bosonic operators with 

bPi, (pk] = Sik, [ipi, i^k] = [4>h 4>k] = 0, the algebra ( |^ may be realized as 



Kij = ipj = ij, Kji = 4)j = 0*, Ki. 



0, 



N) 



0, l<ii,j)<N. 



The corresponding Lax operator ( |9.5D then reads 

m = 



I 



(9.7) 



(9.^ 



which represents a simple lattice version of the vector NLS model, generalizing ( |2.4CI| ) 
to the vector case. 

It would be fruitful to construct and analyze the algebraic structures of the various 
quantum integrable multicomponent models through higher rank ETSA and solve 
them through the nested Bethe ansatz. The quantum Dodd-BuUough-Mikhailov 



model 1 130] might be a good candidate for such analysis. 



9.2 Integrable model with ladder symmetry 



The ladder symmetric models are of immense current interest due to their possible 
connection with the high Tc superconductivity ||149|| . Theoretical models with ladder 
symmetry exhibiting quantum integrability can also be introduced [^. The idea of 
constructing XXX spin-| chain with inhomogeneity and its Bethe ansatz solution 
(described earlier), can be extended to the tensor product of similar Lax operators: 



Lab{X) = + T])U2{X)Lal{\)La2{\ - V) , 



(9.9) 



where Li2{\) = ^ is equivalent to Rrat{^) (|4.32|) associated with the XXX chain. 
Using now the method of Hamiltonian construction of sect. 6 and noticing that 
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Lafc(O) = R^2iPb2i+iPa2iRa2i+i^ wheie R^2i = Rb2ii±v), arrive at the Hamihonian 

= (0 X! -^2i,2i+el(-R^)2i!2j+el + (-R'')2i,2j+el-R2i,2i+el-P2i,2i+el(-R'')2i!2i+el 

22-15*22 + S2iS2i^l) + J2 y^XS2i-lS2i-\-l + 5'225'22+2) 

i i 

+ h Y.iS2^{S2^+l X ^+2) + X 52,) (9.10) 

i 

where Ji = J2 = ffJi^ J3 = —'^V^i ■ Note that the lattice on which the above 

Hamiltonian is defined has ladder symmetric form, which has important physical 
consequences ||149|| (see fig. 9.1a,b). 



The QYBE for Tab = Ti,{X + ri)Ta{X), where Ta{X) is the monodromy matrix of 
the inhomogeneous model constructed from the Lax operators 

Ra^{X) = Ra2i{X)Ra2i+l{X " 7]) (9.11) 

yields trab(Tab{X)) = t{X + ri)t{X) commuting for different A, where t{X) = tr(Ta{X)) in 
turn is the commuting transfer matrix of the inhomogeneous XXX chain. Therefore 
the eigenvalue of tab I >= A(A) | m > is given by A(A) = A°(A + ri)A^{X), 
where A°(A) is the eigenvalue of t{X). Due to the tensor product structure of the 
Lax operator (|9.11| ) the results for the standard XXX spin chain ( [7.20| , |7.21|) are 



generalized to the expression 



A°(A) = n + /?"(A) n " (9-^2) 

j=l - ^ - 2 j=l -^j + 2 



for the eigenvalue along with the Bethe equation 



- br = n 7 \\ \ (9-13) 



where (3{X) = (^) (3;^) • As a result similar to (^2^) we obtain the excitation 
spectrum E.^ = i^lnA(A) |a=o= Eje°(Aj) + e°(Aj +r]), where e°(Aj) = -tjtt- 

j i 

The thermodynamic limit of the result and other details can be found in []77| 

This two-chain ladder model can be generalized to arbitrary k-chains. It is conjec- 
tured that the ladder models with even (odd) number of chains must show presence 
(absence) of mass gap in their excited spectra. Whether such integrable ladder models 
can throw any light on this conjecture is an important open problem. 



135 



9.3 Elliptic quantum group 

The quantum group associated with the trigonometric i?-matrix and generated by the 
QYBE ( |2.9|) has been discussed in detail in this review. In recent years there are some 
results to extend this notion to the elliptic case |T^. Such an algebra Ak^qigli^)) 
may be generated by an extension of the QYBE as 

^^2(6/6)^1(6)^2(6) = ^2(6) Li(6)i?l2+(6/6), (9.14) 

where the quantum indices of L are suppressed, R*^{^) = R^{^, k*^ , q^), k*^ = 
k^q^'^ with -Ri2(0 being the elliptic solutions like ( ^.62[ ). The complex parameters 
qh and ^ may be interpreted as the elliptic nome, the deforming parameter and 
the spectral parameter, respectively. Here for the central element g"'^ = g~' , the 
algebra is called level /-algebra with its simplest representation having I = 0. Impos- 
ing the constraint on quantum determinant detgL{^) = gt the algebra Ak,qisl{2)) 
may be obtained. Note that in contrast to ( p.9|) the relation ( p.l4| ) contains two dif- 



ferent R, i?*-matrices, which makes the introduction of the coproduct structure and 
hence the Hopf algebra difficult. For obtaining the algebra from ( p.l4| ) (for details 

see 0) the matrix elements ^^,'(0, ee' = (++), ( ), (H — ), ( — h) of L{^) should 

be expanded as i^«'(0 = I^nez -^ee'.n^"" • Unlike the trigonometric case this gives 
expansion like 

L+(0 = L{qH) = --- + k-^ Ltr' + + Lt,^ + ■■■ (9.15) 

mixing both +ve and —ve powers of ^ and similarly for L^(^). 

At A; — > 0, the elliptic Reiup —>■ trigonometric Rtrig matrix and L~^{C,) , as seen from 
( |9.15| ), contains only +ve powers of ^. (Similarly L~{C,). contains only —ve powers). 
Consequently the elliptic quantum group reduces to the well known trigonometric 
quantum group. 

Discovering the universal 7^-matrix for the elliptic quantum group, establishing 
its Hopf algebra structure and finding its relation with the quadratic Sklyanin algebra 
discussed in sect. 3 are some challenging open problems. It would also be interesting 
to find physically relevant integrable ultralocal and nonultralocal models related to 
such elliptic algebra in the spirit of the present review. 

9.4 Coupled spin chain and extended Hubbard model 



Using the idea of twisting ( [4.28| ), one can construct an integrable coupled spin chain. 



starting from noninteracting spin chains as 

L{cr, r I A) = F{a, r) {L,,,{a \ A) ® L,,,(r | A)) F{a, r). (9.16) 
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Here Lxxz{,(^ \ A) and Lxxz{t \ A) are related to the independent spin-| operators a, r 
and are equivalent to the Lax operator ( [4 .81) of the XX Z model with Hamiltonian 
(|ti.8|). The twisting operator coupling the two spins is given by 

r) = e'^("'®"'-"'^"') (9.17) 

with the needed property -F12-F21 = /. That makes the coupled system quantum 
integrable with the Hamiltonian 

N 

Hcouple ^ ^ -Pn n+1 

n 

N _ 

= Y.H:L{nn + l) + H::Mn + l). (9.18) 

n 

The interacting parts read as 

H7xzin n + l) = Aalal^^ + /(r^, r^+i)(T+cT-+i + /"^(r^, t^+i)(t~(t^+i 
H::An n + l) = Ar^r^^, + f-\al, a^_,,)r+r-+i + f{al cr^^i)r-r++i (9.19) 

where /(t^,t-^+i) = e'^^'^^+^^+i^ and similarly /(o"^, 0-^+1) = e*^^'^"+'^"+i^ are 
responsible for cubic and quartic spin interactions. 

Recall that the Hubbard model may be given by the Lax operator (|6.27|) involving 



two different spin-i operators. Therefore the same twist ( p. 171 ) can also be applied 



for constructing an extension of the Hubbard model. After transforming to the 
fermion operators C(±)j by (|6.23 , 6.24 ) one can write down the Hamiltonian of such 



an integrable interacting fermion model as 

N 
i 

+ iF^)^,l+lcl)^C(-)^+l + (F+),,,+ic|^),C(_),_i) + U n(+),n(_), (9.20) 

where = e^*^("(±)'+"(±)'+i"^\ 

Note that in the hopping part of the Hamiltonian, one type of fermion is influenced 
by the concentration of the other. Exact eigenvalue solution of these novel integrable 
models through coordinate Bethe ansatz similar to the Hubbard model should be an 



important problem |79 



9.5 Diffusion model as integrable system 

In recent years an important connection between integrable quantum spin chains 



and some one dimensional reaction-diffusion processes has been established |76]. The 
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reaction-difFusion models in the simplest case is described by the probability dis- 
tribution Pn{ti,T2, . . . ,tn I t) of finding certain configuration of molecules and 
vacancies at time t in an open chain of N sites. The value = 1 denotes occupied 
(molecule), while = an empty (vacancy) site k. The master equation represent- 
ing the evolution of probability distribution \ P >= —H \ P > is equivalent to 
an imaginary-time Schrodinger equation, with if as a quantum Hamiltonian. 

If we consider the asymmetric diffusion process involving two-body interactions 
only, then the Hamiltonian may be given by 

N-l 

h=y: {Hu,^iT,lErEii, + {h^)-Er + {h^r.Ei- (9.21) 

k=l 

where due to the existence of the stationary state: < | if = in the matrix 
representation of H the sum of the elements of each column must vanish. Due to 
this crucial property, for the process of partially asymmetric diffusion, we may have 
(iffcfc+i)io — {Hkk+i)oi = —p denoting probability of diffusion to the left and 
to the right, respectively along with (iifc^+i)!]]; = q, {Hkk+i)io = P to ensure the 
required stationary state property. The boundary contributions hi and hjs; describe 
injection (extraction) of particles with rates a and 6 (7 and (3) at boundary sites 
1 and A^, respectively. Therefore (/ii)o = —7 = —{hi)\, (/^i)? = —a = — (/ii)o 
and (/iAf)o = ~P = ~{hM)\, (/^Af)? = —S = — (/iiv)o • After performing an unitary 
transformation, which does not change the average of the observables: < X >= 
J2{t} ^{{'^})Ps{{'''}) , we may rewrite ( |9.21D as the integrable spin chain model, where 



along with the Uq{su{2)) invariant spin Hamiltonian (|6.7| ) additional boundary terms 
appear as 

N _ _ I 

H = -VwE^>n+l + ^n^n+l + T(Q + Q ')(^n^n+l-I) 
n=l ^ 

- ^(g-g-^)(a|-(Tf_,i)+5i + B^ (9.22) 

with Q = y^, where the boundary contributions i?i, contain nondiagonal matri- 
ces cr^ as 

Si = ((a - 7)^? - 2aaf - 27a+ + (a + 7) 

Bn = -^((5-/3)(T^-25g^-V^-2/5gi-^a+ + (/3 + 5)) (9.23) 

The physical consequences of this Hamiltonian related to the diffusion process 
have been analyzed in [131]. However, though the spin model (|9.22 ) can be shown to 
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be integrable ||152|| , due to the lack of the reference state its exact eigenvalue solution 
through Bethe ansatz remains as an important unsolved problem. 



9.6 Spin Calogero-Sutherland model 

In recent years due to relevance in the high superconductivity, interest towards 
models with long-range interactions has been enormously increased. An interesting 
class of integrable models of this kind is the quantum Calogero-Sutherland (CS) 
model and its spin extensions. The Hamiltonian of the spin CS model may be given 



by aniH 



N 

Hcs = Y.P] + '^ E {a^ -aP,k)V{xj-Xk) (9.24) 

j=l l>j<k>N 

with [xj,pk] = iSjk , where the potential V{xj — Xk) = (-^._^^.^)2 for nonperiodic and 
sm^{x -Xk) periodic model. Pjk is the permutation operator responsible for 

exchanging the spin states of the j-th and the k-th particles. In the absence of the 
operator Pij, (|9.24|) turns into the original CS model without spin. 



The spin CS model exhibits many fascinating features, namely its conserved quan- 
tities including the Hamiltonian exhibit Yangian symmetry, the eigenvalue problem 
can be solved exactly using Dunkl operators, the ground state is a solution of the 
Knizhnik-Zamolodchikov equation, the system can be viewed as the free anyonic gas 
related to the notion of fractional statistics etc. ||132|| . The quantum integrability of 



the spin CS model may be proved through an alternative set of operators L and M, 
much in common with the classical case by assuming [Hcs,L] = [L,M] , where L 
is replaced by [Hcs,L] in the quantum case ||133|| . Explicitly, 

Ljk = 6jkPj + ia{l - 6jk)Pjkf{xj - Xk), 

Mjk = 2a{SjkT.Pj''V{xj-Xk)-il-5,k)PjkVixj-Xk)) (9.25) 

with f{xj — Xk) = (^x -xk) nonperiodic and f{xj — Xk) = acot{xj — Xk) 

for periodic model. The conserved set of operators can be constructed as In = 
J2j,kiL^)jk,n > 1 or In = Yjj,k{X°'L'^)jk-,n > 0, with X" being the generators of 
su{n) algebra. Using the crucial property J2j ^jk = Hk^jk = 0, one can show that 
both /„ and commute with the Hamiltonian Hcs-, with I2 = H^s- These conserved 
quantities commute also with the generators of the Yangian algebra. 

Remarkably, at a ^ 00 the Hamiltonian of the CS model ( |9.24D for the periodic 
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case reduces to the well known Haldane-Shastry model ||151|| 



Hhs = E . 2,^'' V- (9-26) 

Jtk {Xj - Xk) 

This discretized long-range interacting spin chain like model seems to be less well 
understood and its Lax operator description difficult to find. 

Using the twisting transformation ( [4.28| ) one can extend the spin CS model in the 



line of to have the deformed Yangian symmetry. The Hamiltonian of such models 
can be given in the form (|9.24|) with the Pj^ operator generalized to QjkPjk, where 
Qi2 = -^12^12-^21^ with the twisting operators Fu and = J2ai3 eaeaa^eaa+eaa^e/?/? • 
The parameters ea = ± can be related to the parity of te particles. Another form of 
generalization can be made through the g-deformation leading to a relativistic spin 
CS model ||134|| . A QYBE type formulation of this model involving i?, L-matrices 



also attempted in |7^, though the related i?-matrix found was a dynamical one. 

Possible description of the spin CS model as a nonultralocal system and the 
formulation of its integrability through braided QYBE, introducing Z and i?-matrices 
following the scheme of this review, would be a challenging problem. 

9.7 Finite-size corrections in integrable systems and related 
CFT 

Statistical systems at critical points should possess conformal invariance. Using this 
fact it is possible to extract the important information about the conformal quantum 
field theory in the scaling limit of the integrable lattice models. Interestingly, from 
the finite size correction of the Bethe ansatz solutions, one can determine the 
CFT characteristics like the central charge and the conformal dimensions. 



Recall that ||135|| the central charge c appears in the central extension of the 
Virasoro algebra 

[L„, Lm] = {n- m)Ln+m + Y2 ~ n)5n+mfl, (9.27) 

the algebra generating the conformal transformations in two dimensions. A is the 
anomalous dimension of the primary fields. For the minimal models M.ppi in CFT, 
the central charge is given by 



6(j9 — p') 



pp 
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with positive integers and finite number of primary fields. For unitary tlieories 
one gets tlie restriction p — p' = 1 yielding 

^ (9.28) 



u u - 1 



with integer p = u > 2 leading to c < 1. 

One may analyze the finite size effect of the Bethe ansatz solutions like ( |7.19| ) 



corresponding to the six-vertex model, but with a seam given by k. Considering the 
coupling parameter q = e*"^, one obtains at the large limit the expression 

^0 = A^/oo - — -c + 0( At) 

for the ground state energy and 

2tt - 1 2-jT ~ 1 

E^-Eo = -{A + A)+0{—) P^-Po = -{A-A) + 0{—) 

for the excited states. Here A, A are conformal weights of unitary minimal models 
and c = 1 — ^j^^jiYJ' = 2, 3, ... is the central charge of the corresponding conformal 
field theory. 

The intriguing relation between the quantum mKdV and the spin-| XX Z chain 
given by equations like ( p.68| , |8l69D frames naturally an worthy problem of finding the 
conformal properties of the quantum mKdV model by suitable choice of k and u. 

9.8 Integrable structure of CFT through quantum KdV 

To capture the integrable structure of CFT, an alternative approach through mass- 
less S-matrix description and based on the quantum KdV model has been proposed 
recently [|75| ,. 



Note that the conformal symmetry of CFT is generated by its energy-momentum 
tensor T(u) = — ^ + J2°^oo L-n^^"'^ , with L„ satisfying the Virasoro algebra (|9.27|) . 



The operators hk-i = ^ jo"" duT2k{u) with T2k{u), depending on various powers 
and derivatives of T{u) represents an infinite set of commuting integrals of motion. 
The idea is to solve their simultaneous diagonalization problem, much in common 
to the QISM for the integrable theory. Remarkably, this is equivalent to solving the 



quantum KdV problem ||136|| , since at the classical limit the field T{u) = —^U{u 



with U{u + 2tt) = U (u) reduces the commutators of T{u) to 

{U{u), U{v)} = 2{U{u) + U{v))6'{u -v) + 6'"{u - v) (9.29) 
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which is the well known Poisson structure of the KdV. 

In the quantum case, defining a Miura transformation through a mKdV like field 
4>{u), the monodromy matrix of the model can be expressed as |73| 



M,(A) 



(9.30) 



with proper normal ordering : :. Here E,F,H are generating elements of Uq{sl{2)) 
and TTj is its 2j + 1 dimensional representation. Interestingly, the operators Mj{X) 
satisfy now the QYBE like (|2.11|) with the trigonometric i?(A)-matrix. The transfer 
matrix Tj(A) = tr.,;.{e^'^^^ Mj{\)) commutes for different parameters A as well as with 
the operators J2A:, ( which are related to the expansion coefficients of Inrj(A)). For 
the values of the central charge 

c=l-3^-— f, n = l,2,3,... (9.31) 
in + 3 

the eigenvalues of Tj{\) was found to satisfy a closed system of functional equations, 
which is sufficient to derive the spectrum. Such equations are obtained from the 
operator relations 

r,(gh)r,(g-h) = 1 + r^._i (A)r^.^i .(A) (9.32) 

For the ground-state eigenvalues the functional equations are equivalent to those of 
the thermodynamic Bethe ansatz equations for the massless ^'-matrix theory associ- 
ated with the minimal CFT M.2,2k+i- 

Since the quantum mKdV model plays determining role in this approach, appli- 
cation of its exact quantum Lax operator as presented in ( p.58| ) might be promising. 



Moreover, the finite-size corrections of its Bethe ansatz solution (|8.68|) might also be 
fruitful in deriving the central charge like ( p.31|) . 



9.9 Other problems 

There are important classically integrable nonultralocal models like complex sine- 
Gordon ||62|, nonlinear cr-model |6^, KdV model etc., whose quantum generalizations 



remain a major open problem. It is desirable to apply the scheme for quantum 
nonultralocal models presented here or its generalizations to these models, 
i) KdV model 

The second Hamiltonian Poisson bracket structure of the KdV field is given by 
(|9.29|) and due to the appearance of higher derivative of (5-function the generalization 



of the scheme should include more number of braiding matrices Z |92|, for incorpo- 
rating the nonultralocality beyond nearest neighbors. 
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ii) Complex sine- Gordon model 

The complex SG model is a relativistic model given by the Lagrangian density 
£ = ^ — m'^i/jij?! , expressed through the complex scalar field i/j with canon- 

ical bracket like {TT^{x),ip{y)} = 6{x — y) and similarly for its conjugate field -0. 
The related Lax operator is given in the form p2| 



— if %Tfl — \ 

+2i{i-ijij)-^ [dQ^ij + d^ij - —{i-ijij)ij\a+ ^ c.c. (9.33) 

where c.c stands for the complex conjugate. Due to the dxil^ term present in the 
Lax operator, the model exhibit nonultralocal property, 
iii) Nonlinear a-model 

It is a relativistic model given by the action S = ^ J d'^x^d^UaY, = 1, a = 
1,2, . . . , N. The Lax operator of the model may be given by [^] 

L = -^^^(AAi + Ao), 4 = 2[n x d^n], (9.34) 
where the Poisson brackets among the na{x) are 

{na{x), nb{y)} = 0, {na{x), TTb{y)} = {6ab - nanb){x)S{x - y), 

{-Kaix), T^b{y)} = (TTa^fc " nbna){x)5{x - y) , T^a = ^O^a- (9.35) 



The presence of the derivative in x in the Lax operator ( |9.34| ) through the term Ai 



makes the model nonultralocal and quantum generalization difficult. 



fig.9.1 



Figure 9.1 Ladder symmetric models 
a) Quantum integrable theoretical model l \9.1^ ) 



b) Ladder structure of Cu2 O3 layers contained in the high superconducting crystals 

m 

10 Concluding remarks 

Recent advancement in the theory of quantum integrable systems, to which the 
present review is dedicated, have shown a remarkable interpenetration between for- 
mal algebraic structures appearing in pure mathematics and the integrable theory 
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of physically interesting quantum models. The quantized algebra and the quan- 
tized braided algebra are found to be the important structures underlying integrable 
systems, where their crucial Hopf algebra property plays a significant role in the 
transition from local to global form of the quantum Yang-Baxter equation. This in 
turn leads to the integrability of both ultralocal and nonultralocal systems. A major 
emphasis is given here to highlight this interplay between abstract algebras in one 
hand and the construction of concrete quantum models on the other. We have tried 
to maintain an elementary level stressing mainly on physical insights and avoiding 
mathematical rigor. The basic aim is to understand the algebraic aspect inherent to 
the quantum integrable systems and to see how these formal structures can be fruit- 
fully applied for constructing key objects in the integrable theory for both ultralocal 
and nonultralocal systems. 

A scheme for systematic generation of the representative Lax operators in ultralo- 
cal models is presented. The approach in a sense is reverse to the standard QISM and 
instead of starting from concrete models or known Lax operators and subsequently 
deriving quantum i?-matrices, it constructs the i?-matrix and the Lax operator by 
using abstract objects like universal 7^-matrix related to the quantum algebra. In- 
troducing spectral parameter through an algebraic procedure one constructs for the 
Hecke algebra class the Lax operator of an ancestor model. Through different realiza- 
tions of the quantized algebra or its limiting Yangian, one builds up known as well as 
new integrable models. A wide range of models as descendant models are thus gener- 
ated from the same algebra and sharing the same quantum i?-matrix. The integrable 
systems are therefore classified into various classes represented by their i?-matrices 
and ancestor Lax operators (see fig 8.1). Due to this relationship at the algebraic 
level, the diverse models like sine-Gordon, XX Z spin chain, derivative Schrodinger 
equation, Liouville model, Ablowitz-Ladik model etc. show remarkably common fea- 
tures at the action-angle level, reflected in their eigenvalue solution through algebraic 
Bethe ansatz. The intriguing question of finding the criteria for the nonlinearity that 
makes a system integrable, seems also to be answered by the specific form of realiza- 
tion of the algebra related to the ancestor models. A general realization in bosonic 
field ipy'^''^ of such ancestor model belonging to the trigonometric class may be given 
by the L operator found in |^ , while the L operator of ||140|| is a similar realization 
in the rational case. 

Lattice regularization of the field models required for tackling ultraviolate di- 
vergences often becomes a nontrivial task, even for deriving integrable models with 
required approximation. The present construction on the other hand yields, exploit- 
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ing the underlying algebra determined by the QYBE, integrable lattice models exact 
in all orders of the lattice constant A. One of the drawbacks of this scheme however 
is that, for wider class of models with nonsymmetric Lax operators belonging to the 
extended trigonometric Sklyanin algebra, the universal 7?.-matrix solution could not 
be found. This should therefore be an open problem to be solved. Similarly general- 
ization to include other algebras like rotational, projective as well as supersymmetric 
algebras would be important directions. 

Due to the recent development in formulating braided extensions of the QYBE, 
the integrable theory for quantum nonultralocal models, which include important 
models like quantum mKdV, nonabelian Toda chain, WZWN model, anyonic models 
etc. can now be presented in a rather unifying way (sec fig. 9.1). The formulation 
is almost similar to the well established theory for ultralocal models, though many 
aspects in this regard need further refinement. Such derivation of known examples 
from a general scheme by particular choices of braiding matrices Z and Z seems to 
suggest also a possibility for their nontrivial extensions. For example, choosing Z ^ 1 
in the reduction for the WZWN model, one is likely to obtain generalizations of the 
known current algebra. 

A distinguishing feature of this review perhaps is its parallel treatment for both 
ultralocal and nonultralocal quantum models emphasizing their algebraic aspects, 
integrability, and interrelations. An extended list of concrete physical models of both 
these classes are provided with details of their constructions and solutions. After 
constructing systematically the representative Lax operators from more fundamental 
level, the next task is to derive the Hamiltonian of the concrete models starting 
from the Lax operators and then to solve the eigenvalues of this Hamiltonian along 
with other conserved quantities. In constructing Hamiltonian we have focused on 
different methods adopted for fundamental and nonfundamental models in concrete 
examples. Likewise, the specialization required for the Bethe ansatz solution for 
handling integrable field models, e.g. SG or DNLS models from that of the lattice 
models like spin chains are highlighted. Similarly, in presenting Bethe ansatz method 
for solving exactly the eigenvalue problems, both its algebraic as well as coordinate 
formulations are described on the same set of models with comparison. This is done 
to bring out the deep relationship between these two fascinating methods and also 
for better understanding them. 

The depth with which the ultralocal systems can be treated is still not possible 
with nonultralocal models in the quantum case. For example we are not able to 
generate yet systematically the Lax operators for nonultralocal models from more 
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fundamental objects like universal R-matrices or ancestor models. Full theory should 
also be able to describe all models with arbitrary braiding and in particular should 
cover well known models like quantum KdV, chiral models, nonlinear cr-models, com- 
plex sine-Gordon model, Calogero-Sutherland models etc. 

Since the subject is evergrowing, many interesting recent developments are tak- 
ing place with yet new problems opening up. We have provided in the last section 
a brief account of some of such recent progress, which can be understood based on 
the material presented in the review. We also indicate in this section various open 
problems and desirable directions stemming from the recent results. 
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